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Summary 

The  whole-body  absorption  of  electromagnetic  waves  by  biological 
bodies  is  strongly  dependent  on  the  orientation  of  the  electric  field 
(E)  relative  to  the  longest  dimension  (L)  of  the  body.  The  highest 
rate  of  energy  deposition  occurs  for  eH  L  orientation  for  frequencies 
(65-70  MHz  for  an  adult  human  being)  such  that  the  major  length  is  ap¬ 
proximately  0.36  to  0.4  times  the  free-space  wavelength.  Peaks  of 
whole-body  absorption  for  the  other  two  configurations  (major  length 
oriented  along  the  direction  of  propagation  (k|l  L)  or  along  the  vector 
of  the  magnetic  field  (H||  L))  occur  for  X  on  the  order  of  4nb,  where  2TTb 
is  the  weighted  average  circumference  of  the  animals. 

Numerical  calculations  for  a  realistic  block  model  of  man  have 
shown  a  fine  structure  to  whole-body  absorption  at  frequencies  higher 
than  the  whole-body  resonant  frequency.  Minor  peaks  in  the  supraresonance 
region  are  ascribed  to  maxima  of  energy  deposition  in  the  various  body 
parts  such  as  the  arm  and  the  head.  For  the  e|1  L,  supraresonant  region, 
an  average  1/f  dependence  of  the  whole-body  absorbed  dose  is  obtained 
theoretically  and  experimentally  to  frequencies  on  the  order  of  4-8 
times  the  resonance  frequency. 

The  effect  of  highly  conducting  (metallic  sheet)  ground  and 

reflecting  surfaces  on  the  electromagnetic  absorption  and  its  distribu- 
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tion  has  been  studied.  For  a  standing  man  model  (e||  L)  with  feet  in 


conductive  contact  with  a  perfect  ground,  the  new  resonant  frequency 
is  roughly  one-half  that  for  free-space  Irradiation.  At  this  lower 
resonance  frequency,  the  mass  normalized  (specific)  absorption  rate  (SAR) 


is  about  twice  that  at  peak  absorption  frequency  for  free-space  ir¬ 
radiation,  and  the  ground  effects  are  lost  even  for  a  3-4"  separation 
from  the  ground. 

Maximum  deposition  of  electromagnetic  energy  in  the  head  region 
occurs  for  head-to-toe  propagation  at  frequencies  such  that  the  head 
diameter  is  approximately  3.0  times  the  physical  cross  section.  This 
value  exceeds  the  ratio  of  nearly  1  reported  previously  for  an  isolated 
sphere  model  of  the  head. 

A  multilayered  planar  model  is  used  to  examine  the  dependence 
of  v/hole-body  power  absorption  on  the  configuration  of  surface  layers, 
e.g.,  skin,  fat,  bone,  muscle,  which  normally  occur  in  biological 
bodies.  It  is  found  that  the  layering  resonance  for  three-dimensional 
bodies  (as  opposed  to  geometrical  resonance)  can  be  predicted  quite 
accurately  by  a  planar  model.  Calculations  for  a  multilayered  prolate 
spheroidal  model  of  man  predict  a  whole-body  layering  resonance  at  1.8 
GHz  with  a  power  absorption  34  percent  greater  than  that  predicted  by 
a  homogeneous  model. 

Behavioral  experiments  with  laboratory  rats  confirm  the  high¬ 


lights  of  the  above  data. 


FOREWORD 


In  conducting  the  research  described  in  this  report,  the 
investigators  adhered  to  the  "Guide  for  Laboratory  Animal  Facilities 
and  Care"  as  promulgated  by  the  Committee  on  the  Guide  for  Laboratory 
Animal  Resources,  National  Academy  of  Sciences  -  National  Research 
Council . 
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ELECTROMAGNETIC  ENERGY  ABSORPTION  AND  ITS  DISTRIBUTION  FOR  MAN 
AND  ANIMALS  AT  DIFFERENT  FREQUENCIES  UNDER  VARIOUS  CONDITIONS 


Objectives 

The  objectives  of  the  project  were: 

1.  To  quantify  the  electromagnetic  power  absorption  and  its 
distribution  for  man  when  subjected  to  radiation  at  dif¬ 
ferent  frequencies  and  under  various  exposure  conditions. 

2.  To  verify  the  important  findings  of  part  (1)  by  exposure  of 
living  organisms. 

The  experiments  for  man  were  performed  with  reduced  proportion¬ 
ately-scaled  models  filled  with  0.9  percent  saline  solutions  and  with 
biological-phantom  materials  developed  to  simulate  electromagnetic 
properties  of  human  tissue.  The  animal  experiments  consisted  of  micro- 
wave-induced  convulsions  and  work  stoppage  in  rats. 

The  purpose  of  the  project  was  to  develop  an  understanding  that 
would  lead  to  projections  for  humans. 

Highlights  of  the  Work  Done  on  the  Project 

* 

A.  Electromagnetic  Energy  Deposition  in  Man  and  Animals 
Free-Space  Irradiation  Condition 

The  condition  that  has  been  studied  the  most  to  date  is  that  of 


A  paper  based  on  these  highlights  will  be  submitted  for  publication  in 
the  Proceedings  of  the  IEEE  (January  1980)  Special  Issue  on  Biologi¬ 
cal  and  Ecological  Effects  and  Medical  Applications  of  Electromag¬ 
netic  Energy. 
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free-space  irradiation  of  single  animals.  The  whole-body  absorption 
of  electromagnetic  waves  by  biological  bodies  is  strongly  dependent  on 
the  orientation  of  the  electric  field  (E)  relative  to  the  longest  dimen¬ 
sion  (L)  of  the  body.  The  highest  rate^  of  energy  deposition  occurs  for 
E|i  L  orientation  (see  Appendix  A),  for  frequencies  such  that  the  major 
length  is  approximately  0.36  to  0.4  times  the  free-space  wavelengtn  (A) 
of  radiation.  Peaks  of  whole-body  absorption  for  the  other  two  configura¬ 
tions  (major  length  oriented  along  the  direction  of  propagation  (kl|  L) 
or  along  the  vector  of  the  magnetic  field  (h1|  L)  have  also  been  reported^ 
for  X  on  the  order  of  4Tib,  where  2'irb  is  the  weighted  average  circumfer¬ 
ence  of  the  animals. 

Using  prolate  spheroidal  and  ellipsoidal  equivalents  of  biologi¬ 
cal  bodies,  theoretical  calculations  have  recently  been  given  in  a 
dosimetry  handbook^  for  frequencies  up  to  and  slightly  beyond  the 
resonant  region  for  the  aforementioned  polarizations  e1|  L,  k||  L,  and 
h||  L.  Numerical  calculations  for  a  realistic  model"*  of  man  (see  Ap¬ 
pendix  B)  have  shown  a  fine  structure  to  whole-body  absorption  at  fre¬ 
quencies  higher  than  the  whole-body  resonant  frequency.  Minor  peaks  in 
the  supraresonance  region  are  ascribed  to  maxima  of  energy  deposition 
in  the  various  body  parts  such  as  the  arm  and  the  head^  (see  Appendix  C 
for  details) . 

For  the  supraresonant  region,  the  E||  L  orientation  has  been 
studied  most  extensively.  An  average  1/f  dependence  of  the  whole-body 
absorbed  dose  is  experimentally  observed^  to  frequencies  f  on  the  order 
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of  1.6  S  times  the  resonant  frequency  f  ,  where  S  is  the  relative 

res  r’  res 

absorption  cross  section  at  the  resonant  frequency. 

Empirical  equations  for  specific  absorption  rate  (SAR)  have  been 
derived^  for  e||  L  orientation  and  found  to  be  fairly  accurate  for  six 
animal  species^  from  25-g  mice  to  2250-g  rabbits. 

These  equations  are  given  in  Table  1. 

The  coefficients  In  these  empirical  equations  for  man  model  are 
approximatley  63  percent  of  the  live  animal  coefficients.  Further 
research  is  needed  to  resolve  this  difference. 

Electromagnetic  Absorption  in  the  Presence 

of  Nearby  Ground  and  Reflecting  Surfaces 

Only  highly  conducting  (metallic  sheet)  ground  and  reflecting 
surfaces®’^  of  infinite  extent  have  been  studied  to  date  (see  Appendix 
D) .  For  a  standing  man  model  with  feet  in  conductive  contact  with  a 
perfect  ground,  there  is  a  drastic  alteration  in  SAR  as  a  function  of 
frequency.  For  e||  L  orientation  the  new  resonant  frequency  is  roughly 
one  half  that  given  by  Eq.  1.  At  this  lower  resonance  frequency,  the 
SAR  is  about  twice  that  at  peak  absorption  frequency  for  free-space 
Irradiation. 

The  nature  of  the  ground  effects  on  SAR  (for  e[|  L  orientation) 
is  such  that  even  a  small  separation^”^  from  ground  (to  break  conductive 
contact)  is  sufficient  to  eliminate  much  of  the  ground  effect  (see  Ap¬ 
pendices  D  and  E) .  For  separations  from  ground  more  than  3-4",  the 
total  energy  deposition  and  its  distribution  are  identical  to  those  for 


Table  1.  SAR  for  man  models  for  conditions  of  free-space  irradiation 


e||  L  Orientation 


Peak  absorption  or  resonant  frequency: 


f  =  11.4/L  GHz 
r  cm 

For  subresonant  region  —  0.5  f^  <  f  <  f^: 


(1) 


SAR  in  mW/g  for  1  mW/cm 
incident  plane  waves 


0.522 


L 


2 


_ cm _ 

t<?eight  in  g 


(2) 


For  supraresonant  frequency  region  —  f^  <  f  <  1.6  f^: 


SAR  in  mW/g  for  1  mW/cm^  _  5.95  ^cm 
Incident  plane  wave  fields  ^®^8ht  in  g 


where  is  the  long  dimension  of  the  body  in  centimeters  and 


S 

res 


0.481  y - - 

\ weight  in  g 


(4) 
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free-space  irradiation  conditions.  Even  for  a  man  model  in  conductive 
contact  with  a  perfect  ground,  the  energy  deposition  in  the  supraresonance 
region  (f  >  2-3  f^)  is  comparable^  to  that  for  conditions  of  free-space 
irradiation. 

Other  orientations  and/or  finite  conductivity  ground  effects  on 
SAR  have  not  been  studied  to  date. 

For  highly  conducting  (metallic  sheet)  reflecting  surfaces  of 
flat  and  90°-corner  types,  enhancements^  in  SAR  by  factors  as  large  as 
27  have  been  observed  for  e||  L  orientation.  Most  of  the  work  to  date 
has  concentrated  on  frequencies  close  to  the  resonance  region.  The  ob¬ 
served  enhancement  factors  are  explained^  on  the  basis  of  antenna 
theory. Indeed,  for  incident  plane  waves  for  e||  L  orientation,  most 
of  the  observed  results  are  as  though  the  target  acted  like  a  pick-up 
half-wave  dipole  with  reflecting  surfaces  in  close  vicinity. 

Finite  conductivity,  finite  size  reflecting  surfaces  and  other 
orientations  have  not  been  considered  to  date.  Also  the  results  for 
frequencies  higher  than  about  8.5  times  the  resonant  frequency  (550  MHz 
for  man)  have  not  been  obtained  even  for  highly  conducting  reflectors. 

Enclosed  structures  such  as  rooms,  etc.,  may  act  as  lossy 
resonators  with  electromagnetic  fields  coupling  in  from  windows.  If 
such  structures  were  to  have  highly  reflecting  walls,  field  enhance¬ 
ments  by  two  to  three  orders  of  magnitude  may  indeed  be  possible.  How¬ 
ever,  since  typically  encountered  walls  are  not  very  reflecting,  it  is 
expected  that  field  enhancements  by  more  than  a  factor  of  5-10  may  not 
occur  realistically.  Research  into  these  structures  is  needed  in  order 
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to  describe  the  nature  of  field  enhancements. 

Head  Resonance 

We  have  recently  identified  a  frequency  region  for  the  highest 
rate  of  energy  deposition  in  the  head.  The  head  resonance^* ® ^ ^  (see 
Appendices  C,  E,  and  F)  occurs  at  frequencies  such  that  the  head 
diameter  is  approximately  one  quarter  of  the  free-space  wavelength. 

For  the  intact  (adult)  human  head,  the  resonance  frequency  is  estimated 
to  be  on  the  order  of  35O-A0O  MHz.  At  head  resonance,  the  absorption 
cross  section  for  the  head  region  is  approximately  3.0  times  the  physical 
cross  section  with  a  volume  averaged  SAR  that  is  about  3.3  times  the 
whole-body  averaged  SAR.  Both  values  greatly  exceed  numbers  reported 
earlier^^*^^  for  spherical  models  of  the  isolated  human  head.  Also 
numerical  calculations^^  (Appendix  F)  using  144  cubical  cells  of  various 
sizes  to  fit  the  shape  of  the  human  head  (340  cells  for  the  total  body) 
give  local  SARs  at  hot  spots  that  are  about  5  times  the  head-average 
values. 

We  feel  that  the  phenomenon  of  head  resonance  may  be  important 
in  the  study  of  behavioral  effects,  blood-brain  barrier  permeability, 
cataractogenesis ,  and  microwave  bioeffects. 

Multianimal  Effects^ 

It  has  been  shown  (see  Appendix  C  for  details)  that  for  resonant 
biological  bodies  close  to  one  another,  the  antenna  theory  may  be  used 
to  predict  the  modification  In  SAR  relative  to  free-space  values.  For 
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two  resonant  targets  separated  by  0.65  to  0.7  A,  the  highest  SAR,  150 
percent  of  the  free-space  value  has  been  experimentally  observed  for 
anesthetized  rats  and  models  of  man  for  e||  L  orientation  for  frontally 
(broadside)  incident  plane  waves.  For  three  animals  in  a  row  with  an 
interanimal  spacing  of  0.65  X,  the  central  animal  SAR  is  observed  to  be 
roughly  two  times,  while  the  two  end  animals  receive  an  SAR  that  is  ap¬ 
proximately  1.5  times  that  for  an  isolated  animal. 

Full  implications  of  the  multibody  effects  on  SAR  are  not 
completely  understood,  even  though  preliminary  experimental  data^  show 
that  similar  enhancements  may  also  occur  for  subresonance  and  supra- 
resonance  regions.  Furthermore,  other  orientations,  irregular  spacings, 
and  nonf ree-space  exposure  conditions  have  not  been  considered  to  date. 

B.  Other  Highlights  of  the  Work  Done  on  the  Project 

Electromagnetic  Absorption  in  a 

Multilayered  Model  of  Man 

A  multilayered  planar  model  is  used  (see  Appendix  G)  to  examine 
the  dependence  of  whole-body  power  absorption  on  the  configuration  of 
surface  layers,  e.g.,  skin,  fat,  muscle,  which  normally  occur  in  biologi¬ 
cal  bodies.  It  is  found  that  the  layering  resonance  for  three-dimensional 
bodies  (as  opposed  to  the  geometrical  resonance)  can  be  predicted  quite 
accurately  by  a  planar  model.  Calculations  for  a  multilayered  prolate 
spheroidal  model  of  man  predict  a  whole-body  layering  resonance  at  1.8 
GHz  with  a  power  absorption  34  percent  greater  than  that  predicted  by 
a  homogeneous  model. 
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Behavioral  and  Thermal  Effects  of  Microwave 
Irradiation  at  Resonant  and  Nonresonant  Wavelengths 

The  details  of  these  studies  are  given  in  Appendix  H. 

Behavioral  and  thermal  effects  of  radiating  an  animal  with  dif¬ 
fering  wavelengths  of  microwave  energy  at  the  same  power  density  were 
investigated  in  the  first  of  two  studies.  Five  Long-Evans  rats  were 
trained  to  perform  a  lever-pressing  task  and  were  rewarded  with  food  on 
a  variable  interval  schedule  of  reinforcement.  Rats  were  individually 

exposed  in  random  order  to  400-,  500-,  600-,  and  700-MHz  CW  radiation 

2 

at  a  power  density  of  20  mW/cm  with  the  long  axis  of  the  rat's  body 
parallel  to  the  vector  of  the  electric  field.  Radiation  at  all  wave¬ 
lengths  produced  rises  of  body  temperature  and  stoppage  of  lever  press¬ 
ing.  The  averaged  rise  in  body  temperature  was  greatest  and  work  stop¬ 
page  was  most  rapid  during  exposures  at  600  MHz.  In  the  second  study, 

six  rats  were  exposed  in  random  order  to  600-MHz  CW  radiation  at  power 

2 

densities  of  5,  7.5,  10,  and  20  mW/cm  while  performing  the  same 

2 

behavioral  task.  Exposures  at  10  and  20  raW/cm  resulted  in  work  stop- 

2 

page,  while  exposures  at  5  and  7.5  mW/cm  did  not.  In  addition,  three 

of  the  rats  were  subsequently  exposed  while  responding  to  600-MHz 

pulsed  radiation  (1000  pps,  3-  or  30-ps  pulse  durations  at  a  peak  power 

2  2 

density  of  170  mW/cm  (averaged  0.51  and  5.1  mW/cm  ).  No  work  stoppage 
occurred  to  pulsed  radiation.  Taken  in  sum,  the  data  show  that  the 
mature  Long-Evans  rat  is  resonant  at  a  frequency  near  600  MHz  while  work 
stoppage  during  short-term  exposures  to  600-MHz  radiation  occurs  at  a 


power  density  between  7.5  and  10  mW/cra 
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Radiators  for  Microwave  Biological  Effects  Research  — 

Waveguide  Slot  Array  with  Constant  Radiation  Intensity 

Microwave  biological  studies  with  large  targets  or  multiple 
targets  require  a  large  working  area  that  provides  near-uniform  exposures 
to  energy.  At  high  microwave  frequencies  where  high  power  sources  are 
not  readily  available,  it  is  desirable  to  use  a  radiator  that  provides 
a  constant  intensity  beam  so  that  the  working  area  may  receive  almost 
all  of  the  radiated  energy. 

Several  synthesis  procedures  have  been  used  for  design  of  a 
longitudinal  shunt  slot  array  having  a  pattern  with  nearly  constant 
radiation  intensity  over  the  beam  width  and  zero  intensity  outside.  A 
beam  width  of  15.5°  was  chosen  as  being  small  enough  for  a  radiated  beam 
to  approximate  a  plane  wave  yet  give  enough  divergence  for  convenient 
chamber  divisions.  Discontinuities  in  the  desired  pattern  force  the 
Fourier-synthesized  pattern  to  have  significant  errors  due  to  overshoot 
and  ripple.  Solutions  developed  on  the  basis  of  FejSr  sums  and  Lanczos 
weighted  least  squares  gave  patterns  with  objectionable  rounding.  A 
new  procedure  based  on  weighted  least  squares  was  developed  for  optimum 
characteristics  and  both  the  new  procedure  and  its  application  to  a 
practical  design  are  described  in  reference  14. 

A  longitudinal-shunt  slot  array  with  reflectors  has  been  designed, 
built,  and  tested.  The  antenna  is  easily  matched  near  the  design  fre¬ 
quencies  of  10.0-10.4  GHz.  The  radiation  pattern  has  an  intensity  within 
±0.5  dB  over  a  total  angle  of  15.6°  in  the  H  plane  and  15.5°  in  the  E 
plane.  Gain  is  19.7  dB  which  is  3.9  dB  above  that  of  a  horn  with  a 
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similar  aperture  area.  By  using  energy-absorbing  material,  it  is  possible 
to  produce  a  pattern  in  which  the  fields  drop  suddenly  from  the  nearly 
constant  values  to  nearly  zero. 

Biological  Phantom  Materials  for  Simulating 

Man  at  Different  Frequencies 

Biological  phantom  materials  were  developed  for  simulating  man 
over  the  frequency  range  13-230  MHz  (including  the  important  resonance 
regions  of  33  and  65  MHz).  These  materials  have  been  used  to  fill  the 
eight  proportionately-scaled  figurines  for  dosimetric  measurements  at 
300,  400,  600,  915,  and  985  MHz.  The  compositions  and  their  measured 
properties  are  given  in  Appendix  A  of  Annual  Report  No.  2. 

Procedures  for  Improving  Convergence  of 

Moment-Method  Solutions  in  Electromagnetics 

The  value  of  numerical  solutions  in  dosimetry  may  be  limited  by 
both  the  accuracy  of  the  model  as  a  representation  of  the  target  animal 
and  the  accuracy  of  the  computations  for  the  model.  Accurate  numerical 
procedures  have  been  developed  for  use  with  crude  models  such  as  the 
sphere,  infinite  cylinder,  ellipsoid,  and  spheroid.  We  have  used 
moment-method  solution  of  the  electric  field  volume  integral  equation 
so  that  an  array  of  cubes  of  different  sizes  and  compositions  may  be 
used  to  better  represent  th.>  shape  of  the  target  animal.  In  order  to 
allow  a  maximum  number  of  cubes  for  increased  detail  in  the  representa¬ 
tion  of  arms,  legs,  etc.,  we  have  used  a  pulse  function  basis  in  the 
moment-method  solution.  Others  have  also  used  pulse  function  solutions 
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for  block  models  in  geophysics  and  in  dosimetry,  but  we  have  addressed 


the  issues  of  both  estimating  and  reducing  the  errors  caused  by  the  use 
of  pulse  functions. 

A  pulse  function  basis  requires  the  assumption  that  the  cells 
are  small  enough  that  the  electric  field  may  be  represented  by  a  constant 
complex  vector  within  each  cell.  Solutions  previously  reported  by 
others  for  a  block  model  of  man  have  a  ratio  of  239:1  for  energy  deposi¬ 
tion  in  one  pair  of  adjacent  cells  at  10  MHz.  We  expect  that  large 
errors  are  present  in  solutions  exhibiting  large  cell-to-cell  variations. 
Our  first  publication  regarding  the  use  of  pulse  functions  is  Appendix  M 
which  establishes  an  upper  bound  on  cell  size.  Such  a  criterion  was  not 
previously  available  for  use  with  moment-method  solutions.  After  defin¬ 
ing  the  source  of  error,  our  next  task  was  to  reduce  the  error.  Ap¬ 
pendix  I  describes  several  novel  procedures  which  we  discovered  for  re¬ 
ducing  the  error  in  solutions  of  2-D  electromagnetics  problems.  We  were 
not  able  to  extend  the  new  procedures  to  3-D  solutions,  but  we  found  a 
different  method  which  is  useful  for  improving  convergence  in  3-D 
problems  which  is  described  in  Appendix  ,1 .  Armed  with  both  our  con¬ 
vergence  criterion  (Reference  15)  and  a  novel  method  of  improving  con¬ 
vergence  (Appendix  J) ,  we  have  proceeded  to  obtain  numerical  solutions 
for  a  realistic  model  of  man  both  in  free-space  and  under  other  condi¬ 
tions  such  as  in  the  presence  of  a  ground  plane  or  flat  or  90°  corner 
reflectors.  Such  applications  are  described  in  other  publications 
referenced  earlier  in  this  report. 
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Conditions  of  Strongest  Electromagnetic  Power 
Deposition  in  Man  and  Animals 
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Abstract — Strongest  power  deposition  for  biological  bodies  is 
found  for  fields  polarized  along  the  longest  dimension  for  frequencies 
such  that  the  major  length  is  about  0.4  times  the  free-space  wave¬ 
length  of  radiation.  Peak  absorption  in  the  presence  of  ground  effects 
is  observed  at  frequencies  about  one-half  as  much  as  for  bodies 
isolated  in  free  space.  At  resonance,  an  effective  absorption  area  3-4 
times  the  shadow  cross  section  is  measured.  Using  biological-phan¬ 
tom  figurines  of  the  human  body,  distribution  of  power  absorption  is 
determined.  Maximum  power  deposition  is  observed  for  the  neck  re¬ 
gion.  Initial  experiments  with  anesthetized  and  dead  rats  have  con¬ 
firmed  this  observation. 

INTUOmXTlOX 

KC’.M’SK  of  the  expanding  uso.-^  of  elootroniagnctic 
(K.M)  radiation,  it  ha.s  bcoome  nm‘.«.'<arv  to  cstah- 
li.sh  tlic  most  hazardous  conditions  of  <‘xpo.surc  to  such 
radiation  for  immans.  In  particular,  it  is  important  to 
know  the  frc(|ucnc,v  regions  of  most  absorption  for  v.arions 
orientations  of  the  bod.v  relative  to  th('  polarization  of 
incident  waves  and  how  the  ab.sorjition  cro.ss  section  varies 
with  frrainency  of  li.M  radiation.  It  is  al.so  imt)orlaut  to 
determine  the  distribution  of  power  absorption  in  the  Ixtdy 
to  a.scertain  locations  of  maximum  di'[)osition  particularly 
close  to  resonant  fnapicncies.  This  knowledge  is  vital  in 
evaluating  and  establishing  radiation  safety  .standards. 

Ua.scd  on  our  [)revious  experiments  [IJ  with  rats  and 
biological-phantom  an<l  saline-fillcfl  prolate  splnaoidal 
bodies,  strongest  whoh'-body  absoridion  was  fo\m<l  for 
K.M  fields  polarizeil  with  electric  field  along  the  long 
dimensions  of  the  bo<Iies.  KM  power  absorption  was 
measured  for  bodies  of  fixed  size  at  various  fnapiencies 
using  a  i)aralli'l-plate  waveguide  as  a  transmission  medium 
of  i)lane  waves  anri  compared  p]  with  data  obtained  from 
fixed-friapiency  free-space  irradiation  exposures  of  bodies 
of  varying  sizes.  'I'lie  salient  features  of  the  experimentally 
f)b.serve<l  results  ((|ualitatively  sketched  in  Kig.  I)  areas 
fc  illows. 

1)  'I’he  £  [I  //  orientation  is  found  to  be  the  iuo.st 
absorbing,  and  |1  K  the  least  absorbing,  with  the  con¬ 
figuration  k  11  L  onlv  slightly  more  ab.sorbitig  than  the 
H\\  /.  orientation.  'I'he  vectors  E,  H,  and  k  are  along  the 


Frequency  (arbiiraTv  units) 


EiK.  1  Tyniciil  RK  absorption  curves  for  animals  and  bodies  of 
prolaU'-sphcroidal  shape  obtaineil  with  a  |)arallel-plalc  wavt>- 
gui<|p.  Symbols  with  arrows  above  are  rppre.senle<i  by  boldface 
italic  in  text . 

electric  and  mtignetic  fields,  and  along  the  diri'ction  of 
propagtifion,  resirectively.  and  I,  is  along  the  major  length 
//  of  the  body.  The.se  exiM'rimental  results  tire  in  (piali- 
tative  agreement  with  ri'cent  theoretical  results  [d]  on 
power  absorption  in  homogeneous  los.sy  prolate'  spheroids, 
(iuantitafive  comptirison  has  not  been  po.s.sible  hecau.se  of 
the  long-wavelength  tkK<0.(>)  approximation  of  the 
theory.  Hy  contrast,  the  ('xperinii'ntal  results  pertain 
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mostly  to  fretpieneies  on  the  order  of  resotiant  and  higher 
fretiut'iieies  (klj>  1.0). 

2)  The  fre(|Ueneies  of  peak  absorption  occur  in  the 
reverse  order.  .Maximum  ab.sorption  for  E  [J  /,  occurs  at 
the  lowest  fretpieney  with  kl.  -  ‘IrL'k  on  the  order  of 
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TABLE  I 

Actital  and  Eqdivalent  Prolate  Spheroidal  Dimensions  por 
Anesthetized  Rats  op  Dipperent  Weights 


U*t|ht  of 
tb*  Rat 

In  CroM 

tat  Langth 

L  froa  Noaa 
to  taaa  of 

Tail 

(• 

- 1 

Naaaurod 

Clrcuafaronca 
of  th«  AOlBol 

llMlBal 

Clrcuafaraaea* 

2«b 

In  CM 

Aapact  tatlo 

L/2b  of  Mul- 
valant  Pcolaia 
Syharoid 

96 

15.5 

11.5 
(body) 
mm.  9.0 
(haad) 

10.80 

4.50 

158 

18.0 

n.5 

(body) 

MR.  4.0 
(haad) 

12.82 

4.10 

2«l 

21.0 

15.0 

(body) 

MR.  10. 0 
(haad) 

15.12 

d.W 

19Q 

15.5 

18.0 

(body) 

MR.  11.0 
(haad) 

17.70 

4.» 

•  This  is  calculated  from  the  relatioDship  2/3  —  W. 


1.0  to  1.6  (L/Xs:!0.16  to  0.25).  Peak  absorption  for 
k\\  L  and  H  \  \  L  orientatioas  occurs  at  successively 
higher  frequencies,  with  kL  for  these  configurations  on  the 
order  of  L/2h  where  2h  is  the  length  of  the  equivalent 
prolate  spheroidal  body*  along  the  minor  axes  (see  Table 
I  for  numbers  illustrative  of  90-4(X)-g  rats).  These  values 
of  kL  for  resonance  in  the  different  orientations  are  ob¬ 
tained  approximately  by  requiring  that  for  a  sphere  [4], 
[5]  the  shortest  circumference  for  the  lossy  creeping  waves 
launched  at  the  center  of  the  shadow  plane  be  0.5  X. 

:i)  In  the  frequency  region  immediately  below  reson¬ 
ance,  absorption  diminishes  rapidly  with  frequency.  An 

”  dependence  similar  to  that  for  spheres  is  observed. 

4)  There  is  an  excellent  correlation  of  the  results  ob¬ 
tained  from  the  parallel-plate  waveguide  with  the  free- 
•spacc  irradiation  results  Ql  j,  [2]].  The  only  exception  is 
that  maximum  absorption  in  the  free-space  case  occurs  for 
body  size  of  about  0.4  X,  i.e.,  at  approximately  twice  the 
frequency  observed  in  the  parallel-plate  waveguide  situa¬ 
tion.  The  reason  for  this  may  be  that  while  a  body  isolated 
in  free  space  is  required  to  be  approximately  0.4  X  for 
resonance,  the  same  condition  is  met  in  the  presence  of  a 
ground  plane  (such  as  in  the  parallel-plate  waveguide) 
with  only  a  0.2-X  body  and  its  image  in  the  ground  plane 
acting  together. 

CORRKLATION  WITH  LETHALITY  RESULTS 

A  study  [6]  has  just  been  completed  at  Walter  Reed 
Army  Institute  of  Research,  Washington,  I).  C.,  on  the 
lethality  of  100-12.5-  and  380-420-g  rats  and  25-3.5-g  mice 
when  expiKsed  to  microwave  radiation  at  710,  9H5,  1700, 
24.50,  and  3000  MHz.  The  choice  of  frequencies  was 
dictated  by  the  anechoic  chamber  facilities  available  at 
Walter  Reed  Laboratories.  The  time  to  convulsion  was 

'  The  dimension  26  for  sn  equivalent  prolate  spheroid  is  calcu¬ 
lated  from  the  expression  2/3  wb>Lp  »  W,  where  W  is  the  weight 
of  the  body  and  p  is  the  averam  mass  density  of  nearly  1  g/cm*. 
The  nominal  circumference  2ir6  so  calculated  has  been  found  to 
correspond  well  to  the  average  circumference  actually  measured  for 
rats  of  weights  varying  from  90  to  400  g  (see  Table  1). 


measured  for  E\\L  and  H  \  \L  configurations.  At  each 
frequency,  eleven  animals  were  studied  for  each  of  the  two 
orientations  picked  to  correspond  to  the  maximum  and 
minimum  absorbing  configurations  of  Fig.  1.  The  lethality 
data  correlate  well  with  our  results  [1],  [2]  obtained  from 
the  modeling  experiments.  In  particular,  the  following 
points  of  agreement  are  noteworthy. 

1)  The  time  to  convulsion  in  the  ^  II  ^  orientation  is 
alw’ays  shorter  than  that  for  the  H\\L  orientation.* 

2)  On  the  low-frequency  side  of  resonance,  the  time  to 
convulsion  rises  very  rapidly,  corresponding  to  a  fast 
diminishing  absorption  in  this  region  (see  Fig.  1). 

3)  The  minimum  time  to  convulsion  for  ^7  ||  L  orienta¬ 
tion  occurs  for  kL  4.4.  From  Table  I,  the  nominal  L/2h 
for  rats  of  different  weights  is  on  the  order  of  4.4.  The 
lethality  results  are  consequently  in  agreement  with  the 
parallel-plate  waveguide  results  for  resonant  absorption 
at  kL  L/2b. 

4)  For  the  E  \\  L  configuration,  minimum  time  to  con¬ 
vulsion  is  obtained  at  frequencies  such  that  the  length 
from  the  animal’s  eyes  to  the  base  of  the  tail  is  approxi¬ 
mately  0.4  X.  This  is  in  agreement  with  the  free-space 
experiments  described  in  the  following  (.see  Fig.  2)  and 
also  in  p  J 

SCALING 

To  evaluate  the  validity  of  scaling,  free-space  heating 

oeriments  were  repeated*  at  710  MHz  for  saline-filled 
prolate  spheroidal  bodies  (L/2h  =  6  corresponding  to 
humans)  of  major  length  varying  from  7.6  to  25.4  cm. 
The  temperature  increase  for  a  5-min  exposure  to  100- 
mW/cm*  radiation  is  plotted  in  Fig.  2.  From  the  amount  of 
energy  absorbed,  the  relative  absorption  coefficient  [73  S 

*  The  only  exception  to  this  observation  is  40O.|g  rats  at  2450  and 
3(X)0  MHz.  For  these  animals,  such  frequencies^  are  4-5  times 
higher  than  the  resonant  frequency  of  the  E  jj  L  orientation.  A 
considerably  reduced  total  power  deposition  under  these  circum¬ 
stances  may  well  be  the  reason  for  this  anomaly. 

'Experiments  were  previously  done  at  1700  MHz  with  smaller 
sise  prolate  spheroidal  bodies  with  L  of  4.8,  6,  7.2,  9.6,  and  14.4  cm 
(see  (11). 
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Fig.  2.  TemiierAture  rise  in  s&line-filled  prolate  spheroids  with  an  aspect  ratio  I//26  of  6.  Frequency,  710  MHz;  ir> 
radiation  time,  S  min;  incident  power  density,  100  mW/cm‘;  polarization,  £  ||  L. 
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Fig.  3.  Relative  absorption  coefficient  5  for  different  size  prolate  spheroids  {L/2b  »  6  in  each  case).  Polarisation; 


is  calculated*  and  is  plotted  in  Fig.  3.  Maximum  absorp¬ 
tion  is  observed,  as  before  [1  ],  [2J,  for  bodies  of  major 
length  on  the  order  of  0.38-0.48  X,  and  whole-body  deposi¬ 
tion  of  about  3.5  times  as  much  as  that  called  for  by  the 
shadow  area  is  measured.  This  verifies  experimentally  the 
concept  of  scaling.  Hodies  scaled  down  in  all  dimensions 

*  The  relative  absorption  coefficient  S  is  defined  as  the  absorotion 
cross  section  ditrided  ny  the  physical-shadow  cross  section.  Power 
absorbed  in  milliwatts  =  (4180  X  temperature  increase  of  the 
body  X  weight  in  namsl/iiradiation  time  in  seconds.  This  divided 
by  the  incident  field  intensity  in  milliwatts/square  centimeter  gives 
the  absorption  cross  section. 


result  in  peak  absorption  frequencies  scaled  up  by  the  same 
factor  and  vice  versa. 

CORE  TEMPERATURE  OF  RATS  EXPOSED 
IN  DIFFERENT  ORIENTATIONS 

In  order  to  dramatize  the  polarization-dependent 
variability  of  the  EM  hazard,  250-g*  Wistar  rats  were 
exposed  to  l(K)-mW/cm*  fields.  The  rectal  temperature  of 
the  animal  was  measured  under  irradiation,  using  the 

'  This  particular  weight  was  selected  to  observe  the  L  0.4-X 
condition  for  mammum  absorption  in  the  free  space. 


_  .93  _ 


1 


_ L- 
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Tlae  in  minutes  after  the  onset  of  radiation 


(0 


Fig.  4.  Core  temperature  of  250-g  rate  exposed  to  710-MHz  radiation,  (a)  Orientation  £  1 1  L;  incident  power  density, 
100  mW/cm’.  (b)  Orientation  AT  ||  A;  incident  power  density,  100  mW/cm’.  (c)  Orientation  H  ||  incident  power 
density,  100  mW/cm’.  (d)  Orientation  E  ||  incident  power  density,  50  mW/cm*. 


relatively  nonperturbing  liquid-crystal  fiberoptic  tem¬ 
perature  probe  [8].  The  digital  voltmeter  reading  as  a 
function  of  time  was  recorded  and  is  given  in  Fig.  4(a), 
(b),  and  (c)  for  the  three  orientations.  The  core  tempera¬ 
ture  read  off  the  calibration  chart  of  the  probe  is  marked 
alongside  the  curves  in  Fig.  4.  For  an  incident  power 
density  of  .'iO  mW/cm*  4(d)],  the  power  deposition 
in  the  £  ||  /./  orientation  was  still  much  higher  than  that 
for  the  other  two  configurations  at  100  mW /cm*.  The  in¬ 
crease  in  animal  temperature  as  a  function  of  time  for  the 
two  power  densities  in  the  E\\L  orientation  is  plotted  in 
Fig.  .5.  F  or  this  configuration,  the  relative  absorption 


coefficients  S  of  3.8  and  4.26  are  calculated  from  the 
initial  slopes  of  the  temperature-rise  curves  for  100-  and 
.“iO-mW/cm*  incident  field  intensities,  respectively.  This 
may  be  compared  to  S  on  the  order  of  3.5  for  resonant-size 
saline-filled  prolate  spheroids  (F'ig.  3). 

DISTRIBUTION  OF  POWER  DEPOSITION— 
“HOT  SPOTS” 

To  determine  the  distribution  of  power  deposition,  an 
18.4-cm-tall  figurine  of  the  female  body  has  been  used  as  a 
mold  to  form  a  cavity  in  styrofoam.  The  size  of  the  figurine 
has  been  selected  such  as  to  allow  resonance  (Fig.  3)  at 
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Irradlarion  tt«c  in  alnutea 

FiR.  5.  Core  temperature  increase  of  25(Vr  Wiatar  rats  exposed  to 
710-MHz  free-space  radiation  (E  ||  1.  orientation). 


FiR.  6.  PhotORraph  of  the  parallel-plate  radiation  chamber. 


710  .MHz  for  frcc-spacc  radiation  (/y/X~0.44)  and  al  a 
froriuoncy  about  half  a.s  much  in  the  parallcl-platc  radia¬ 
tion  chamber.  In  ord('r  to  .simulate  ground  effects  on  lOM 
power  ab.sorption,  a  larger  parallel-plate  radiation  chamber 
(i'ig.  (>)  has  been  de.signed  and  fabricated.  For  the 
chamber,  the  central  working  area”  con.sists  of  a  copper 
plate  of  fid.o-cm  width  separate<l  by  a  2r).4-cm  clearance 
from  the  ground  [)late  of  1 17-rm  width.  The  overall  trans- 
mi.s,sion  length  of  l!).S  cm  is  occupied  by  two  symmetrical 
tapt'red  sections  of  (il-cm  axial  length,  the  central  working 
area  of  fil-cm  length,  and  the  7..'j-cm  end  .sections  con¬ 
nected  to  I J( S/U  input  and  output  coaxial  connectors. 

I'Vom  FM  field  theory,  a  body  reduced  by  a  factor  of  0  in 
all  dimensions  may  be  used  to  obtain  RF  abs<)rption 


*  The  uniformity  of  fields  in  the  c.entral  working  area  was  estab- 
lished^by  heating  a  Ht-cm-long  .saline-fillerl  prolate  spheroid  in  the 
E  II  //  orientation.  For  the  spheroid  placer!  at  the  center  of  the 
platen,  a  temperature  increase  of  1 .2.3°C  was  measiirnl  for  an  inci¬ 
dent  field  intensity  of  49.10  mW/cm*  at  .'100  MHz.  By  comparison, 
a  temiierahire  increase  of  I.I.^’C  at  tS.SS  mW/cm’  was  observed 
for  the  same  spheroid  placed  10.2  cm  off  the  center  line. 


characb-ristics  of  thi‘  full-.scale  body,  provide.d  an  ir¬ 
radiation  fre(|uency  scaled  up  by  a  factor  of  /3  is  used. 
However,  it  is  neee.s.sary  that  the  comidex  p<‘rtiu'tivity 
(fi  —  jajuitn)  represented  by  the  .sealed-down  model 
correspond  to  the  value  at  the  lower  frerpieney  charac¬ 
teristic  of  the  wholi’  body.  With  these  .scaling  |)recautioiis, 
the  distribution  of  [xrwer  deposition  in  the  rerluced  sc.iie 
model  is  identical  (though  tin'  magnitudes  are  higher  by  a 
factor  of  0)  to  that  of  the  full-scale  body.  Since  the  pro¬ 
jected  fre(]uency  of  maximum  ab.sorption  for  adult 
humans  is  ..n  the  order  of  ().'>-7.a  .MHz  in  free  space  and 
about  half  :i.s  much  in  the  presence  of  ground  effer-ts.  the 
biological-phantom  material  that  is  used  should  have  the 
dielectric  properties  of  humans  in  that  fre()uency  range. 
A  composition  [tl]  of  the  material  with  d.2f)  percimt  Nadi, 
.S.74  percent  Super.st\iff  (obtained  from  Wliamo  Manu¬ 
facturing  dompany,  San  flabriel,  California),  and  S7.0 
percetit  water  has  a  measiirc'd  dielectric  constant  of  60  and 
a  conductivity  of  4.2!)  mho/m.  For  u.se  at  a  frequency  of 
IkOO  MHz,  this  corre.sponds  toe,  =  66  —  j22-'>.('>.  Ry  com¬ 
parison,  for  mu.scle,  skin,  and  tissues  with  high  water 
content,  t,  =  KK  —  ./2.50  at  .50  MHz  [4].  Becau.se  of  the 
relative'  c.losene'ss  of  the  two  dielectric  constants,  the 
mixture  of  this  compo,sition  was  used  to  fill  the  figurine¬ 
shaped  cavity  for  experinK'nts  in  the  parallel-plate  radia¬ 
tion  chamber.’  Temperature  under  irradiation  w.as  re¬ 
corded  using  th(>  liquid-crystal  tenrperaturc  probe  [s,^. 
From  the  increase  in  temperature,  the  ab.sorbed  powe'r 
de'n.sity  in  milliwatts/gram  was  calculated  from  thi' 
rel.-ition.ship  [7]  (4180  X  temperafaire  increase) /irradia¬ 
tion  t  ime  in  seconds.  The  absorption  coefficient  a,”  defined 
as  milliwatts/gram  of  absorbed  power/milliwatts /stjuart' 
centimeter  of  incident  field  intensity  for  various  parts  of 
the  body,  has  been  calculated  and  is  given  in  Table  II.  For 
E  II  fj  (irientation,  the  strongest  intensity  of  power  deposi¬ 
tion  is  observed  in  the  neck  area  of  the  body.  For  the 
den.sity  of  power  absorption  this  is  followed  by  the  shins, 
the  thighs,  the  chest,  the  eyes,  and  the  pudendal  region,  in 
that  order.  In  Fig.  7  the  thermographic  record  of  the 
temperature  before  and  after  I  min  of  free-space  irradia¬ 
tion  with  a  field  intensity  of  100  mW/cm’  at  710  MHz  is 
shown.  1’he  pattern  of  power  deposition  bears  a  remarkable 
resemblance  to  the  distribution  obtained  in  the  parallel- 
plate  radiation  chamber. 

In  order  to  obtain  a  quantitative  information  on  the 
power  depiKsition  under  free-space  irradiation  conditions, 
a  1 2.1 -cm  resonant-size  (L/\  ~0.4)  figurine  was  used  at 
9.S.5  MHz.’  A  field  intensity  of  100  mW/cm’  was  used  to 

’  For  free-space  irradiation  experimente,  the  composition  of  the 
material  used  is  9.2,5  percent  NaCI,  8.2.5  percent  Superatuff,  and 
82. .5  percent  water.  This  mixture,  having  a  measurer!  dielectric, 
constAnt  of  66  and  a  conductivity  of  10.09  mho/m,  corresponds  to 
a  c  'mplex  relative  permittivity  r,  of  66  -  f2.'>.')  when  used  at  710 
MHz. 

'  The  parameters  8  and  a  are  relaterl  in  that  8  =  (o  X  weight  of 
the  Ixxly) /(shadow  area  of  the  body).  For  animals  and  ot.her  Ircxlies 
of  irregular  shape,  .S  is  a  more  difficult  [larameter  to  calculate  be¬ 
cause  of  the  difficulty  in  evaluating  the  .shadow  area  of  these  bodies. 

’  K.  I,.  Hunt  of  the  Dmartment  of  Microwave  Research,  Walter 
Reed  Army  Institute  of  Research,  Washington,  D.  C.,  participaterl 
in  these  experiments. 
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TABLE  II 

Measured  Distribution  OF  Power  Deposition  in  an  18.4-cm-Tai.l 
Biological-Phantom  Human  Figurine  (Parallei,-Plate 
Waveguide  Radiation  Chamber) 


Frequency 

1 

1  Pudendal 

HHz 

l./» 

Neck 

Shin 

ly.. 

[  Region 

£||  L  orientation 

230 

0.141 

2.893 

240 

0.142 

3.381 

0.788 

0.135 

220 

0.16b 

1.927 

300 

0.184 

5.097 

2.236 

0.840 

0.442 

320 

0.146 

3.162 

130 

0.213 

6.261 

1.86} 

JflO 

0.2H 

4.  168 

1.148 

14S 

0.242 

6 .'  i  3  J 

lij}  I.  orientation 

193 

0.242 

0.094 

fill  1.  orientation 

100 

0. 184 

0.1 10 

Note:  Symbols  with  arrows  above  are  represented  by  boldface 
italic  in  text. 


(b) 

Fig.  7.  The  therino|(raphic  recording  of  the  central  section  of  an 
I8.4-em  biological-phantom  fifpirine  before  and  after  free-spaee 
plane-wave  irradiation  for  1  min  at  100  mW/cm’  at  710  MHr, 
(L/X  ~  0.44).  (a)  Before  exposure  to  radiation,  Jt.5°C  baseline; 
3°C  full  scale,  (b)  After  irradiation  in  the  E  ||  L  configuration; 
waves  incident  on  the  front  aide  of  the  body,  O’C  baseline;  10°C 
full  scale. 

spiHHl  Up  the  heating  process  and  thus  reduce  the  heat 
diiTusion  to  ciRiler  areas  of  the  body.  Temperature  under 
irradiation  was  recorded  using  the  liipiid-crystal  ternjK’ra- 
ture  probe.  Due  to  the  extremely  large  energy  deposition 


rate  for  the  neck  region,  expo.sure  times  on  the  order  of  !)  s 
were  used  when  measuring  temperature  inerea.ses  for  this 
part  of  the  body.  Irradiation  times  of  lD-20  s  were  used 
for  measurements  in  other  [larts  of  the  body.  The  measured 
rates  of  tem|)erature  ri.se  are  used  to  calculate  the  a’s  for 
various  parts  of  the  body  in  a  I.?.")-!!!  man  and  these  are 
given  in  Table  III.  The  patterns  of  pow’er  deposition  ob¬ 
tained  from  both  the  parallel-plate  waveguide  radiation 
chamber  and  free-spaee  irradiation  are  very  similar.  For 
the  E  II  L  resonance  condition,  the  very  maximum  power 
depo-sition  is  observed  for  the  neck  region  of  the  body. 
Under  both  free-spaee  and  parallel-jilate  waveguide 
irradiation  conditions,  an  a  parameter  of  aiiproximately 
20  times  the  whole-body  average  is  ob.served  for  the  neck 
region  of  the  body. 

A  point  of  concern  is  the  much  lower  values  observed  in 
the  [larallel-plate  radiation  chamber  as  compared  to  the 
values  meiLsured  under  free-spaee  irradiation.  'This  may  be 
seen  by  comparing  the  values  in  Tabli'  II  to  those  of 
column  of  Table  HI.  The  fact  that  a  larger  scaling  factor 
/S  is  used  for  the  12.1-em  figurine  than  for  the  liS.4-cm 
figurine  should  account  only  for  a  factor  of  1..T2  in  going 
from  Table  II  to  Table  HI.  Somewhat  ri'duced  value's  in 
Table  1 1  will  also  be  expected  because  of  lower  field  in¬ 
tensities  and  the  resultant  heat  diffusion,  particularly 
from  hot  spots.  In  order  to  create  a  more  satisfactory 
representation  of  the  real-life  ground  effects,  a  monopole- 
above-ground  radiation  has  been  designed  and  is  eurrentl.v 
under  construction  at  our  laboratory.  The  new  radiation 
chamber  ip  sketched  in  Fig.  8.  I'he  radiator  consists  of  a 
45°  corner  reflector  in  conjunction  with  a  quarter-wave 
monopole  above  ground.  Thi'  gain  of  such  an  antenna  has 
been  calculated  [10]  and  is  tabulated  in  Table  I\\  For  the 
antenna-to-corner  spacing  of  0.65  X,  the  calculated  radia¬ 
tion  pattern  of  the  anti'nna  system  is  plotted  in  I'ig.  9.  The 
field  intensity  at  a  far-field  distance  of  6  ft  for  8.50  W  of 
radiated  power  may  be  as  high  as  60  mW/em’*.  Such  field 
infensities  are  adeipiate  for  the  iilanned  experiments.  The 
radiation  chamber  is  being  (^instructed  with  the  provision 
for  alterations  for  use  at  different  freiiuencies  in  the  800- 
1000-MHz  band.  Such  a  radiation  chamlx'r  is  not  easily 
amenable  to  wide-band  whole-bodv  absorption  nn'asure- 
ments  such  as  have  been  possible  in  the  parallel-plate 
transmissumdine'Hsing'a  netANork  analyzer.  The  monopole- 
above-ground  radiation  chamber  will,  however,  offer  a 
real-life  simulation  of  ground  efTeets  on  IlM  power  absorp¬ 
tion  in  man  and  animals. 

A  significant  result  of  the  ('xperiments  with  biological- 
phantom  figurines  is  that  a  high  power  depo.sition  in  the 
neck  region  is  observed.  Power  absorption  coeffieient  a  = 
6.7<)  for  the  neek  region  for  E  [j  /.,  orientation  is  over  17 
times  larger  than  the  average  value  for  saline-filled  pro¬ 
late  spheroids  of  the  same  overall  dimen.sion  (see  Table  V 
for  whole  body  or  average  a  for  a  19-em  saline-filled  pro¬ 
late  spheroid  at  different  frixiuencii's) .  For  fc  || //  and 
H  II  L  orieiifafions,  the  rate  of  heating  is  minimal  (Table 
II)  and  this,  for  the  neek  region,  is  a  factor  of  40  to  50 
times  smalh'r  than  the  rale  for  the  E  ||  L  configuration, 
loir  thi'se  orientations  the  results  once  again  are  in  agrei'- 
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TABLE  III 

Rates  of  Energy  Deposition  in  Various  Parts  of  the  Body 
Measured  Data  and  Projected  Rates  of  Energy  Deposition  in  a  1.75-m  Man 


of 

cbo  lody 

985  HBa  Raaonant  Slaa  (12.1  ca:  L  o.4  *) 
Blologleal-Phaotoa  Plgurln* 

1.75 

Kan.  67.9  KHt 

Rata  of  Taaparatura  , 
lacraaaa  at  100  ■W/oa^ 

(U/ksB)/(«W/ca^) 

Rata  of  Energy  Dapoaltlon 

For  Incident 
Flelda  of 

10  aU/ca^ 
W/kga 

For  Incident 
Field#  of 

50  aU/ca^ 
V/kga 

■ 

Uppor  pare 

128.3 

89.38 

6.16 

61.62 

306. 10 

LoMOT  part* 

«S.2 

45.42 

3.13 

11.31 

156.57 

Shlaa 

9.5 

6.62 

0.46 

4.56 

22.80 

Thlgha 

10. a 

7.52 

0.52 

5.19 

25.95 

Cbaat 

8.64 

6.02 

0.41 

4.15 

20.75 

Pudaadal 

lagloB 

3.2 

2.23 

0.15 

1.54 

7.7 

*yaa 

0.9 

0.63 

0.043 

0.43 

2.15 

Az«* 

? 

7 

7 

7 

1 

Note:  E  II  L  resonance  condition  (L/X  =  0.4). 

•  The  rate  of  temperature  rise  for  the  lower  part  of  the  neck  is  considerably  lower  because  of  its  vicinity  close  to 
the  torso  which  provides  a  large  heat  drain. 


-10  "jB 


Fig.  9.  Calculated  radiation  pattern  at  the  ground  plane  of  a  mono¬ 
pole  in  a  45°  comer  reflector. 


TABLE  IV 

Calculated  Field  Intensities  Due  to  a  Quarter-Wave 
Monopole  in  a  45°  Corner  Reflector 


Dlatanca  to  the 
Corner  of  the 
Reflector 

Driving  Power 
lapedance 

Calculated  r.aln 
of  the  Antenna 

Field  Inten#1rv 
at  R  •  6  feel 

a 

I 

Ohaa 

aU/ca^ 

0.55 

16  4  J56.5 

75 

62.4 

0.6 

23.5  *  ]64.6 

79.3 

66.0 

0.65 

3?  ♦  365.0 

82.7 

68.9 

0.7 

47.5  4  j6l.4 

73.1 

60.9 

0.75 

59.5  4  J59.7 

70.87 

59.0 

Note:  Transmitter  Power  »  350  W. 
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Absorption  Parameter  a  or  a  1!)-cm  Saline-Killed  Prolate 
Spheroid  as  a  Fl  nction  or  Kreoi  ency 


Frequency 

HKt 

Nuab«r  of 

LA 

{•v.r.gO 

260 

1 

0. 132 

.2943 

260 

3 

0.163 

.  3097 

IQQ 

\ 

0.V9Q 

.2989 

}20 

2 

0.203 

.3214 

360 

3 

0.229 

.4089 

380 

‘ 

0.241 

.4908 

393 

‘ 

0.231 

.  3?86 

(b) 


Fi(S.  10.  Thermographic  reconiinK  of  the  central  section  of  an  18,4- 
cm  biological-phantom  figurine;  free-.space  irradiation  of  710 
MHz  (temperature  of  the  body  before  irradiation  ~~7  i°C).  (a) 
After  irrailiation  in  the  k  \\  1^  configuration  for  10  min  at  100 
mW/cm’;  baseline  H^C;  'd'C  full  scale,  (b)  After  irradiation  in  the 
H  II  A  configuration  for  10  min  at  100  mW/cm’;  baseline  8°C; 
3°C  full  scale. 

nicnt  with  those  olitaiiied  from  free-space  irardiation.  For 
up  to  10  rnin  of  expo.sure  to  710-.MHz  free-siiace  (ields  of 
intensity  100  mW/ein*.  negliKible  modifications  of  the 
temperature  pattern  were  observed,  and  this  is  illu.strated 
in  FiR.  10. 

The  thermographic  record  of  the  tem|)erature  liefore 
and  after  1  min  of  free-space  irradiation  {£  ||  />  orienta¬ 
tion)  at  24.j0  MHz  is  shown  in  Fig.  11.  .Ma.\iiTHim  heating 
at  ‘JdoO  .MHz  (T/X  =  l..">),  a()))roximately  il.7.')  times  the 
resonant  freiiuem^v,  is  still  ob.served  in  the  neck,  tlu'  anus, 
and  the  legs,  ,\s  ex|)ected,  the  temperature  rise  is  con¬ 
siderably  lower  than  that  for  710-MIIz  radiation.  ,\l.so.  (he 
torso  is  heati'd  hardly  at  all,  which  is  a  feature  distinct  from 
the  power  distribution  at  the  re.sonant  fre(|Uency  of  710 
M  Hz  f  Fig.  7 ) ,  The  ((uantitative  values  of  a  mea.sured  with 


(b) 


Fig.  11.  The  thermographic  recordings  of  the  central  section  of  an 
I8.4-cm  biological-phantom  figurine  before  and  after  free-space 
plane-wave  irradiation  for  ]  min  at  100  mV\'/rm*  at  2450  XHlz 
(£/X  =  1.50).  (a)  Before  exposure  to  radiation, ^8. 0°C  baseline, 
1°C  full  scale,  (b)  After  irradiation  in  the  E  ||  L  configuration; 
waves  incident  on  the  front  side  of  the  body;  lO'C  black  level; 
1°C  full  .scale. 

th<‘  li(|uid-crystal  temperature  jirobe  for  a  .somewhat 
srntiller  .size  figurine  are  given  in  Tabit'  VI.  The  hottest 
parts  of  the  body  are  tht'  iit'ck,  the  legs,  and  the  arms. 

The  large  power  deposition  observed  in  (he  neck  region 
of  the  figurines  has  prom|tted  us  to  look  for  similar  efTect.s 
in  animals.  Initial  meastin'meiits''’  conducted  with  the 
liquid-crystiil  probe  implanted  .subcutancoti.sly  in  the  neck 
region  of  anesthetized  HSS-g  long  Kvans  rats  have  shown 
a  =  0,47!)  at  dtiO  MHz  as  against  a  of  O.dS  for  the  deep 
rectal  region.  For  dead  animals,  on  the  other  hand,  the 
corresponding  numbers  are  a  =  O.So.u  for  the  neck  region 
and  a  =  0.231  for  the  core  of  the  animal.  Further  experi¬ 
ments  are  currently  in  progress  to  ascertain  the  nature  of 
neck  heating  and  its  biological  implications.  Since  reports 
in  the  literaturi'  suggest  that  continued  low-level  exposure 
to  microwave  irradiation  in  both  humans  [ll]  and  ani¬ 
mals  [123  produces  fatigue  and  hyiioactivity,  the  po.ssi- 
bility  exists  that  this  may  be  linked  to  the  heating  in  (he 
neck  region.  The  heating  in  the  neck  area  may  produce 
changes  in  circulating  blood  hormones  such  as  thyroxine. 
In  future  experiments,  blood  thyroxiiu'  levels  will  be 
examined  by  serum  thyroxine  and  protein-bound  iodine- 
level  assay  techniques, 

SUMM.ARY 

The  Rl'"  power  deposition  is  found  to  vary  significantly 
with  orientation  and  with  frequency.  The  strongest  ab¬ 
sorption  is  found  for  waves  polarized  along  the  long  di- 

“  These  ex|icrjmeiil.s  arc  being  performer!  by  J.  D’Andrca  of  our 
lalmratory. 
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.  TAHI.E  VI 

£  !!  /„  Okf-Resonance  Condition  (/^/X  =  ().9sr)) 
Measi'red  Data  and  Proje«ted  Kates  of  Enkrgv  Deposition  in  a  1.7')-.m 


2<.^0  HMt 

17.1  c-n  Btol--^t>di-l*hdnl->ci  Flxui 


h: 


Region  r>t  Rat*-  of  T«-Bper^turt- 


'  (W/kgm)/laWfca  ^ 


1 


J.^0 
».  70 

I 


!  Thighs 
•  Chest 


0.42^ 

O.i 


O.Of- 

0.02 

O.Ul 


:n.  iJ.m  ^  t'.f  lii.'idei.: 


d  h-.  ,  n..'. 

.'1  i  11.(1', 


O.dfl 

o.:i 


Man 


iiicnsion  of  bodies  at  frciiucncics  such  tliaf  thi'  major 
length  is  approximately  0.4  times  tlu'  wavi’length  of 
radiation  for  bodies  in  free  spare.  I’eak  absorption,  in  the 
pre.senee  of  ground  I'fTeets,  is  observed  at  fre((ueneies  about 
one-half  as  much  as  bodies  isolated  in  free  siiace.  A  .sig¬ 
nificant  result  of  the  experiments  with  biological-phantom 
figurines  is  that  a  high  rate  of  energy  depo.sition  is  ob¬ 
served  for  the  neck  region.  Power  ab.sorption  coefficient  a 
for  the  neck  region  for  E\  \  L  orientation  is  about  20  times 
larger  than  the  whole-body  average.  This  may  be  bio¬ 
logically  significant  and  is  being  investigated  further. 
Initial  experiments  with  anesthetized  and  dead  rats  have 
indeed  demonstrated  a  higher  power  deposition  in  the  neck 
region  than  in  the  core  of  the  animals,  f'or  figurin(‘s  for 
kljL  and  ffHE  orientations,  the  rate  of  heating  is 
minimal,  and  for  the  neck  region  is  a  factor  of  40  to  50 
times  smaller  than  the  rate  for  the  E  l!  L  resonance  con¬ 
dition.  The  frecpiencies  of  maximum  ab.sor|)tion  for  adult 
humans  is  projected  to  be  on  the  order  of  0.5-75  .MHz  for 
free-space  conditions  and  about  half  as  much  in  the 
presence  of  ground  effects. 
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Numerical  Calculation  of  Electromagnetic 
Energy  Deposition  for  a  Realistic 
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Abstract — Numerical  calculadons  of  ahsorited  eneri^y  deposition  have 
heen  made  for  a  Mock  model  of  man  that  is  defined  with  careful  attention 
given  to  the  biometric  and  anatomical  features  of  a  human  beini;.  (  aku 
lated  post-resonant  absorption  and  distribution  of  energy  depirsition 
through  the  body  have  better  agreement  with  esperimental  results  than 
prcvloSB  cakulatioas  made  using  less  reallstk  models. 
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I.  ISIROIM  I  IION 

1  HAVI’  used  numerical  methods  with  a  realistic 
m.'del  to  calculate  the  deposition  of  electromag¬ 
netic  energy  in  man.  A  block  model  of  man  was  chosen  to 
allow  niaximum  freedom  in  defining  bolh  the  shape  and 
content  of  the  model.  We  have  used  the  mvtment-method 
solution  of  the  electric-fielii  integral  equation  with  a  pulse 
function  basis  and  della  functions  for  lesltng  11),  |2).  1  he 
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Fig,  1.  Improsid  mudel  uf  niiin 


folkiwing  sectuins  of  this  paper  tonlain  hrief  i.li.-'.i. riplii 'iis 
of  the  improsed  model  and  luimerisai  ^leIhl'd^.  ,i^  well  .i' 
a  comparison  of  the  calculations  ssuh  olhei  niimencal  and 
experimental  results. 

II  f)l  s(  Kit’iios  111  nil  Miilii  I 

We  have  given  careful  attention  to  hinmeiiic  and 
anatomical  diagrams  in  ilefmmg  the  shape  and  innient  of 
a  model  having  improveil  realism  X  total  nl  INI)  cubical 
cells  of  varuHis  si/es  was  used  to  obtain  a  be-.!  fil  uf  the 
contour  on  diagrams  uf  the  5()th  percentile  'laudaril  man 
(3|.  Si/es  and  placement  uf  the  celK  aie  shuui,  m  I  le.  1. 
Height  of  the  mudel  is  I  ts 

Partitioning  of  bune.  fal.  skin.  musi.le.  lung  tissue,  air. 
heart,  brain,  kidnev.  liver  .invl  spleen  ihiuiieh  the  cells 
was  done  with  the  .lul  ut  an.iluiuKal  crus-,  sevliuns  14). 
Whole-bodv  volume  fractiuiis  uf  e.ii  h  iivsue  iv[se  are 
agreement  with  pubhsheil  values  I  he  vutunie-vveighted 
complex  permittivitv  of  each  veil  w.is  fuuml  living,  proper¬ 
ties  reported  in  the  literature  |('|  |.'|, 

We  have  forced  ;>  plane  of  svmmelrv  between  the  left 
and  right  bode  halves  in  selecting  the  si/e.  plaicment.  and 
Volume-weighted  complex  permittivities  of  the  velN  I  he 
symmetrv  has  been  utih/ed  so  that  the  vi/e  of  the  m.iirix 
used  in  the  numeric.il  volutions  is  deteimiiud  bv  on  n-iK 
rather  than  the  total  of  lf<ll  cells  Ifie  svmmeliv  m.iv  onlv 


MVS 

be  ulih/ed  fur  incident  licUb  in  which  buih  /.  .old  A  are 
con'amed  in  the  pi  me  of  scmmclr'.. 


III.  Soilllosol  nil  \1  M  Kl\  1,1.11  .v  1  11  ).v 

Since  I'lie  plane  of  svmmetrv  is  used  and  there  are  three 
unknown  field  cumpuiieiits  m  each  cell,  the  matrix  is  270 
bv  2'’<f  complex,  which  currcspunds  to  a  system  of  s40 
simuilaneous  equations  in  540  real  unknovsiis.  Soiuliun  of 
the  matrix  equation  is  complicated  by  the  matrix  being 
both  full  and  asvinmctnc  I  tih/ing  svmmctiv.  or  cells 
iiaving  different  si/es  ui  |iei  nuitiv  iiies.  would  e.ich  be 
suflicient  ic' dcstii's  svmmetrv  of  ihc  m.tlrix 

We  have  used  the  successive  uv  errel.ix.ilion  iStiRl 
method  j4|  fui  uer.itive  solution  of  the  iiiainx  equ.itii'ri  It 
is  km'wn  that  a  ni.itiix  having  ,it  least  one  complex 
eigenv.ihie.  as  m  the  present  cise.  will  luit  allow  conver¬ 
gence  wilti  SOR  unless  all  cieeiivalues  h.ive  le.o  p.iits  wnh 
magnitudes  less  than  .iuiiv  (  oiiv ei ceiic e  ot  S()R  h.is  the 
further  restriction  tii.P  tlie  leiaxalaui  p.M.imete:  must  he 
less  th.ii!  .1  critical  value  th.ii  depends  up.m  the  m.itrix  |U| 
We  have  found  that  ^onvereence  woh  S(1R  requires  a 
rciaxatau!  factor  less  ih.ui  umiv  .it  lo  MM/.  ,ind  the 
critical  value  decre.ises  witli  uk  re.ismg  (reipiencv  until  it 
reaches  zero  near  on  \lfi/  (  onveiceiice  with  SOR  is  not 
possible  at  frequcncpc's  evceeding  '*o  \1)i,  wilii  ihe  present 
model  of  man.  sci  we  n.oe  been  ‘  'ixid  g.  use  Oaiiss 
Jordan  and  other  ni'niteoilive  mclhi'd'  j  loj  l:i  the  present 
problem  when  SOR  mav  be  used  in  I'iace  ol  nomteralive 
methods,  the  computer  time  o  tvpicdlv  reduced  bv  about 
a  facte'r  of  four. 

We  li.ive  oidercd  numbering  p'‘  .elN  m  the  model  in 
such  .1  w.iv  that  the  diffeietice  helween  the  numbers  for 
tw.>  ceils  is  .Ml  ii'.c  rcMsiiie  tunciioii  oi  the  di'l.ince  between 
th.e  cells.  We  h.ive  foimed  llie  m.tlnx  using  .m  arras  of 
It  ,»  1  blocks  for  each  cell  p.iir  which  is  msidc  out  friuii  tlie 
structure  used  previcuislv  1 1  1 1.  Ihe  two  ch.mges  pist  de¬ 
scribed  have  c.iused  the  m.ilnx  ti'  luive  a  r.ipid  decre.tse  in 
m.ignitucle  of  the  elements  wnh  increasing  disiancc  from 
the  di.ieonal.  I  i>r  all  liccpiencies  used  so  far  with  the 
realistic  model  I'f  m.ui  (lb  laK)  MH/i.  the  iinlv  matrix 
elements  h.ivmg  apprec  i.ible  m.igmiude  .ire  cv'nt.nned  in  a 
band  .iboul  Ihe  diagiui.d,  which  .nnounls  lo  .ipproxi- 
nialelv  lb  percent  id  the  eiuiic  m.ilnx,  Ihe  ipproxim.ition 
c>f  forced  b.mdmg  m.o  he  .ippiied  l.i  such  ,i  matrix  to 
allow  either  a  large  reduc  tion  in  cosi  of  compul.itu'ns  i>r  a 
significant  iikicmsc  in  number  of  cells  |I2],  We  have  not 
vet  implemented  forced  b.indmg,  but  we  esiim.ile  that  the 
method  will  allow  enough  cells  so  ih.it  reii.ible  cnlcula- 
lions  max  be  i"  ide  foi  ,i  rc.ihsiic  model  of  man  liv 
frequencies  as  high  .o  24‘'()  Mil/ 

l\  (  1  isv  I  1(1 ,1  M  I  (  I ivsini  R  V novs 

Nurieric.il  solutums  using  ,i  pulse  function  basis  can 
oi.lv  be  cxicl  it  there  is  uo  v.iri.ition  ol  /  within  e.ich  cell 
ol  the  niocle.  S.'liilions  h.ouie  l.irge  diflerences  in  calcu- 
l.iled  cneiL'v  depoMiion  t.u  .uli.icenl  celK  must  h.ive  signi- 
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fisanl  variation  ot  E  vviihin  the  cclh,  which  suKireslc  lack 
''f  ci'iiv creence.  I’leviouslv  reported  solutions  for  block 
models  of  man  have  a  ratui  of  2dh  ;  1  for  energs  deposi¬ 
tion  ;n  one  p.ni  of  adiacent  cells  at  10  MH/  |  I  j.  I  he 
arrangenient  and  different  sues  of  cells  in  the  present 
moUel  cause  the  maximal  ratio  of  energv  deposition  for  a 
pair  of  adjacent  cells  to  be  8  :  1  at  the  same  frequency. 

If  d  IS  the  side  of  a  cubical  scii  and  k  represents  the 
magnitude  of  the  complex  propagation  vector  within  the 
cell,  then  substantial  sanation  of  the  electric  field  must 
occur  witnin  the  cell  if  A  A  ■  \  h  [13)  I’htis  the  pulse 
function  .ipptoximation  mas  mu  he  justified  at  freqiien 
cies  above  2(kl  \1H/  with  the  new  model,  suggesting  an 
increasing  eiior  .it  higher  frequencies.  Volunie-weighling 
of  the  complex  jierniil'is  its  w  ithin  each  cell  mas  also  he 
unjustified  at  frequencies  above  21X)  MHz'  though  the 
lesults  app  ai  reasonable  to  at  least  500  Mfl/. 

.A  numerical  solution  usinu  a  pulse  funetiim  basis  re¬ 
sults  in  a  suieie  value  representing  i  withiii  each  cell.  It  is 
possible  to  use  the  /,  values  to  calculate  1,  2u/f  £'*  for 
each  cell,  and  to  use  a  voitime  average  to  estimate  the 
speciiic  absorption  r.ile  (.h.AR)  l.arge  numbers  <>f  cells 
must  be  used  in  order  to  find  aeeur<ile  values  of  S.AK  b> 
Mich  a  (vroeedure,  lor  example,  "ur  calculations  of  the 
SAR  ot  ,1  12-eni  nuisde  cube  at  I  Mil/  show  an  ernu  of 
percept  wiih  k  cells.  2<i  percent  with  27  cells,  and  20 
percent  w  ith  k4  cells.  !  he  delta  functions  used  for  testing 
e’lfoiee  the  miegral  equation  ,il  the  center  of  each  cell  so 
tfi.it  the  calcui.ited  /.  >  a  hits  are  most  representative  of  the 
ceil  centers  Inspection  of  the  solutions  suggests  that  the 

c.i:  value-  ot  /  hzive  .ippremahlv  less  erri'r  than  occurs 
IP  ihe  .S  \  K. 

l;  lliere  is  mucti  variation  ol  /  witnin  a  scattering  body, 
tlici!  even  if  vve  had  exact  v.ihies  of  /.  .it  the  cell  centcis, 
.ipl'reei.ihie  errot  wiuild  be  expected  in  the  ealciilated 
S.AR  Ihe  5.  sU|sposeil!'.  exact,  local  values  of  energv 
deposition  copstnute  sample-  of  an  olherwi-e  unknown 
population  so  that  if  large  eell-Io-eell  variation  is  present, 
-lati-l.'cal  p;  > 'cedures  mav  he  u-ed  to  c.ilc  uhile  eonlidenee 
limits  for  values  of  wtn'le-hodx  or  part-hvids  S.AR  (14| 

I  or  man  at  kn  Mil/  with  /,  /.  and  A  front-to-baek.  the 
XO-percent  eontidence  limits  for  whole-hods  S.AR  arc 
•  12  2  percent  id  the  mean  for  the  realistic  model  of  man. 
and  -  14.4  [lercent  of  the  mean  for  calulalions  reported 
earlier  having  somewh.il  gre.ile;  eell-Io-cell  variation  |lj. 

We  have  found  th.ii  .iccuraev  i-  improved  bs  using  a 
three-dimensional  inlerpolanl  with  the  /.  values  initially 
calculated  for  each  cell  lo  .leeounl  for  some  vif  the  sana¬ 
tion  of  /,  within  each  cell  |I51  I  nlinear  and  tnquadralic 
interpolants  h.ive  hiuh  been  used  to  estim.ite  the  vanalu>n 
>  f  /-  f'Ctweeti  ihe  ceii  center-  Ihe  mterpol.iiu  is  integralcd 
in  c.ilcui.iiine  the  S.AR,  I  or  the  12  cm  musile  cube  at  1 
Mil/  the  SAR  f.nitid  with  2’  celK  ,ind  itie  iriqii.idraiie 
tnlerpol.iril  have  eiior  viunp. liable  with  the  calcul.ition 
bascvl  on  (i4  cells  without  the  nileipolan!  Ihe  meieasc  in 
cosi  clue  lo  use  id  the  inlerpolanl  is  .ihoul  I  percent  All 


values  of  eneres  absorption  in  iho  pa(ier  have  been 
calculated  using  inlerpolanl- 
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f'lg,  2  shows  the  s.ilues  »»!  vv  holc-ht 'civ  S.AR  CLilctzlated 
for  man  in  free  sp.ice  with  1.  /  and  an  incident  inlen-il. 
of  1  inW  cm"  Points  on  l  ig  2  lepieseni  the  results  of 
numerical  si'lutiotf.  foi  ihe  realist, e  model  of  man.  bolh 
with  and  without  the  use  v'f  an  inlerpolanl.  Ihe  homoge¬ 
neous  approximation  ot  iwo-ihirds  ihe  complex  perniiiliv- 
ils  of  musc'lc  w.is  used  in  all  c.ik  illations  foi  the  figure  It 
IS  seen  in  Fig.  2  thal  the  foicliomil  correction  which  is 
made  bx  the  mtcrpiilant  increase-  .il  I'agh  ! requeue le- 
w'here  ineieasing  erreu  i-  expected. 

Values  of  whole-bodv  SAR  belovv  re-oiiaiKe  in  I  ig  2 
are  tvpiealK  wi'hiu  lb  peicent  of  values  c-iLuialed  fiU 
prcvlate-spheioidal  and  ellip-oid.ii  model-  u-ing  the  -.ime 
homogeneous  appic'Ximatu'n  of  complex  permiffiv  ;tv  ||n  ' 
Fixperunenta!  results  have  shown  ihal  the  SAR  toi  !. 
polarization  is  inveiseiv  prc'jvi'riional  to  frec|uencv  for 
frequencies  Irom  abou"  I  tvi  n  nines  tluit  of  resi'nar.ce  [  i'l. 
(18)  I'or  small  .uiiin.ils.  such  a-  liic  rat,  the  high  value  ol 
the  re-  >nanl  frequenev  allow-  a  reduced  magnitude  toi 
the  relative  pcrnultiv il\  of  li  — ue  -o  that  |sosl-iesonanl 
calculations  for  prolale-spheiancfa!  models  mav  be  m.ide 
using  the  extended  bouiid.irv  cotidilum  nielhod  l!b|.  Such 
.'esulls  shi'w  the  expeiimeiilally  obscrveci  1  f  behavior, 
but  t)ie  ealeulalums  there  cannot  be  direclls  applieil  lo 
man.  .A  1  (  dependence  o‘'  the  post-i esonant  S  AR  is  also 

expected  troin  antenna  theory.  Ihe  .inticipatcd  I  f  post- 
resonant  behavi  u  i-  evi.leni  iti  1  ig,  2  I  ather  nunieiic.ii 
soluiu'ns  suggested  a  sigmlicaiilix  f.istei  rcul  oif  between 
I  '  '  and  !'  ■’  I '  i 

Fig.  3  alscT  shows  ihe  whole-bodv  S.AR  of  man  in  fiee 
space  vvith  E  I  and  aii  ineidenl  inten-itv  of  1  mW  cm'. 
The  solid  curve  in  Fig  3  was  obiained  bv  spiining  ilie  2(' 
values  c'htained  in  numeneal  solulioiis  foi  ihe  luimoge- 
neous  realistic  nu'del  of  man.  I’oinls  on  Fig  3  represeiil 
experiinenlal  values  c'biained  for  human-shaped  figurines 
filled  with  saline  or  biological  phantom  mixtures  ( I  I’l,  1 1 S 
The  numeneal  solutions  suggesi  thal  ihe  resonant 
frequenev  is  77  Mil/  which  is  soniewhal  highei  lh;in  ihc 
value  of  t)K  7)  Ml  1/ deteriniiicd  cxpciiiiienlalls , 

An  impc'ttant  cc'iilribution  thal  can  be  expeeled  fii'iii 
numerical  solutions  for  .i  realistic  model  ot  man  is  the 
distribution  of  energy  deposition  ihiough  ihe  model.  Pro¬ 
late-spheroidal  and  ellipsoidal  nii'dels  e.iiinot  be  expected 
to  provide  an  aecuralc  description  of  the  distribulion  ol 
energv  deposition  since  the  dislribulion  is  strongls  depen¬ 
dent  upon  detiiiled  geomeirs.  When  the  homogeneous 
.tpproximation  of  complex  permittivity  is  used.  Ihe  disiii- 
bution  of  energy  dep<>siliori  in  the  realistic  model  of  man 
IS  in  gi'od  agreement  with  that  found  expennientalls  tor 
homogeneous  models  |17|.  |1X]  I.ible  I  cauit.nns  numeri- 
eal  and  experiment, il  values  for  the  distnbulum  of  energv 
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exajk'eraied  neck  con'-lPclc  m  Neck  lieatm^  fo'  die  do’K  ■-  ai'oi-it  ler' 
'imes  [hdi  c'i’  ihe  expected  f-T  mai; 
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Fig.  4  L.ocal  SAR  values  (exalts  kx;  per  mLL'  cm'i  ■  !CK’  for  htniioge- 
ncous  mode!  of  man  with  verlic  d  pola  ization  at  80  MU/, 

deposition  for  man  near  resonance  in  free  space.  The 
deposition  of  energy  in  the  regions  of  Ihe  e\e.  heart,  thigh, 
and  calf  has  been  severely  underestimated  in  previous 
calculations  made  with  less  realistic  models  ( 1 1,  ( I  1 1. 

When  the  inhomogeneous  complex  pennitlivitics  .ire 
used  with  the  model,  a  change  of  less  th.in  ?  percent 
typieallv  occurs  in  the  whole  hiuiy  .ixer.ige  S.AR,  but  the 
distribution  of  energv  depostitioii  is  ch.inged,  f  igs.  4  and 
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fig  *■  I'Va:  S.-VR  xalucs  xc  j  x-'  ■”,'xk  .m-i  ■  ioi  'nn-nm,, 

geru'i'us  nK'dei  '.'l  man  aith  xcrtix.ii  ■'  .I'nmi:.'!;  ,ii  8"  Mfi/ 


r  lilasirate  ’he  ci.stdhuiion  of  ihs  tped  eiiergx  in  in.in  .it 
yo  \1H/  in  free  -p.ice  [<■:  itic  iioir,. igc'ieoii-,  and  inhoino- 
geneous  'ni'Jeis,  re'.pe-.ti\e:\  One  liilferer.ce  n  ',11.1’  tile 
inhomogeneous  model  has  reduced  .liisx  .-I'tuni  in  regions 
with  High  bone  cx'ntenl.  l  arlier  workers  have  not  repoiied 
results  ivT  inhoiiiogencMUs  models  ii.i'.ing  a  celi-bx-ced 
representation  of  the  dielectric  properties 

1  L  I  )Nl  I  l  SI'  IVs 

Nuniericai  calcuLnix'ns  of  absxirbpj  energx  depk's'tiv'ii 
have  been  made  for  .1  model  of  man  dial  h.is  been  made 
with  careful  altentum  to  both  biometric  and  anati.'mical 
diagrams.  ,\  udai  of  I  SI)  cells  was  used  in  the  model, 
which  IS  a  relativelv  large  number  f'ir  numerical  solutions 
but  facilitates  accurate  mi'dehng. 

Calculated  post-resonant  absorption  and  distributu'n  I'f 
energy  deposition  through  the  body  have  been  found  to 
have  better  agreement  with  experiiuental  results  than  pre- 
vuius  calculations  made  using  less  lealistie  models. 
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l  ine  struelure  in  the  whole-body  resonanc  curve  lor  radio  rreguency  energy  deposilion  in  man  can  be 
aiiribuied  to  pan-body  resonances.  As  lor  head  resonance,  which  occurs  near  7.'l)  MH/  in  man.  the  ah 
sorplise  cross  section  is  nearly  three  times  the  physical  cross  section  ol'  the  head.  The  arm  has  a  prominent 
resonance  at  150  MH/.  Numerical  solutions,  antenna  theory,  and  experimental  results  on  animals  have 
shown  that  whole-body  energy  deposition  may  be  increased  by  50  percent  or  more  because  of  multiple 
bodies  that  are  strategically  located  in  the  field,  tmpirical  equations  forSARs  are  also  presented  along 
with  test  data  for  several  species  of  laboratory  animals.  Barbiturate  anesthesia  is  sufficiently  disruptive  of 
thermoregulation  that  STs  of  colonic  temperature  yield  energy  dose  values  in  several  mammals  that  com¬ 
pare  quite  favorably  with  those  based  on  whole-body  calorimetry. 


I.  INTRODUCTION 

Some  of  the  experimental  results  on  absorption  of  RF 
energy  that  have  been  obtained  from  animals  and  from 
models  of  man  have  recently  been  evaluated  in  the  light  of 
numerical  calculations  on  an  improved  model  of  man 
[Hagmann  el  al.,  19771,  Two  new  effects  that  were 
predicted  by  numerical  solutions  and  were  confirmed 
experimentally  are  part-body  resonances  and  multibody  ef¬ 
fects.  Both  are  described  in  this  paper  for  the  first  time. 
New  empirical  equations  for  the  SAR  for  EIIL  orientation 
are  also  given  along  with  test  data  obtained  from  several 
vertebrate  species. 

:.  PART  BODY  Rf  SONANtT  S 

Figures  1  and  2  show  pan-  and  whole-body  SARs 
calculated  for  the  improved  model  of  man  at  various  fre¬ 
quencies  in,  respectively,  free  space  and  in  electrical  contact 
with  ground,  when  the  vector  of  the  electric  field  is  parallel 
with  the  moUers  long  axis  (Fill.)  and  the  direction  of  pro¬ 
pagation  (k)  IS  from  the  ventral  to  dorsal  aspect.  Fine  struc¬ 
ture  m  the  whole-body  SARs  is  due  to  part-body 
resonances. 

In  previously  reported  experiments  (with  figurines)  on 
whole-body  SARs  tor  man,  an  anomalous  iii'  r.-asc  in  the 
rate  ol  energy  deposition  was  observed  in  the  legion  of  470 
Mil/  lor  Kill,  and  kill,  orientations  [Gandhi  ei  at..  1977). 
I  heoretical  studies  of  a  lO-cm  isolated  spherical  model  of 
the  human  head  have  shown  that  its  first  geometrical 
resonance  occurs  near  4,50  MH/  [Jotnes  and  Spiegel.  1974). 
Our  numerical  calculations  (Figures  I  ajid  2)  for  the  im¬ 
proved  model  of  man  (Figure  .1)  with  Kill,  have  shown  that 

(  I'pyrijihi  lilt'  Aincrk'ai)  <  icophYMcal  {  'nuin 


resonance  of  the  intact  human  head  occurs  at  a  frequency 
near  350  MHz.  For  head  resonance,  numerical  solutions 
give  the  S  parameter  (ratio  of  absorptive  cross  section  to 
physical  cross  section)  of  the  intact  human  head  as  3.00  for 
man  in  free  space  and  2.66  for  man  on  a  ground  plane.  The 
magnitude  of  resonant  absorption  of  the  intact  human  head 
is  sufficient  to  cause  the  experimentally  observed  increase 
in  whole-body  SAR.  The  numerically  obtained  head- 
resonance  frequency  of  350  MHz  is  in  agreement  with  the 
experimental  data  presented  in  this  paper,  but  is  somewhat 
lower  than  the  450-  to  470-MHz  values  obtained  from 
earlier  experimental  data  and  is  also  lower  than  the  first 
geometrical  resonance  of  the  10-cm  isolated  sphere.  Further 
work  is  in  progress  to  resolve  these  differences. 

Joines  and  Spiegel  |1974)  calculate  that  the  ratio  of  the 
absorptive  cross  section  to  the  physical  cross  section  of  the 
isolated  human  head  is  approximately  1,1  at  geometrical 
resonance  near  450  MHz  and  1 .4  at  a  second  resonant  peak 
that  occurs  near  2. 1  GHz,  if  appropriate  allowance  is  made 
for  the  inhomogeneous  structure  of  a  multilayered  model. 
The  present  studies  indicate  that  the  magnitude  of  the 
geometrical  resonance  of  the  intact  human  head  is  more 
than  twice  that  of  the  isolated  human  head  for  cither  of  the 
previously  identified  rcsonanecs. 

Table  1  gives  experimental  values  of  SARs  of  figurines 
filled  with  biologically  simulating  materials  with  kill,  pro¬ 
pagation  from  head  to  toe.  For  these  experiments,  ac¬ 
curately  scaled  figurines,  25.4,  33.0,  and  40,6  cm  in  length, 
were  exposed  in  free  space  to  radiation  at  2450  MHz.  The 
whole-body  dose  was  obtained  by  a  Thermonctics  Model 
2401 -A  gradient-layer  calorimeter.  The  heads  of  the 
figurines  were  attached  to  torsos  with  a  saline-soaked  cloth 
(to  maintain  conductive  contact),  which  permitted  easy 
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FREQUENCY  CMH2) 

Fig.  1 .  Part-body  SARs  for  a  homogeneous  model  of  man  in  free  space;  EH,,  k  ventral  to  dorsal,  incident 

power  density  of  1  mW  cm 

separation  during  subsequent  measurements  of  head  SARs.  loss  by  evaporation  for  the  duration  of  its  placement  in  the 
To  determine  the  energy  dose  in  the  head  after  separation,  calorimeter. 

the  part  of  the  polyurethane  mold  that  held  the  biologically  We  believe  that  the  significance  of  the  head’s  resonance, 
simulated  head  was  wrapped  in  Saran  Wrap  to  prevent  heat  as  shown  in  Table  1,  should  be  extended  to  studies  of 
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FREQUENCY  (flHZ) 

Fig  2  Part-body  SARs  for  homogeneous  model  of  man  in  contact  with  a  ground  of  infinite  conductivity, 
Eli.  k  ventral  to  dorsal;  incident  power  density  of  I  mW  cm  '. 
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1  iji,  ?.  All  inipnncil  model  ol  man 


iniriials  F  or  animals  with  a  head  similar  in  si/c  lo  that  ot 
'he  lahoraiory  rat.  resonance  is  near  the  commercially 
'ignificani  trequency  ol  2450  MH/  The  head's  resonance 
Aas  not  obtained  in  earlier  numerical  solutions  where  the 
'Htio  ol  head  lo  whole-body  SAK  wa'  reporied  as  1.07  at 
hX)  MH/  and  1  oy  at  4(Xl  MH/  \(  hen  el  al..  1076]. 

The  Ireguency  ol  masiinal  eiieigy  deposition  in  the  arms 
s  approximately  150  MH/  when  the  arm's  leiiyrth  is  nearly 
lall  a  tree-space  wasclength.  Al  resonance,  the  arm  has  an 
.ibsorptivi  cross  section  that  is  2, .5  times  the  ptivsical  cross 
seel  ion 

'  Mt  1  iiHinry  i  tin  is 

It  'las  been  shown  that  the  energy  absorbed  m  man  it 
leson.incc  corresponds  to  the  energy  received  m  a  hall 
wave  lipnie  ol  length  apptoximalely  equal  lo  tin  Height  of 
the  nan  \(,aru/hi  el  ai.  1976,  I977|  Antenna  theory  has 
been  ised  lo  calculate  the  effective  area  (or  gain)  per  dipole 
loi  .1  pair  ol  thin  dipoles  in  hrtradsidc  (sagittal)  conligiira 
non  iK  II..  k  perpendicular  lo  the  line  joining  thi  dipole's 
.enters),  l-tiinre  4  gives  the  variation  of  effective  area  as  a 


function  of  spacing  lor  half-wave  dipoles.  Figurt  5  gives  tin 
variation  of  efteenve  area  as  a  function  of  Irequency  to- 
dipoles  of  length  l  75  meters  al  a  fixed  spacing  of  I)  65  ' 
Standard  numerical  procedures  have  been  used  to  cornputi 
the  required  values  of  mutual  impedance  of  Jipoles  {Hake' 
ei  al.,  19621 

Figure  4  stiows  i)iat  if  two  thin  hall-wa'  dipoles  an 
tangent,  each  w  ill  receive  approximately  one  hab  the  energv 
it  weiuld  receive  i|  isolated  in  free  space  I'o  separation' 
greater  than  iwe'oi  qiiee  tree-space  wavelenglli  the  amen 
has  have  little  coupling  sO  that  each  has  at  elleclive  area 
approximating  the  isolated  tree-space  value  Ai  a  spacing 
of  0.65  A.  each  dipole  will  receive  about  50  percent  more 
energy  that  11  would  it  it  vvere  isolated,  f  igiire  5  shows  that 
al  a  fixed  spacing  ol  (1.65  x.  two  dipoles  with  lengths  of  I  .'5 
meters  will  have  an  enhancement  of  abon  50  pciceni  m 
received  energy  ovei  a  fairly  wide  band  ol  frequencies 
rather  than  only  at  the  frequency  for  which  they  are  half¬ 
wave  dipoles  A  sene'  of  oscillations  and  then  a  roli-off  are 
seen  at  the  high-lrequency  end  in  Figure  ‘  Since  antenna 
theory  ha'  been  shown  to  predict  accurate  SARs  Iqi  tiumaii 
mcidels  and  lor  animal'  vvuh  and  without  nearby  reflectors 
\Gaiidhi  el  al.,  |y77|,  q  is  anticipated  that  the  variation  of 
SARs  (o!  the  Fill,  virieiilation  lor  near-resonant  larget-  as  a 
function  ol  spacing  will  be  given  by  Ihv  curve  of  F  igure  4 
The  variation  ol  SAK  wrth  irequency  lo’  biological  targets 
ihal  are  sepaTalcd  a  distance  ol  0.65  A  would  similarn  be 
given  by  the  curve  in  Figure  5. 

Numerical  solution'  have  been  obtained  lo'  iik 
multibody  effect  using  the  improved  model  ol  mar 
{Hagmann  el  ai.  I977|  I  he  solutions  have  been  obtaines 
for  the  confieiiration  ol  two  men  standing  side  by  side  h.i 
ing  Kill,  and  k  directed  venirally.  l.acF  of  symmeiiv  t'c 
ween  the  left  and  right  sides  in  each  model  has  necc’ssuaiec. 
the  use  ol  180  independeni  cells  in  llu  lorniulation  o'  h', 
matrix.  Since  there  are  three  unknown  held  componen:'  'o' 
each  cell,  the  matrix  i'  540  by  540  complex,  whish  .  o- 
responds  10  a  system  ol  1080  simuliaiieoii'  cqiiahoi'  '■ 
1080  real  unknowns.  Tola',  coniptiie-  nine  'or  each  iiiai: 
body  problem  was  19  lioms,  25  iiiiiiiiies  .viih  dcvluau.l 
usage  ol  a  PDF’-IO  digital  compute:  Weak  diag.'ii.il 
dominance,  which  is  inheien:  111  'he  iiiairix,  .■omribiiie  10 
the  conditiomiig  lo  allow  reasonable  rourid->'lt  erioi'  "iili 
such  a  large  matrix.  Nomteraiive  main'  solutions  wen  re 
quired  lor  convergence. 

The  data  points  in  l  iguie  4  and  '  lepieseru  numvnal 
solulivins  of  calculated  enliaiicenteni  11  S.-XR  lot  huiiian 
models  vlite  to  two-body  ellecis  Flu  minietic.il  solutions 
are  reasonably  ciinsisietil  win  anlcnii,-  iheorv  hu  111.11.  .it 
resonaiivt  as  is  show  n  in  Figim  4,  but  sigiiilicam  dcpaiunc 
hour  atuemia  iheorv  o  seer  al  snbiesonam  hequeiKiC' 
(I  igure  'l  All  c.ilcul.il ions  I'.ised  on  .iiileiin.l  llieoiy  have 
been  zero  order;  that  is.  the'  h.ive  no  soiiccnon  h'l  non 
sinusoidal  current  disti  ibulions  Sigiiiliv.ini  hist  .iiid 
sccond-oiiler  coriecliotis  arc  leqimed  wiien  'he  aiileim.ic 
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TABl.l:  I .  lixpcrimenial  SAR  for  kill,  propagation  f  rom  head  to  toe.  Phantom  models  of  man.  Incident 

power  density  =  10  mW  cm  ^ 
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arc  not  \ery  thin  and  the  hall  'englh  is  much  diri'erenl  from  MHz.  The  whole-body  SAR  is  45.4  percent  greater  than 

X/4  lA'/ng.  1956),  Thus,  the  i/able  lateral  extent  of  the  that  for  man  in  free  space.  In  general  the  body  halves 

model  of  man  is  thought  to  ne  responsible  for  the  dif-  closest  to  each  other  have  somewhat  less  enhancement  than 

I'erence  between  numerical  solutions  and  antenna  theory  at  the  more  distant  body  halves,  but  the  asymmetry  is  reversed 

low  frequencies.  in  the  arms. 

Some  asymmetry  is  present  in  the  distribution  of  A  greater  enhancement  in  the  SAR  by  a  factoi  as  large  as 
abiorbed-energy  density  between  the  two  body  halves  in  all  2.5  is  anticipated  for  the  inner  targets  trom  antenna  theory 

numerical  solutions  obtained  for  the  two-body  effect,  for  three  or  more  bodies,  each  of  which  is  separated  0.65 

Figure  6  illustrates  the  calculated  local  SAR  values  for  one  from  the  closest  neighbors.  Table  2  gives  the  results  ot 

of  two  models  of  man  when  the  separation  between  the  preliminary  experimental  tests  of  the  multibody  effect  with 

models  is  0.65  X  =  2.53  m  at  the  resonant  frequency  of  77  anesthetized  Long-Evans  rats.  The  configurations  used  in 
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Fig.  4.  Variation  of  gain  enhancement  as  a  function  of  spacing  of  a  Fig.  5.  Variation  of  gain  enhancement  with  frequency  for  a  pair  of 

pair  of  half-wave  dipoles  in  broadside  (sagittal)  configuration.  1.75-me»pr  dipoles  in  broadside  (sagittal)  configuration  at  a  spac- 

Daia  points  represen  calculated  enhancements  of  SARs  in  humeri  ing  of  0.65  X.  Data  points  represent  calculated  enhancements  of 
models  due  to  two-body  effect.  SARs  in  human  models  due  to  iwo-bodv  effect. 
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Fig,  6.  l.otal  SAR  values  (mW  g  '  per  mW  cm ')  x  100  for  one  of 
two  homogeneous  models  of  man  having  a  cenler-to-center  spac¬ 
ing  of  0,65  \  at  77  MH/.  Incident  power  density  of  1  mW  cm 
distance  between  closest  arms  of  the  two  models  =  0,51  X. 


Fig,  7.  Geometry  used  in  measurement  of  enhancement  of  SARs 
due  to  proximity  of  other  biological  bodies,  (a)  two-animal  ex¬ 
posure;  (b)  three-animal  exposure. 


preciably  greater  enhancement  in  SAR  than  is  obtained  in 
either  animal  of  a  two-animal  exposure.  Experiments  ai 
several  frequencies  with  varying  interanimai  spacing  need 
to  be  done  before  the  full  implications  of  multibody  effects 
on  the  SAR  are  understood. 


the  tests  are  shown  in  Figure  7.  Note  that  for  each  of  the 
animals,  an  experimental  frequency  considerably  higher 
than  its  resonant  frequency  was  used,  confirming  the  fairly 
wide-band  nature  of  the  multibody  effect.  Note  also  that 
the  central  animal  tn  a  three-animal  ensemble  has  ap- 


4.  EMPIRICAl-  EQUATIONS  FOR  \X  HOI  F-BODS 
AVERAGE  SAR  UNDER  I  REE-Sl’At  F  IRRADIATION 

The  su'  -sonant  frequency  dependence  [Gandhi  ei  a/.. 
1977)  ans  observed  1//  dependence  in  the  suprareso- 
nant  region  have  been  used  to  develop  the  empirical  equa- 


I  ABI  E  2.  Experimcntal^cnhanccmcnl  m  SAR  due  to  proximity  to  other  biological  bodiev.  Anesiheti/cd 
I  ong-F.xans  rats  with  F.lll,  at  2450  MM/.  Gentcr-to-ccntcr  interanimai  separation  =  0.65  X  (  -  7  %  cm). 
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TABLt  3.  Empirical  equations  for  SARs  in  human  models  and  in  laboratory  animals.  Free-spacc  irradia¬ 
tion.  KBI.  polarization,  resonant  frequency,/.  =  11.4  1 .  ‘  GHz, 
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lions  in  Table  3  for  ihe  whole-body  average  SAR  for 
human  models  and  for  laboratory  animals  with  E  polariza¬ 
tion.  While  the  numbers  of  biologically  simulating  models 
of  man  are  the  same  as  those  reported  earlier  [CondhietoL, 
1977),  the  coefficients  for  S.„  and  SARs  of  laboratory 
animals  are  59  percent  higher  and  are  derivable  from  anten¬ 
na  theory  for  the  capture  cross  section  of  0.38  X  dipoles  by  a 
procedure  similar  to  that  of  Gandhi  el  at.  [19761.  The  re¬ 
vised  higher  values  for  SARs  are  in  good  agreement  with 
experimentally  determined  values  (Table  4)  for  several 
animal  species.  For  the  data  in  Table  4,  the  animals  were 
anesthetized  (sodium  pentobarbital,  45  mg  kg  ')  to  prevent 
movement  and  to  maintain  thereby  a  fixed  orientation  and 
also  to  limit  pharmacologically  the  normal  thermo¬ 
regulatory  functions  of  the  animal  [Panhoff  el  at.,  1977). 
The  experimental  SAR  was  calculated  from  the  increase  in 
colonic  temperature  (cloacal  temperature  in  the  case  of 
lizards  and  birds)  after  three  minutes  of  free-space  irradia¬ 
tion  at  100  mW  cm For  the  500-g  rat  and  the  120. 3-g 
dove,  freshly  euthanized  animals  were  used  and  the 
calorimeter  was  employed  to  determine  the  whole-body 
dose;  the  procedure  was  similar  to  that  outlined  by  Phillips 
el  al.  (1975).  The  points  to  note  in  Table  4  are;  (1)  The 
SARs  from  the  empirical  equations  for  laboratory  animals 
(for  the  supraresonant  region  shown  in  Table  3)  correlate 
well  with  experimental  measurements  on  several  species 
that  range  in  mass  from  18.8  to  2245  g,  and  in  length,  from 
8  to  44  cm.  For  these  animals,  a  whole-body  SAR  varying 
by  a  factor  of  36.1:1  is  observed  at  50  mW  cm  '  for 
2450-MFlz  irradiation.  (2)  For  ring  doves  there  is  a  signiFi- 
canl  difference  between  experimental  values  of  SARs  and 
those  calculated  using  the  empirical  equation.  The  reason 


for  this  discrepancy  is  not  clear.  Subtracting  the  mass  of  the 
plumage  from  the  whole-body  mass  of  the  birds  increases 
the  values  calculated  via  the  empirical  equation  by  a  few 
percent  but  does  not  explain  the  observed  difference  for 
this  avian.  (3)  Reliable  correlations  between  whole-body 
calorimetric  values  on  freshly  killed  animals  and  SARs 
based  on  measurements  of  colonic  temperature  of 
anesthetized  animals  confirm  the  rapid  hemodynamic 
dispersion  of  heat  in  living  animals.  Previous  measurements 
of  field  intensities  needed  for  comparable  times  to  convul¬ 
sion  [Gandhi  and  Hum,  in  preparation]  at  different  fre¬ 
quencies  has  also  demonstrated  that  whole-body  integral 
dose  and  dose  rates  are  important  parameters  in  the  study 
of  living  animals. 

S  fONfl  IJSIONS 

Numerical  solutions  based  on  an  improved  model  of  man 
have  demonstrated  that  the  deposition  of  energy  at 
supraresonant  frequencies  has  a  fine  structure  that  is 
ascribable  to  resonance  of  body  parts.  Failure  to  allow  for 
the  detailed  structure  of  the  body  of  man  through  use  of 
prolate  spheroidal,  ellipsoidal,  and  other  idealized  models 
would  cause  the  part-body  resonances  to  be  overlooked. 
The  enhancement  of  electromagnetic  energy  deposition  in 
the  head  al  head-resonant  frequencies  may  be  important  in 
future  studies  of  the  blood-brain  barrier  and  of  animal 
behavior.  Numerical  solutions,  antenna  theory,  and  ex¬ 
perimental  results  with  animals  have  shown  that  energy 
deposition  is  altered  by  multibody  effects.  A  maximum  of 
about  50  percent  enhancement  in  SAR  is  found  with  two 
bodies,  but  significantly  greater  enhancement  is  possible 
with  three  or  more  bodies,  interbody  spacing  being  critical. 
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I  All!  I  4.  Whok-  body  S \H-  -n  mW  g  '  lor  V.H.  ot 

'0\crvi’  sptv;cs  ('!  lahoja'or>  animaN  itrudiaicd  by  24'0  MU/ 
energy  ii.  :rce  ^pace.  I  reshly  -‘uthu.n/ed  aniiiiaK  and  a  incriiu)- 
nelKs  Model  2401 -A  gfadieiu-lay er  ealoriincfcr  were  U‘cd  lo 
measure  ilie  5(K)  g  ral  ai'.d  ihe  l2().3-g  do^e.  Other  measuremenl'. 
are  based  on  anesiheli/ed  animals. 
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Antenna  theory  and  the  observed  i  //  dependence  at 
supraresonant  frequencies  have  beer  u  ed  to  develop  em¬ 
pirical  equations  for  the  whole-bod,  tiAR  for  models  of 


i.ian  and  for  laboratory  animals  with  the  vector  of  the  t 
field  parallel  to  Ihe  body's  long  axis.  The  empirical  SAR 
values  arc  shown  to  be  in  good  agreement  with  the  ex¬ 
perimental  values  obtained  for  several  animal  species. 
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Numerical  calculation  of  electromagnetic  energj  deposition  in  models  of  man  vsith 

grounding  and  reflector  effects 

M.  J.  Ha^munn  am!  O.  P.  Gamthi 
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tin  se.iitering,  it  is  necessary  lo  treat  all  parts  of  ihescai- 
■  T  tie  ''ody  I'l-I  the  images  as  discrete  cells.  If  pulse  func 
"  ,  s  .  '  s  bu'is,  it  IS  nece".iry  lhai  each  s'.itnolunn 

-in;  ''  cm  tha:  ilic  electric  field  that  il  contains  iiiav 
h'-  .is-"n..,l  ,  ristant. 

mo:’-  :  , -ethod  o'liitinii  with  pulse  functions  uill 

'ehbe  ’It's'  ..i-liimn  matiix  /,  which  contains  the  ’Uikm'vrii 
r  I  stris  ti.clil  -s.thicv  i-i  s.ich  cell,  to  the  ci-himn  malii\  / 

.vhich  sonisiiits  val  I' s  of  the  incidem  electric  field  in  each 
,r  !'.  'hro.igh  ,i  ’siy.iatc  matrix  A  by  the  equalion 

-1/’.  /.  (it 

1  I  mav  he  soKed  iliiecllv.  but  a  redusOoil  ol 

malt  s  a/s  '-'ill  illovv  the  solution  lo  be  louitd  vsilh  .eretllc 

‘  :  ■  •  V  !•'  ]  I'v  ;  t,!’  \  tti-  r ..  .(H  <  1; .  mchI  I  njoii 

iXU.H  f.Mi.-:  “'t  111:  <Xi 


v^lierc  J . .  /.  arc  fu  h>  t'i  '  a..  -  ' 

ar*--  . i •  •  rc.  u  ;  . 

liKiiiix. 

(  oiiiprc"ioii  r;a\  be  k>hiaini\l  lutiini.'  m!;v  'Mc 

/’/  -araiilafKOii'  c^uaiioM-  ric.-J  i  ■  I'.  '  q:  : 

(•';  mj!  Niilui.iia  f2i  ni  il  \mi'  .  s  ■»; 

siimihanc";  N  eqaaia  . in  i!ib  .»->/  is.’  .•  .i. 

/:. 
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A 


£- 


(4) 


\^here  A  is  an  m  by  m  matrix  formally  defined  by  A'  = 
first  m  rows  of: 


A 


/. 


(5) 


In  implememing  compression,  the  computation  of 
matrix  *1  is  made  by  making  ri  modified  passes  ihn'ugli 
the  loops  that  would  normally  be  used  to  aecumiilaie  the 
scattering  matrix. 

We  have  used  compression  with  image  theory  to  treat  a 
body  with  one  plane  of  symmetry  above  :■  ground  plane,  a- 
is  shown  in  f  igure  1 .  The  field  relationships  used  for  a  bv'dy 
with  one  plane  of  symmetry  in  front  of  a  flat  lefleelor  are 
shown  in  I'igure  2.  It  is  aUe  possible  to  use  ccunpression 
with  image  theory  to  treat  a  body  in  a  C'erner  rencetor,  it 
the  corner  reflccior  has  a  dihedral  angle  equal  to  180'  .V 
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hfijs  ui! 

,  ',nu’  plan’,  o:  niiv.s-',’ 

in  tront  of  a  flat 

tcncvtor 

Kll/.  kiiT 

where  V  is  an  integer,  field  relationships  used  for  u  body 
with  one  plane  of  symmetry  in  fiom  o',  a  'X)  cornet  reflec¬ 
tor  are  shown  in  Figure  2, 


lx 


-  42  - 


J 


(.k(il  M)IN(,  AM)  kl  1  1  U  lOk  I  I  1  1  I  IS 


1  here  are  l\so  types  of  symmetry  relationships  used  for  K 
sallies  m  the  three  figures:  (1)  Tsso  points  loealed  sym- 
melrieally  in  respeet  to  a  plane  ol  symmetry  of  the  body 
that  eontains  both  k  and  F.'  will  neeessarily  have  K  values 
that  ate  mirror  images  so  that  the  two  halves  are  in- 
disiiiiguishable  under  relleetion.  (2)  Two  points  loeated 
s\  inmetneally  m  respeet  to  a  plane  of  symmetry  of  the  body 
thai  eontains  k  and  is  perpendieular  to  E  will  neeessarily 
base  I,  values  that  are  negated  mirror  images,  so  that  if  the 
values  are  exehanged,  they  will  eorrespond  to  the  solution 
lor  opposlielv  direeted  K  .  This  relationship  is  also  found 
lor  ivvo  points  loealed  symmelriealiv  in  respeel  to  a  ground 
plane  or  lo  a  flat  reHeetor  as  reguired  for  symmelrieal  mo- 
lion  ol  a  ehaige  and  its  image. 

An  aliernaiive  lo  the  use  of  image  iheory  with  eompres- 
sion  1'  10  use  speeial  (ireen's  funetions  in  plaee  of  the  free- 
spaee  Cueen's  lunelion,  Speeial  Cireen's  funetions  are 
available  that  could  be  used  to  account  for  scattering  ob- 
jeeis  loealed  near  targets  having  various  shapes  and  even 
finiie -ondueiivilies  17'<;i,  1'I71;  Buno.s,  1966|.  We  have  not 
vet  pursued  the  use  of  sueii  Green's  funetions. 

fhe  results  of  ealeulations  for  grounding  and  reflector 
effeels  are  described  in  the  next  sections. 

3.  GROUND  EFFECTS 

All  ealeulations  of  ground. ng  effects  have  assumed  that  a 
man  is  standing  on  or  above  a  perfectly  conducting  ground 
plane  that  is  infinite  in  e.xtent.  The  incident  field  is  verti¬ 
cally  polari/.ed  with  k  directed  from  the  vental  to  the  dorsal 
aspect  of  the  man.  Internal  fields  in  the  model  are  the  same 
as  those  in  one-half  of  a  double-man,  which  consists  of  the 
model  and  its  image  in  free  space.  A  180-eell  model  of  the 
50th-pcreentile  standard  man  has  been  used  for  the  calcula¬ 
tions  \Haiimann  el  a!..  1977). 

i  he  ealeululed  resonani  frequency  ol  man  standing  on  a 
ground  plane  is  47  Mil/  as  compared  lo  77  Mfl/  for  man  in 
free  space.  1  he  speeilie  absorption  rale  (SAKt  of  man  on 
lire  ground  plane  al  ihe  reduced  resonani  frequency  is  32..“' 
pereenl  greaier  ihan  Ihe  S.AR  for  man  m  free  space  al  the 
free-spaee  resonani  frequency.  At  lower  frequencies,  the 
S.AR  ol  man  standing  on  a  ground  plane  may  be  increased 
considerably  from  ihal  ol  the  same  model  in  free  space, 
fills  Increase  results  Irotn  these  freqiieticies  beitig  closer  to 
resonance  for  the  grouttded  titan,  forevamplc,  at  It)  Mfl/ 
ihe  S.AR  of  man  sianding  on  a  ground  plane  is  0.016.3  W 
kg  '  per  niW  cm  ,  which  is  about  seven  limes  that  found 
lor  Ihe  same  model  in  Iree  space. 

figures  4  and  5  illusirate.  respeclively.  the  disiribulion  ol 
absoi  bed  energy  in  man  al  10  Ml  1/  in  free  space  aiul  in  eon- 
lacl  with  the  ground  plane.  Ihe  etihaneemeni  of  local 
ctietgv  deposition  due  to  the  groundliig  effect  is  approv 
uu.iielv  ,1  factor  of  60  in  the  .ire.i  ol  Ihe  heel 

I  ignie  6  shows  ihe  pail  body  .ind  whole  bodv  SAR  lor 
Ihe  model  ol  nian  wiih  dillerenl  spaeings  from  the  ground 


fig.  4.  local  SAk  values  (\V  kg  ’  per  m\V  cm  ')  x  UK)  for 
iiihoimigeneous  model  ol  man  in  tree  space  at  11)  MH/. 


plane  al  10  MH/.  figure  6  is  eonsisleni  w  iih  ihe  experimen¬ 
tal  observation  on  models  that  a  small  separation  lhal 
breaks  electrical  eonlaet  with  the  ground  plane  is  sufficient 
to  eliminate  much  of  ihe  grounding  effecl  el  al., 

1975).  Several  ealeulalions  made  for  ihe  grounded  resonani 
frequeney  of  47  MH/  show  a  fall-off  in  magniiude  of 
grounding  effects  with  increasing  distance  from  the  ground 
plane,  which  is  similar  lo  ihe  rcsulls  al  10  MH/  (I  igiire  6). 
.A  model  composed  of  292  cells  was  made  by  dividing  each 
of  the  six  lowest  cells  in  each  leg  of  ihe  1. SO  cell  model  into 
eight  identical  cubical  cells.  Six  solutions  made  for  the 
292-eell  model  and  duplicated  with  the  ISO-cell  model  sug¬ 
gest  that  the  fall-off  of  the  grounding  effect  with  increasing 
distance  front  ihe  ground  plane  may  he  slightly  slower  lhan 
is  found  with  the  180-cell  model. 

Figure  7  shows  part-  and  whole-body  SARs  for  the  model 
of  man  standing  on  a  ground  plane  at  various  frequencies. 
The  part-body  resonances  shown  in  Figure  7  have  also  been 
observed  in  calculations  for  the  same  model  of  man  in  free 
space  [Gandhi  ei  al..  1977al.  The  enhancemei  in  absorp¬ 
tion  of  energy  due  to  the  presence  of  a  ground  plane  is 
greatest  in  the  leg  and  least  in  Ihe  head.  Figure  8  shows 
whole-body  S.ARs  for  Ihe  model  of  man  both  in  free  space 
and  when  standing  on  a  ground  plane  at  various  frequen- 
eies.  (irotindine  effects  are  most  prominent  at  low  frequen- 
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(  iii.  5,  1  oail  SAK  wiliiC'  (W  Kf  pi-t  i;)\\  ^iii  )  •  UK)  loi 
iiilioiiiuyciictni^  iiiotU'l  v'l  ni;in  ^utKiuii’  on  a  iitounj  plane  al 
III  MU/ 

L'ies.  Above  a  frequene>  of  about  2(K>  MH/,  pan-  anJ 
whole-body  SAR  values  show  little  dependence  on  the 
presence  of  a  ground  plane. 
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All  calculations  of  rcllecloi  cllecis  base  assumed  that  a 
man  is  standing  in  Iroiil  ol  a  leflectoi  that  is  perlecilv  con 
ducting  and  infinite  in  evtenl.  T  he  incidcni  field  is  seiiicalls 
polarized  with  k  directed  from  the  ventral  to  the  dorsal 
aspect  of  the  man.  One  image  is  required  for  a  flat  reflector 
and  three  are  needed  for  a  90  corner  reflector. 

The  enhancement  of  SAR  due  to  reflector  effects  is 
approximately  the  same  as  the  ratio  of  c.  fective  area  of  a 
half-wave  dipole  with  a  reflector  to  that  of  an  isolated 
dipole  in  free  space  [Gandfii  e(  a/.,  1977b].  Antenna  theory 
has  been  used  to  perform  numerical  calculations  made  for 
the  model  of  man.  Standard  procedures  have  been  used  for 
computation  of  the  required  values  of  mutual  impedance  of 
dipoles  [Baker  and  LaGrone,  1962). 

^igurcs9-12  show  the  enhancement  of  energy  absorption 
as  calculated  for  the  model  ol  man  with  reflector  cITccis 
and  the  corresponding  gain-cnhanccmcni  factors  calculated 
for  thin  dipoles  from  antenna  theory.  Values  of  enhance¬ 
ment  for  the  whole-body  asciage  SAR  have  been  used  in 
preparing  figures  9- 1 2.  but  the  curses  lor  pan-body  cncigy 
absorption  are  nearly  identical,  which  inslicales  that  the 
distribution  of  energy  ihrougtt  the  model  is  almost  iiiilepen- 
dent  of  the  presence  of  a  reflector,  fsperimcnlal  salues  m 
Figures  9-12  are  for  models  of  man  and  were  reported 
earlier  [Gandhi  ei  al.,  1977b]. 

f  igures  9  and  10  show  the  sariaiion  of  enhancement  ot 
energy  absorption  as  a  funelis'ii  ol  spacing  from  a  rcflecioi 
al  the  frcc-space  resonant  ITs'qucncy  ol  In  f  ieun  '> 

the  numerical  values  frsr  the  model  of  man  ;iie  inaike.il' 
below  the  curse  for  thin  dipoles  al  small  values  ol  sepa'.i 


DISTANCE  EFsOM  GROUND  PI  .ONE  CM. 


I  If.',  h  VV  link-  lioily  S.ARs  in  man  toi  iliflcreni  sp.n.ings  liom  the  gumiul  .it  to  Mil/.  I-  HI  .  k  scnii.il 
dorsal,  meidcni  posset  densjis  ol  I  inW  cm 
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Fig.  7.  Pan-body  SARs  for  a  homogeneous  model  of  man  on  a  ground  plane;  EIL,  k  ventral  to  dorsal, 
incident  power  density  of  I  mW  cm'*. 


tion.  Antenna  theory  requires  that  the  gain  be  small  when 
the  spacing  is  not  much  greater  than  the  thickness  of  the 
dipole,  but  such  variation  is  missed  in  the  thin  dipole 
approximation. 

Figures  1 1  and  12  show  the  variation  of  enhancement  of 
energy  absorption  as  a  function  of  frequency  for  fixed 
ratios  of  spacing  to  wavelength.  The  antenna-theory 
calculations  were  made  for  a  thin  dipole  with  a  length  of 
1.75  meters,  which  is  identical  to  the  height  of  the  model. 
Figures  1 1  and  12  show  that  the  reflector  effects  are  found 
over  a  wide  range  of  frequencies.  Values  for  the  model  of 
man  have  a  pronounced  roll-off  at  higher  frequencies, 
which  is  not  seen  in  calculations  for  the  thin  dipoles.  In 
Figures  1 1  and  12  it  can  be  seen  that  the  rate  of  roll-off  of 
the  SAR  at  high  frequencies  is  consistent  with  experimental 


;p*  -5^ 
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Fig.  8.  Whole-body  SARs  for  homogeneous  model  of  man:  EIL.  k 
ventral  to  dorsal,  incident  power  density  of  1  mW  cm  '. 


results  that  are  obtained  from  models  of  man.  Antenna 
theory  requires  that  the  gain  decrease  at  higher  frequencies 
when  the  effective  area  of  the  dipole  images  is  blocked  by 
the  physical  cross  section  of  the  dipole,  but  such  variation  is 
missed  in  the  thin-dipole  approximation. 

All  gain  calculations  have  been  based  on  zero-order 
antenna  theory;  that  is,  they  have  no  correction  for  non- 
sinusoidal  current  distributions.  SigniFicant  first-  and 
second-order  corrections  are  required  when  the  antennae 
are  not  very  thin  and  the  half-length  is  much  different  from 
X/4  [King,  1956).  The  sizable  lateral  extent  of  the  model  of 
man  is  thought  to  be  responsible  for  the  difference  between 
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Fig.  9.  Variation  of  gain  enhancement  with  spacing  for  a  half-wave 
dipole  in  front  of  a  flat  reflector.  Squares  represent  numerical 
solutions  for  enhancement  of  SARs  in  man  at  77  ,V1Hz  (L-  \  = 
0.45).  Solid  circles  represent  experimental  values  for  anatomically 
scaled  figurines  at  L/X  =  0.417. 
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Fig.  10.  Variation  of  gain  enhancement  with  spacing  for  a  half¬ 
wave  dipole  in  front  of  a  90°  corner  reflector.  Squares  represent 
numerical  solutions  for  enhancement  of  the  SAR  in  man  at  77 
MHz  (L/\  =  0.45).  Solid  circles  represent  experimental  values  for 
anatomically  scaled  Figurines  at  L/\  =  0.417. 

numerical  solutions  and  antenna  theory  at  low  frequencies 
In  Figure  12. 

The  numerical  solutions  indicate  that  for  frequencies 
near  resonance,  the  enhancement  in  SAR  due  to  a  reflector 
is  approximately  equal  to  the  enhancement  in  gain  of  a 
half-wave  dipole  with  the  same  reflector  configuration.  Ex¬ 
perimental  studies  on  models  of  man  have  also  shown  that 
the  enhancement  of  SAR  is  equal  to  the  enhancement  in 
gain  of  a  dipole,  but  the  gain  must  be  calculated  for  the 
Tinite  experimental  reflector  rather  than  for  an  infinite 
reflector  having  the  same  dihedral  angle.  For  certain  com¬ 
binations  of  length  and  width  of  the  reflecting  sheets,  a 
flnite  comer  reflector  may  exhibit  supergain  -  a  gain 
significantly  greater  than  tlut  for  infinite  reflecting  planes 
[Coitany  and  Wilson,  1958;  ScMkunoff  and  Friis,  19521. 


Fig.  1 1 .  Variation  of  gain  enhancement  with  frequency  for  a 
1.75-meter  dipole  0.75  \  from  a  flat  reflector.  Squares  represent 
numerical  solutions  for  enhancement  of  SARs  in  man.  Solid  circles 
represent  experimental  values  for  anatomically  scaled  figurines. 
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Fig.  12.  Variation  of  gain  enhancement  with  frequency  for  a 
1.75-meter  dipole  1.5  \  from  a  90°  corner  reflector.  Squares  repre¬ 
sent  numerical  solutions  for  enhancement  of  SARs  in  man.  Solid 
circles  represent  expenmental  values  for  anatomically  scaled 
figurines. 

Examples  of  the  effect  of  such  super^n  combinations  are 
seen  in  the  experimental  values  in  Figures  10-12. 

Figures  13  and  14  were  copied  from  an  experimental 
study  of  finite-size  comer  reflectors  [Coitany  and  Wilson, 
1958).  In  both  figures  the  angle  of  the  aperture  was  ad¬ 
justed  for  maximum  gain  at  a  particular  distance  of  the 
dipole  from  the  reflector.  For  reflector  lengths  greater  than 
about  0.5  X  (Figure  13),  the  aperture  angle  is  90®  and 
distance  of  the  dipole  from  the  reflector  is  1.5  X.  For  reflec¬ 
tor  lengths  greater  than  1.0  X  (Figure  14),  the  aperture  angle 
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Fig.  13.  Contours  of  constant  gain  for  a  corner-reflector  antenna 
with  the  dipole  in  third  position  for  various  reflector  sizes.  For 
lengths  greater  than  -  0.5  \.  the  aperture  angle  is  90°  and  the 
distance  of  the  dipole  from  the  reflector  is  1.5  \  [Cotranv  and 
Wilson,  1958). 
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I  u  14  (  il  n.iisl.iii:  1,11  ,1  ^oiiKT-rcl'liclor  anu-nna 

■Alt.’,  ihc  JiptiK  in  ■L’t.t'iitl  jVAitiun  S’l  '.ttiiou'  fcllc^'ioi  sizes.  Ia>i 
is’iii’liis  jiiealer  ll,ai,  1,0  A.  ihe  aperture  aiijtle  is  65  and  the 
iJisiaris'e  nl  the  dipole  Ironi  the  rellesior  is  1  2  A  (6’«r/ozii  ami 
Hilui/i.  ld5X|, 

is  fs5  ami  ilisUiiK'c  oV  llie  dipole  I  tom  the  reHceH'r  is  1 

Note  111  I  leiiie  1.’  lluti  ior  piopei  eoiiihiiialioiis  ol  Iciittlli 
stud  Aidihs  o!  Ihe  rerieeliiit.'  phiiies,  the  eiihtuieciiieni  in 
S\K  o\ei  ilie  Tree  spaee  \aliie  may  be  ^0  lo  45  pereem 
ioealei  ili.iii  llui!  U'l  inliniie  relleelois.  (  t>i  swell  swperjtain 
soiiihiiitilioiis  ilie  S,'\K  Ilia)  be  eiiluiiieed  (see  f  Igure  15)  by 
as  niueti  as  14  lo  14. dB  (;i  laeloi  ol  25  lo  2K)  rather  ilian 
12. y  dB  (19.5)  aniieipated  for  infinite  refleeiors,  or  (he  fae- 
Uii  ol  16  that  is  proje’eted  from  eonsideratiori  of  Ihe  loeal 
fields  in  the  sianding-wase  paltcrn  of  a  90  eorner  rellee- 
tor.  ,As  IS  shown  in  f  igure  Id,  (he  dimensions  2  A  x  d.5Aol 
the  relleelors  used  for  the  cxperimenlal  results  of  l  igures  10 
and  12  eorrespond  to  one  of  the  supergain  combinations. 

Supeigaiii  eonihinalions  of  length  and  widih  are  also  pre- 
seni  lor  a  65  eorner  reflector  (see  figure  14).  Antenna 
ealeulations  with  image  theory  are  restricted  to  eorner 
reflectors  having  dihedral  angles  of  '  vhere  \  is  an 

integer.  As  seen  in  I  igure  14,  gain  c.  nt  is  not 

limited  to  such  special  angles,  and  we  expect  ttia,  .  .nifieant 
enhaneements  of  the  SAK  will  oeeur  for  a  wide  variety  ol 
reflectors, 

s  (  ()N(  I  t  SIONS 

Image  theory  has  heeti  used  lo  treat  grounding  and 
reflector  effects  so  that  all  ealeulations  are  limited  lo  a 
perfectly  conducting  ground  plane,  or  to  reflectors  with  an 
infinite  extent.  F  xperimenlal  results  suggest  that  reneelors 
of  finite  si/e  mav  provide  signifieanllv  more  or  less 


enhancement  in  energy  deposition  than  that  ealeulaled  lor 
infinite  leneciiirs.  flic  efleei  of  finite  eondueiiviiy  on 
gtouitding  or  reflector  cfleets  with  models  of  man  has  not 
yet  been  evaluated  either  iheoretieally  or  experimentally, 

Ihe  decrease  in  resonant  Irequetiey  and  alteration  in 
energy  depv'siiion  due  to  grounding  may  be  explained  by 
cvinsidering  a  double-lenglh  man  in  tree  space.  If  the  model 
ol  man  had  a  cros-  section  independeiil  of  height  and  a 
complex  periiiilliv  ity  independent  of  frequency,  the 
giounded  resonant  Irequeiicy  would  be  one-half  the  Iree- 
space  lesonant  frequency,  but  a  ratio  of  11.611)  o  ohiained 
for  Ihe  model  that  has  been  used. 

Ihe  dependence  i>f  rel’leclor  el  feels  on  spacing  and  Ire- 
iiuency  are  in  agieemeni  with  the  gam  enhancemeni 
ealeulaled  for  dipoles  by  antenna  ifiei'iy  The  enhancemeni 
in  energy  depc'siiion  due  lo  lel'lecloi  ellecls  becomes  small 
lor  small  value'  ol  separation  or  lor  high  frequencies, 
which  effects  agree  '-vith  antenna  theory  itniy  if  linile  width 
ol  Ihe  dipole  is  considered. 
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ABSTRACT 

A  realistic  human  model  and  Improved  numerical  methods  have 
been  \ised  for  calculation  of  deposition  of  electromagnetic  energy. 
Unlike  earlier  solutions  both  the  average  absorption  and  the  dis¬ 
tribution  of  absorbed  energy  within  the  model  are  in  good  agree¬ 
ment  with  experimental  measurements  made  using  phantom  models.  The 
distribution  of  absorbed  energy  is  frequency-dependent  and  may  be 
explained  in  terms  of  resonance  of  the  various  body  parts.  Numer¬ 
ical  solutions  for  man  near  a  ground  plane  and  near  reflectors  are 
presented  for  the  first  time.  At  10  MHz  the  specific  absorption 
rate  (SAR)  of  man  standing  on  a  ground  plane  is  about  seven  times 
that  for  man  in  free  space.  Multibody  effects  have  been  predicted 
from  antenna  theory  and  observed  with  experiments  using  anesthetiz¬ 
ed  rats.  For  two  resonant  targets  separated  by  0.65  X  an  increased 
SAR,  170  percent  of  the  free-space  value,  has  been  observed. 
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DESCRIPTION  OF  MODEL  AND  NUMERICAL  TECHNIQUES 

Solutions  in  electromagnetics  are  facilitated  by  choosing  a  model  of 
simple  geometry.  Then  an  exact  analytical  solution  may  be  possible,  or  a 
numerical  solution  may  be  expressed  in  terms  of  easily  constructed  whole- 
domain  basis  functions.  We  have  rejected  such  methods  since  emphasis  was 
placed  on  a  realistic  model  of  man.  When  simple  geometries  are  not  possible, 
one  can  always  attempt  a  solution  by  "brute-force"  using  moment  methods  with 
a  subsectional  basis.  We  have  used  moment  method  solution  of  the  electric 
field  Integral  equation  with  a  pulse  function  basis  and  delta  functions  for 
testing^ . 

A  total  of  180  cubical  cells  of  various  sizes  were  used  to  obtain  a  best 
fit  of  the  contour  on  diagrams  of  the  50th  percentile  standard  man^.  Sizes 
and  placement  of  the  cells  is  shown  in  Fig.  1.  Anatomical  cross  sectic  3^’** 
were  used  in  partitioning  bone,  fat,  skin,  muscle,  lung  tissue,  air,  heart, 
brain,  kidney,  liver,  and  spleen  through  the  cells.  Whole-body  volume  frac¬ 
tions  of  each  tissue  type  are  in  agreement  with  published  values.  Properti.i, 
reported  in  the  literature  for  the  tissue  types®"^  have  been  used  to  calcu¬ 
late  the  volume-weighted  complex  permittivity  of  each  cell. 

Previously  reported  solutions  for  models  of  man  have  a  ratio  of  239.1 
for  energy  deposition  in  one  pair  of  adjacent  cells  at  10  MHz®.  The  arrange¬ 
ment  and  different  sizes  of  cells  in  the  present  model  cause  the  maximum  ratio 
of  energy  deposition  for  a  pair  of  adjacent  cells  to  be  8:1  at  the  same  fre¬ 
quency.  The  pulse  function  approximation  is  invalid  if  conditions  force  sig- 
nificant  variation  of  E  within  the  cells  and  a  large  difference  between  cal¬ 
culated  energy  deposition  in  adjacent  cells  suggests  the  existence  of  such 
variations. 

If  A  is  the  side  of  a  cubical  cell  and  k  represents  the  magnitude  of  the 
complex  propagation  vector  within  the  cell,  then  substantial  variation  of  the 
electric  field  must  occur  within  the  cell  if  kA  >  '‘/s  ^ ,  Thus  the  pulse  func¬ 
tion  approximation  may  not  be  justified  at  frequencies  above  200  MHz  with  the 
new  model  suggesting  an  increasing  error  at  higher  frequencies.  Volume 
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weighting  of  the  complex  permittivity  within  each  cell  may  also  be  unjusti¬ 
fied  at  frequencies  above  200  MHz  though  the  results  appear  reasonable  to  at 
least  500  MHz. 

A  numerical  solution  using  a  pulse  function  basis  results  in  a  single 
value  representing  E  within  each  cell.  It  is  possible  to  use  the  E  values  to 

•>  -t-jf 

calculate  1/2  aS  *  E  for  each  cell  and  use  a  volume  average  to  estimate  the 
specific  absorption  rate  (SAR)  .  Large  numbers  of  cells  must  be  used  in  order 
to  find  accurate  values  of  SAR  by  such  a  procedure.  For  example,  our  calcu¬ 
lations  of  the  SAR  of  a  12cm  muscle  cube  at  1  MHz  show  an  error  of  37%  with 
eight  cells,  26%  with  twenty  seven  cells,  and  20%  with  sixty  four  cells.  The 
delta  functions  used  for  testing  enforce  the  integral  equation  at  the  center 
of  each  cell  so  that  the  calculated  ?  values  are  most  representative  of  the 
cell  centers.  Inspection  of  the  solutions  suggests  that  the  local  values  of 
E  have  appreciably  less  error  than  occurs  in  the  SAR.  If  there  is  much  vari- 
ation  of  E  within  a  scattering  body  then  even  if  we  have  exact  values  of  E  at 
the  cell  centers  appreciable  error  would  be  expected  in  the  calculated  SAR. 

We  have  found  that  accuracy  is  improved  by  using  a  three-dimensional  inter- 
polant^®  with  the  E  values  initially  calculated  for  each  cell  to  account  for 
some  of  the  variation  of  E  within  each  cell.  Trilinear  and  triquadratic 
interpnlants  have  both  been  used  to  estimate  the  variation  of  E  between  the 
cell  centers.  The  interpolant  is  integrated  in  calculating  the  SAR.  For  the 
12cm  muscle  cube  at  1  MHz,  the  SAR  foxind,  using  twenty  seven  cells  with  the 
triquadratic  interpolant,  has  error  comparable  with  the  calculation  using 
sixty  four  cells  without  the  interpolant.  The  increase  in  cost  due  to  use 
of  the  interpolant  is  about  one  percent.  All  values  of  energy  absorption  in 
this  paper  have  been  calculated  using  interpolants. 


DISTRIBUTION  OF  ABSORBED  ENERGY  —  PART  BODY  RESONANCES 

The  principal  contribution  from  the  free  space  calculations  is  the  dis¬ 
tribution  in  energy  deposition.  Figure  2  shows  the  part-body  and  whole-body 
SAR  for  the  homogeneous  model  of  man  with  two-thirds  the  complex  permittivity 


of  muscle  with  E 


L  and  k  directed  front-to-back.  Calculated  values  of 


i 


FREQUENCY  (MHZ) 


Figure  2. 

Part-body  SAR  for  homogeneous  model  of  man.  Incident  Intensity  = 


1  mW.cm' 


whole-body  SAR  are  typically  within  20%  of  values  found  for  prolate-cphecoidal 
and  ellipsoidal  models*^.  Unlike  earlier  numerical  models  the  distribution 
of  absorbed  energy  within  the  model  is  in  good  agreement  with  that  found  ex¬ 
perimentally  for  homogeneous  phantom  models For  free  space  Irradiation  with 
this  polarization  near  resonance,  the  local  absorption  in  the  legs  and  neck  is 
considerably  higher  than  the  whole-body  average  while  the  torso  has  less  than 
average  absorption. 

\^hen  the  inhomogeneous  complex  permittivities  are  used  with  the  model,  a 
change  of  less  than  two  percent  occurs  typically  in  the  whole-body  SAR,  but  a 
more  significant  change  occurs  in  the  distribution  of  energy  deposition. 

Figures  3  and  4  illustrate  the  distribution  of  absorbed  energy  in  man  at  80 
MHz  in  free  space  for  the  homogeneous  and  inhomogeneous  models,  respectively. 
One  difference  is  that  the  Inhomogeneous  model  has  reduced  absorption  in 
regions  with  high  bone  content. 

In  previously  reported  experimental*^  results  on  whole  body  SAR  for  man, 
an  anomalous  increase  in  the  rate  of  energy  deposition  was  observed  in  the 

A  ^  A 

region  of  470  MHz  for  k  ||  L  and  E  j|  L  orientations.  This  has  now  been  iden¬ 
tified  as  the  first  (geometrical)  resonance  frequency  of  the  head.  Continu¬ 
ing  experiments  have  given  the  absorption  cross  section  for  the  head  region 
as  large  as  3.1  times  the  physical  cross  section,  perhaps  on  account  of  re¬ 
flections  from  the  nearby  torso.  At  the  head  resonance  frequency,  an  SAR  4.5 
times  the  average  value  for  the  rest  of  the  body  has  been  experimentally  ob¬ 
served. 

Higher  order  resonances  based  on  a  multilayer  formulation  may  exist  for 
the  head  region  at  higher  frequencies.  Initial  numerical  results  show  one 
such  resonance  at  frequencies  of  the  order  of  2000  MHz.  It  is,  however,  an¬ 
ticipated  that  the  overall  absorption  cross  section  at  these  frequencies  may 
not  be  as  large  as  that  for  the  first  resonance,  where  an  enhancement  factor 
of  3.1  has  been  observed  for  the  head  region. 

The  frequency  for  maximum  energy  deposition  in  the  arms  is  approximately 
150  MHz  (see  Fig.  2)  with  an  absorption  cross  section  that  is  2.3  times  the 
physical  cross  section. 


GROUND  EFFECTS 


Experimental  observations  of  ground  and  reflector  effects have  been 
reported  previously  but  numerical  methods  have  not  been  available  to  explain 
the  observations.  We  have  used  image  theory  to  reduce  the  problem  of  man 
above  ground  or  in  front  of  reflectors  to  a  multibody  problem  in  free  space. 
Symmetries  have  been  used  to  reduce  the  matrix  size. 

All  calculations  of  ground  effects  have  assumed  that  man  is  standing  on 
or  above  a  perfectly  conducting  ground  plane  that  is  infinite  in  extent.  The 
Incident  field  is  vertically  polarized  with  k.  directed  front-to-back  on  the 
man.  Internal  fields  in  the  model  are  the  same  as  those  in  one  half  of  a 

double-man  consisting  of  the  model  and  its  image  in  free  space. 

The  resonance  frequency  of  man  standing  on  a  ground  plane  is  one  half 
that  for  man  in  free  space.  The  SAR  of  man  on  the  ground  plane  at  the  re¬ 
duced  resonant  frequency  is  within  2  percent  of  the  SAR  for  man  in  free 

space  at  the  free-space  resonance  frequency. 

At  10  MHz  the  SAR  of  man  standing  on  a  ground  plane  is  0.0163  watts/kg 
per  mW/cm*,  which  is  about  seven  times  that  found  for  the  same  model  in  free 
space.  The  enhancement  in  SAR  due  to  the  ground  effect  is  found  since  the 
frequency  of  10  MHz  is  much  closer  to  the  grounded  znan  resonance  frequency 
(=35  MHz)  than  to  the  free-space  resonance  frequency  (=70  MHz) . 

Figures  5  and  6  illustrate  the  distribution  of  absorbed  energy  in  man  at 
10  MHz  in  free  space  and  in  contact  with  the  ground  plane,  respectively.  The 
enhancement  factor  in  local  energy  deposition  due  to  the  ground  effect  is 
about  60  in  the  heel  area. 

Figure  7  shows  the  part  body  and  whole  body  SAR  for  the  inhomogeneous 
model  of  man  with  different  spacings  from  the  ground  plane  at  10  MHz.  A  small 
small  separation  that  breaks  electrical  contact  with  the  ground  plane  is  suf¬ 
ficient  to  eliminate  much  of  the  ground  effect. 
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Figure  6.  2 

Local  SAR  values  (watts/kg  per  mW/cm  )  x  100  for  inhomogeneous  model  of 
man  standing  on  a  ground  plane  at  10  MHz. 
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REFLECTOR  EFFECTS 

All  calculations  of  reflector  effects  have  assumed  that  man  is  standing 
in  front  of  a  reflector  that  is  perfectly  conducting  and  infinite  in  extent. 
The  Incident  field  is  vertically  polarized  with  k  directed  front- to-back  on 
the  man.  One  image  is  required  for  a  flat  reflector  and  three  are  needed  for 
a  9C°  comer  reflector. 

At  65  MHz  the  computed  SAR  is  A. 87  tines  the  free-space  value  when  man 
is  0.1875  X  in  front  of  a  flat  reflector  and  16.6  times  the  free-space  value 
when  man  is  1.5  X  in  front  of  the  axis  of  a  90°  comer  reflector. 

Antenna  theory  may  be  used'^  to  calculate  the  ratio  of  effective  area 
of  a  half-wave  dipole  with  a  reflector  to  that  of  the  dipole  in  free  space. 
Such  ratios  are  within  16%  of  the  above  calculated  factors  of  enhancement  of 
SAR  due  to  reflector  effects. 

For  certain  length-to-width  ratios,  the  experimentally*^  observed  en¬ 
hancement  in  energy  deposition  is  30  to  40%  higher  than  that  anticipated 
from  antenna  theory  and  calculated  numerically  for  reflectors  of  infinite 
dimensions.  This  phenomenon  is,  once  again,  in  agreement  with  experimentally 
obtained  antenna  gains  for  finite  size*"*  comer  reflector  antennas.  Signi¬ 
ficantly  enhanced  rates  of  energy  deposition  are  projected  for  all  kinds  of 
comer  angles  (not  just  the  values  corresponding  to  180° /n,  where  n  is  an 
integer)  and  for  reflector  lengths  and  widths  that  are  no  more  than  a  frac¬ 
tion  of  a  wavelength  at  the  resonance  frequency. 

MULTIBODY  EFFECTS 

Driving  point  inpedance  values  for  a  broadside  array  of  two  half-wave 
dipoles*^  have  been  used  to  prepare  Fig.  8  which  shows  the  variation  of  ef¬ 
fective  area  per  dipole  with  spacing.  If  the  two  dipoles  are  tangent,  each 
will  receive  approximately  one  half  the  energy  it  would  receive  if  isolated 
in  free  space.  For  separations  greater  than  two  or  three  free  space  wave¬ 
lengths,  the  antennas  have  little  coupling  so  that  each  has  an  effective  area 
approximating  the  isolated  free  space  value.  At  a  spacing  of  0.65  X,  each 
dipole  will  receive  about  56%  more  energy  than  it  would  if  it  were  isolated. 
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SPACING  IN  WAVELENGTHS 


It  has  been  shown' ^  that  the  energy  absorbed  by  man  at  resonance  corresponds 
to  the  energy  received  by  a  half-wave  dipole  of  length  approximately  equal  to 
the  height  of  the  man.  Mutual  impedance,  and  hence  driving  point  impedance, 
values  are  dependent  upon  the  shape  of  the  antenna  elements,  but  we  may  ex¬ 
pect  that  the  variation  of  SAR  with  spacing  for  two  men  near  resonance  is 
similar  to  Fig.  8. 

Preliminary  experiments  using  animals  have  confirmed  the  existence  of 
multibody  effects.  We  have  used  anesthetized  adult  rats  with  E  | |  L  ori¬ 
entation  and  frontal  (broadside)  incidence.  In  tests  with  480  ±  50  gram 
rats  at  600  MHz  with  an  Incident  intensity  of  100  mW/cm^,  isolated  animals 
had  an  average  temperature  Increase  of  0.675°  C/mln.  For  two  animals  placed 
0.65  X  apart,  the  monitored  animals  had  an  average  temperature  increase  of 
1.14^  C/min.  A  70  percent  increase  in  SAR  was  caused  by  the  presence  of  the 
second  animal. 

Antenna  theory  suggests  similar  enhancements  in  SAR  for  targets  in  the 
subresonance  and  supraresonance  regions,  also  for  spacing  on  the  order  of 
0.65-0.7  X.  Furthermore,  a  greater  enhancement  in  the  SAR  by  a  factor  as 
large  as  2.5  is  anticipated  for  the  inner  targets  with  three  or  more  ele¬ 
ments. 
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Head  Resonance:  Numerical  Solutions 
and  Experimental  Results 

MARK  J  H ACIM  AN  N.  \ii  stm  k.  ii  1 1 .  OM  R  (lANl)lll.  nil  ou  .  ii  l  i ,  JOHN  A  [)'  ANDRhA.  am) 
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Abairact — We  have  used  numerical  solutioas  and  experimenis  wiih 
phantom  models  of  man,  and  experiments  with  the  Ding  Kvans  rat  to  show 
the  existence  of  head  resonance,  greatest  absorption  in  the  head  reyjon  of 
man  occurs  at  a  fre<)uency  of  about  37f  Mfiz.  Absorption  is  stronger  for 
wave  propagation  from  head  to  toe  than  it  is  when  the  electric  field  is 
parallel  to  the  long  axis.  The  highest  absorption  cniss  section  for  the 
human  head  Is  projected  to  be  approximately  .1.!i  times  Its  physical  cross 
section. 

I  iMROni  i  IION 

E:  EI.AVE  previously  reported  numerK.il  solutions 
for  the  deposition  of  eleelromagnetic  energy  m  a 
realistic  model  of  man  yvhich  shoyved  the  existence  of 
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resonances  for  btidy  parts  such  as  the  head  and  arms,  as 
well  as  for  the  whole  body  1 1  ].  I  urther  wiirk  has  been 
done  regarding  head  resonance,  since  yve  believe  the  phe¬ 
nomenon  may  be  important  in  the  study  of  behayu'ral 
effects,  bkiod  brain  barrier  permeability,  cataractogenesis. 
and  other  microyyavc  buieffects. 

Our  results  show  that  the  first  resonance  of  the  int.ici 
human  head  oecurs  at  a  frequency  of  abinit  .''“'5  Mil/  and 
has  an  .S'  parameter  tratio  of  absvirption  cross  section  to 
physical  ^r>'ss  section)  of  about  ,7 N  for  incident  plane 
yvaves  propagating  from  head  to  toe  E„irlier  calcui.itions 
for  a  spherical  model  of  the  isolated  human  head  shv'wed 
a  geometrical  resonaiwe  with  .V  =  1  1  near  4.70  MH/.  and  ,i 
second  resonance  yxith  .S' =  I  4  occurring  near  2  1  ("ill/ 
when  allowance  is  made  for  the  inhomogeneous  snucture 
by  using  a  multilayered  model  |2].  T  he  results  of  our  study 


0018-9480/ 79/0900-08(W$lX)  7?  .  1474  ICFE-: 


-  64  - 


IJ  n  1  S''*-  M  II*  '  '*•  **<  i‘  '<*  i'*''  •"  t  i  1 11  '  'H  ‘i  ‘  i  1  I  J  ••  '.  II  ‘  '  i  i’  ;  )  M  H  '  f< 


ihal  the  .ihsoriilii'ii  iii  the  he.nl  ie‘e:.'i;  !' 
hepeiiJeiit  cn  the  pie^eitee  nf  tlte  u-^t  ..f  the 

In  tile  teh.'U.ne  seal<>n^  .h  ihi'  paper  'Ai-  preseiit  ne\\ 
dai.i  ieeaikii'H:  leM.naiiee  pI'  tile  iiUai.1  iiiiaiaii  he  iri  .'!■ 
t  laieJ  aMtie  hat';  n  nr.eniai  ^aiiitiaii'  .,a>i  phani'.a.i  iii-  .1 
cK  fat  man  m  'ael!  a'  e\pei  i  men  lal  levu!!'  na  the  lane 
!  all-  iai 


II,  \  l  \ll  Kli  M  St  il  I  I  a  ea 

1  le  1  thu'lr.-.te'  a  feali'tie  nualei  .>1  m  m  a  hi.. !:  -ae  .ae-.l 
in  eai  her  eampntaiiaim  ^lla\\ln,e  1k\h!  ;ea..name  I  1 1.  I  lie 
nunie!  u  a'  ahia  nu\l  h\  hm lie  t.  'la!  a|  I  M  U  ahie a  1  i.  elN  .  .| 
eara'U'  si/e'.  te  al'i.im  .i  best  f.t  at  ihr.  I'.nt.'tit  ■  m  >1:  i- 
eranm  '.tie  .''lllh  peieenlile  ^talul,Ul.i  m.ii.  j  ']  Ali  imineti- 
eai  Nakiturli''  'ep.'ileP  m  ttU'  paper  ate  tai  a  in.'<ie!  lia'.ni'e 
ni.'fe  Jeta:'ei!  m.Hleline  >'l  the  lUak  .tin!  iie.ni  .i-  .h..\\ii  in 
F  :e'  -  a.ii'.i  '  Kept eAen l.iti. m  at  the  lesl  at  the  baiis  ia  .m 
nIiiu'.;!  t!  hie.  1,  \aie  that  tlie  imiiibe;'  .'f  -ell-.  Inm  heeii 
ilierea'ej  fraiii  12  ta  144  in  tiie  lie,a.l  'eeiaii.  ai!..l  Irani  X  ta 
.^2  in  tile  neek  le.eian  far  a  ahaie  barK  lai.il  ai  .'4(iiellA 

tA e  haNc  uvOel  a  nianient  niethail  saluiian  <‘1  the  eleeita 
field  inieierai  equ.ilian  with  a  ptilse  funetiaii  basm  and 
delta  lunetians  Im  tes'in.i;  [Ij  and  an  ii.lerpalatil  fallaw up 
far  iinpras etiieni  at  eatnereenee  |4|,  Sinee  there  are  three 
unknawn  field  er'nipanent.A  per  eell.  e\en  iisiiie  one  pi, me 
at  s\niinetr>  af  the  niadel  the  matrix  is  ■'1(1  •  5  1(»  eoniplex, 
■.ehieh  earrespands  ta  a  sssteni  af  102(1  siimilt, means  eqiia- 
tians  in  1020  real  unknawns,  Namter.itive  ni.ilnx  salii- 
tians  are  required  sinee  the  matrix  has  ,il  least  ane  earn- 
plex  eigenxrilue  havine  a  real  part  with  mae.nitude  execeel- 
ing  units  [?|,  A  talal  af  17  h  af  dedieated  usage  af  a 
PnP-|o  digital  eamputer  ssas  required  fsir  etteh  solution. 
I  artunatels.  sseak  siiagaiial  daminanee  which  is  inherent 
in  the  ni.itnx  formulation  contributes  to  the  eandiiiomng 
ta  allass  reasonable  round-off  errors  with  such  a  huge 
matrix  M/e. 

lags  4  ,i;id  .s  show  the  heati  ,md  sshale-badx  specific 
.ibsarisliaii  i  .iie  ( S  A  R )  l  alculaled  far  the  ness  model  of 
m.in  m  free  sp.ice  ssith  .m  iiKident  plane  wave  h.ixing  two 
different  palari/atians.  Note  Ih.il  he. id  lesoiiaiiee  is  more 
pranaiincesl  sshen  the  propagation  seetur  is  parallel  to  the 
hods  .ixis  (k  I.)  with  prop.igation  from  head  to  toe  than  it 
is  when  the  eleetrie  field  is  p.irallel  to  the  long  axis  ( /.  /,» 
[.xperiment.il  il.ita  far  the  l  ong  [  siiiis  r,il.  slescribed  i  i 
.Section  l\  of  this  p.iper.  ,iiso  suggest  ih.il  head  resonance 
IS  stronger  far  A  /  th.m  for  i  1.  polari/atian.  Some 
experimental  d.it.i  obt. lined  using  biologic.il  phantom 
models  I'f  man  .ire  presented  far  com|s.irison  m  the  next 
section  of  this  p.iper. 

[  igs  (i  and  7  ilhisir.ile  the  distribution  of  .ibsorbed 
energy  in  the  he.ul  region  near  he.id  reson.mce.  I’he 
presence  of  ,i  hot  spat  exists  near  the  center  of  the  head 
which  Is  .lisa  seen  in  the  distrihiitian  ,it  the  geomelrieal 
resonance  for  ,i  sphere  (2|.  I  he  distribution  shown  m  lags. 
6  and  7  is  only  approximate  since  a  homogeneous  model 
was  used  in  .ill  ealeulations. 
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Fig  2  Delailcd  nu*dcling  of  iho  nokk  tnd  head,  fruni 


I  ig  3  Ociailcvl  modeling  of  the  neck  and  head,  side  view 
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f  t’  Head  ,i:'d  'A  fi.  'ds  absDrp'u'r.  ft)r  u  .-.c  pri^pagatioi'  fr-.*m 

'liM'i  f''f  /  M’lfi'  '«•  f\n.k.  iritidcfil  «  10  mW  \fn'  0>n<'n 

'r.irkcr*'  arc  'uj-ru-fi'-a'  sttluth-ns  and  s‘'iid  marker'  arc  for  phan 
:  ■'!  :;i-  'dels  ■  d  nan 


f  ig  b  l  ocal  SAR  valuer  walls  per  kilogram  for  lOmW  'em'  imidcnt 

fields  Inner  laser  of  tells.  ^^0  Mff/.  k  heati  lo  Uk".  f  front  to  back, 
whoie-bods  average  ==  0  ft  n2  head  average -»  2  014 


III.  Kxi>l  kl\tl  M  \i  Ri  SI  I  I  ^  !  ;  »K  Pil \M  I  |\! 

I  S  (  II  \1  ,\s 

We  have  u.scu  uccuratelv  scale.!  fieurines  [fil  ■.>(  2*'.?-. 
25.4-.  .V'.O-.  anJ  4(i6-cni  leni:tii  .uii  hiclucical  phanloni 
mixture^  P).  descnhc/i  in  l  .ihle  I.  !i>  njodel  ihe  exposure 
of  man  lo  plane  sxaves  in  free  sp.ae  .Ml  sallies  of  ah- 
sorbed  dose  were  measured  iisiiiu  a  I  hermonelK  s  model 
240I-.A  graslienl  laser  caiorinu-ter  l-.‘.eh  fipiinne  had  the 
head  allaehed  |o  the  torso  b>  a  l.ixer  of  s.tline-soaked 
cit'th,  thus  mamtaininu  conduetoe  esTil-iet  but  aliosvitig 
easx  separabilux  foi  measurement  of  energs  deposition  in 
the  head,  for  determination  the  absorbed  dose  in  the 
head,  after  separation,  the  part  of  the  pols  urethane  mold 
holding  the  biological  phanloni  matcn.il  sh.tped  in  the 
fi'rm  of  the  head  was  wrapped  in  Saritii  Vs  rap  to  limn  loss 
of  water  b\  evapu'ration. 

fable  II  gives  the  experimental  sallies  of  SAR  for 
figurines  containing  biological  phantom  mixtures  and  ir- 
radiated  with  propagation  from  head  to  toe.  A  plot  s'f 
data  from  fable  11  on  lag  5  shosxs  good  agreement 
betsseen  Ihe  numerical  solutions  and  experiments  ssiih 
phantom  models  of  man.  The  values  in  Table  II  suggest 
that  for  propagation  from  head  lo  toe,  the  magmluile  iT 
absorption  m  the  head  region  near  head  resonance  is  not 
strongly  dependent  on  orientation  of  the  f,  vector 

Both  the  numerical  solutions  in  Section  II  .tml  the 
phantom  figurines  base  used  homitgeneous  models.  The 
next  section  of  this  paper  presents  experimental  data 
obtained  with  Ihe  I  ong  I  vans  rat 

IV  HXIMRIMISIM  Rl  SI  I  IS  KIR  lilt  l.OSt.  Tv  ASS 
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Table  111  gives  the  experimental  values  of  S.AR  for 
several  tests  made  with  both  freshiv  killed  and  anesthe¬ 
tized  long  Tv.ins  rats  The  number  of  measurements  [n) 
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head  res  in.mce  fur  ,i  mcijium-si/eii  i.ii  .liiiiili!  uieiii  rie.ir 
the  lest  frequenis  of  2450  Mil/. 

1  he  anesiheti/ed  nils  were  civen  a  J''-iii!: .  ki'  ilosi.'  .a 
sodium  pentohai  hitol  to  faeihtate  use  of  a  tixesl  orieiila- 
tion  .iiid  limit  therino-reizulatorv  fimitiHiis.  (  aloriniel-K 
measurements  lAere  used  to  lietermme  the  absorheil  dose 
in  ihe  freslih  killed  rats.  )  hese  nie.isurenieTiis  Asere  made 
using  a  I hcrnioneties  model  24(11 -A  Seebeck  envelope 
gradient  layer  ealorimeter,  I  he  dose  was  also  determined 
b\  measurement  of  the  rate  of  increase  in  rat  eolonie 
tempcr.ilure  anil  brain  temperature,  using  iKpiid  ervstal 
temper, ilure  optical  fiber  (I.COf  )  probes,  in  the  anesihe- 
fi/ed  rats  f  or  brain  temperature  the  I COI-  probe  was 
impl.inted  in  a  Inphine  liole  ,1  mm  posterior  to  the  Hregnia 
cranial  suture.  2  mm  l.iteral  to  the  midline  iranial  suture, 
and  6  K  mm  below  the  lop  surface  of  Ihe  cortex 

Ihe  results  in  lable  III  show  that  heaii  resonance  is 
more  pronounced  for  A  /.  than  for  /  /  orientation 

w-'hieh  was  also  noted  vsiih  Ihe  numene.il  solutions  Note 
that  the  ratio  of  hc.id-to-whole-bods  heal  eonleni  is  less 
for  the  anesthetized  r.il  ih.in  for  the  de.ul  rat.  Hlooil 
circulation  rnav  be  rerluemg  the  tel.iltse  inaemlude  of 
heating  in  the  head  region. 


\.  Si  MMAKS  \M)  I  OSl  I  I  SIOSS 

kk  e  ha\e  used  three  methods  to  study  head  resonance: 
numerical  sohuions,  experiments  with  phantom  models  of 
m.i'i.  .ind  expenmenls  with  the  l.ong  I'xans  rat.  .Ml  three 
appro.uhes  show  the  existence  of  head  resonance.  Ihe 
phenomenon  appears  to  haxe  a  greater  magnitude  fr'r 
ware  prop.ig.iluin  fixmi  head  to  toe  than  for  /■.'  /,  orienta¬ 
tion  I  or  the  former  orientation  the  head  absorption  cross 
section  as  high  as  ,b5  times  the  physical  cross  section  is 
projected 

The  numerical  solutions,  which  are  supported  by  the 
experimental  results  for  phantom  models  of  man.  suggest 
that  the  absorption  is  much  sirAUiger  than  would  be  pre¬ 
dicted  using  a  sphere  to  model  the  isolated  head.  The 
strong  dependence  on  polarization  wimld.  of  course,  also 
be  missed  using  a  spherical  imuiel. 

kk'e  believe  that  the  enhanced  absorption  in  the  head 
region  max  m.ike  head  resonance  significant  in  the  study 
of  behavioral  effects,  blood-brain  barrier  permeability, 
caiaractogenesis,  and  other  microwave  bioeffects. 

Rin  Ki  \<  ts 

(l|  (1  P  (iiindhi  and  M  J  Hafimann,  "Some  recent  results  on  deposi 

turn  ot  electromastnelic  energs  in  animals  and  in  minlels  id  man,' 
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Electromagnetic  Absorption  in  a  Multilayered 

Model  of  Man 


PETER  W.  BARBER,  membkr,  ieek,  OM  P.  GANDHI,  i  ellow.  ikee.  MARK  J.  HAGMANN, 
MEMBER,  IEEE,  ANO  INDIRA  CHATTERJEE,  STITIENT  MEMBER.  IEEE 


Abstract-\  multilayered  planar  model  is  used  to  examine  the  depen¬ 
dence  of  whole-body  power  absorption  on  the  configuration  of  surface 
layers,  e.g„  skin,  fat,  and  muscle  which  normally  occur  in  biological 
bodies.  It  is  found  that  the  layering  resonance  for  three-dimensional 
bodies  (as  opposed  to  the  geometrical  resonance)  can  be  predicted  quite 
accurately  by  a  planar  model.  Calculations  for  a  multilayered  prolate 
spheroidal  model  of  man  predict  a  whole-body  layering  resonance  at 
1.8  GHz  with  a  power  absorption  34  percent  greater  than  that  pre¬ 
dicted  by  a  homogeneous  model. 

Introduction 

ECENT  interest  in  quantifying  both  the  hazardous  and 
potentially  beneficial  effects  of  nonionizing  electromag¬ 
netic  (EM)  radiation  on  man  has  been  the  impetus  for  a  great 
deal  of  experimental  and  theroretical  research.  Of  particular 
interest  are  calculations  to  determine  the  relationship  between 
incident  power  density  and  the  resulting  absorbed  power  due 
to  whole-body  irradiation  of  man.  Recent  theoretical  methods 
that  have  been  used  include  a  perturbation  approach  ( 1 1  useful 
for  analyzing  homogeneous  ellipsoidal  models  and  an  integral 
equation  method  (2|  for  analyzing  homogeneous  prolate 
spheroidal  models.  The  absorption  characteristics  of  inhomo¬ 
geneous  irregular  shaped  models  constructed  of  cubical  sub¬ 
volumes  have  been  studied  by  volume  integral  moment  method 
techniques  [3],  [4).  In  these  whole-body  irradiation  studies, 
the  power  absorption  is  determined  as  a  function  of  the  angle 
of  incidence,  polarization,  and  frequency  of  the  incident  EM 
wave.  Usually  the  results  are  shown  as  a  plot  of  specific  ab- 
sorpuon  rate  in  watts/kilograin  versus  frequency,  where  an 
incident  power  density  of  1  mW/cm^  is  assumed.  Typically,  for 
a  given  angle  of  incidence  and  polarization,  (he  power  absorp¬ 
tion  increases  relatively  rapidly  with  frequency  to  a  resonant 
peak  and  then  slowly  decreases  to  an  asymptotic  high-frequency 
limit. 

There  is  another  interaction  which  takes  place  which  has  not 
been  considered  in  recent  whole-body  calculations.  Both 
homogeneous  and  inhomogeneous  models  have  not  accounted 
for  the  EM  interaction  due  to  the  surface  layering  of  biological 

Manuscript  received  March  13,  1978:  revised  I  ehruary  1979  This 
work  was  supported  in  part  by  the  United  States  Army  Medical  Research 
and  Development  Command  under  Contract  DAMD-I  7-74-(  -4(192 
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bodies,  eg.,  the  layers  of  skin,  fat.  and  muscle.  An  investigation 
of  these  effects  is  the  subiect  of  this  paper. 

.Absorption  Epeects  Due  to  LAVi  RiNt; 

Early  calculations  to  assess  the  biological  significance  of  EM 
radiation  used  planar  models  consisting  of  skin,  fat,  and  semi- 
infinite  muscle  layers  |5l  Later  a  two-layered  (fat  and  muscle) 
spherical  model  was  used  lb] .  These  earlier  investigations  are 
summarized  in  (7|.  Both  studies  shtmed  that  the  absorption 
characteristics  are  critically  dependent  on  the  number  and 
thickness  of  the  surface  layers.  More  recent  whole-body  calcu¬ 
lations  have  not  considered  the  effect  of  surface  layers.  One 
reason  for  this  is  due  to  the  unavailability  of  multilayered  ana¬ 
lytical  solutions  for  any  geometry  other  than  the  plane,  infinite 
circular  or  elliptical  cylinder,  and  sphere.  In  the  case  of  nu¬ 
merical  approaches,  such  as  the  volume  moment  method 
technique,  practical  limitations  on  the  number  of  cells  which 
may  be  used  prohibit  consideration  of  layering  effects. 

A  six-layered  sphere  has  recently  been  used  to  model  the  iso¬ 
lated  head  [8l-[10| .  Of  particular  interest  are  results  in  [9] 
and  (I0|  which  compare  the  frequency  dependent  absorption 
in  six-layered  spheres  with  that  m  corresponding  homogeneous 
spheres.  One  of  the  power  absorpiion  results  in  |9)  has  been 
recalculated  and  is  shown  in  Eig.  1  as  a  plot  of  absorption 
efficiency  versus  frequency.  This  and  subsequent  multilayered 
sphere  calculations  were  made  using  a  Mie  theory  computer 
program  based  on  tf  ■  mathematical  development  in  |8|  .  The 
absorption  cfficr.  s  the  total  power  absorbed  divided  by 
the  powe*-  ,id  .  i  the  geometrical  cross  section.  An 
absorplio  i.-  lency  ater  than  unity  indicates  that  the  in¬ 
fluence  ol  the  body  on  the  incident  wave  extends  beyond  the 
geometrical  boundary.  The  absorption  efficiency  for  the  homo¬ 
geneous  brain  tissue  sphere  increases  with  frequency  until  the 
maximum  is  reached  and  then  slowly  decreases.  The  amplitude 
and  frequency  of  the  resonance  are  dependent  upon  the  size, 
shape,  and  dielectric  characteristics  of  the  model.  This  max¬ 
imum  absorption  condition  can  be  called  the  t’cnmctrical 
resonance.  The  absorption  efficiency  for  the  six-layered  model 
has  an  additional  resonance  which  is  due  to  the  impedance 
matching  effects  of  the  surface  layers  and  at  this  resonance 
the  absorption  efficiency  is  almost  30  percent  greater  than  that 
at  ihe  geometrical  resonance.  This  second  maximum  absorption 
condition  can  be  called  the  layering  resonance. 
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Fig.  2.  Absorption  characteristics  of  a  six-layer  plana:  lodel.  The  five 
outside  layers  are  the  same  as  those  in  Fig.  1 .  The  innermost  layer  is 
serniinfinite  brain  material. 


FREQUENCY,  GHZ 

Fig.  1.  Absorption  characteristics  of  the  six-layered  10-cm  radius  spher¬ 
ical  head  model  of  Joines  and  Spiegel  (9|.  The  five  outside  layers 
(and  thickness)  ate  skin  (0.15  cm),  fat  (0.12  cm),  bone  (0.43  cm), 
dura  (0.1  cm),  andCSF  (0.3  cm).  The  radius  of  the  inner  brain  sphere 
is  8.9  cm. 

The  most  interesting  feature  shown  in  Fig.  I  is  that  the  pres¬ 
ence  of  layering  enhances  the  absorption  of  the  spherical  model. 
Another  feature  is  that  the  geometrical  and  layering  resonances 
appear  to  be  independent  of  one  another,  i.e.,  the  geometrical 
resonances  of  the  homogeneous  and  layered  spheres  are  almost 
identical,  while  the  layering  resonance  of  the  six-layered  sphere 
appears  merely  as  an  enhancement  of  the  absorption  over  that 
of  the  homogeneous  sphere. 

Calculations  for  a  homogeneous  and  six- layered  planar  model 
have  been  made  for  the  same  tissue  thicknesses  and  dielectric 
characteristics  as  were  used  to  obtain  the  results  in  Fig.  1 .  The 
calculations  were  made  using  conventional  planar  techniques 
[11].  The  brain  tissue  core  in  the  spherical  model  was  assumed 
to  extend  to  infinity  in  the  planar  model.  These  results  are 
shown  in  Fig.  2.  The  planar  model,  which  does  not  exhibit  a 
geometrical  resonance,  does  however  show  the  same  layering 
resonance  as  is  observed  in  Fig.  1  for  the  layered  spherical 
model.  The  location  of  the  layering  resonance  and  the  en¬ 
hancement  of  the  absorption  efficiency  due  to  the  layering  are 
almost  identical  for  the  planar  and  spherical  models.  Taking 
the  ratio  of  the  layered  to  homogeneous  absorption  efficiencies 
for  the  planar  model  in  Fig.  2  jnd  multiplying  by  the  homo¬ 
geneous  sphere  result  in  Fig.  1  gives  the  dashed  curve  in  Fig.  3, 
a  predicted  absorption  for  the  layered  sphere.  Comparing  this 
predicted  absorption  to  the  actual  absorption  shows  that  the 
layering  enhancement  of  the  planar  model  can  be  applied  to 
the  homogeneous  sphere  solution  to  predict  the  absorption  in 
the  layered  sphere  to  within  5  percent,  both  in  absorption  ef¬ 
ficiency  at  a  given  frequency  and  in  the  location  of  the  reso¬ 
nant  peaks. 
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Fig.  3.  Actual  and  planar  predicted  absorption  characteristics  of  the  six¬ 
layered  spherical  head  model  of  Fig.  1 . 


There  are  two  obvious  questions  that  neeu  to  be  answered: 
Why  are  the  geometrical  and  layering  resonances  essentially 
independent  in  Fig.  1 ,  and  why  does  the  planar  model  predict 
the  layering  resonance  of  the  head  model  with  such  accuracy? 

To  answer  the  first  question,  consider  the  geometrical  reso¬ 
nance  of  the  head  model  in  Fig.  I .  The  resonance  of  the  homo¬ 
geneous  sphere  is  dependent  orr  the  radius  and  constitutive 
makeup  of  the  sphere.  For  high-loss  dielectric  bodies  of  the 
type  considered  here,  the  peak  absorption  occurs  when  the 
ratio  of  circumference  to  free-space  wavelength  (2to/X  =  ka) 
is  about  unity.  For  the  homogeneous  sphere  in  Fig.  1,  reso¬ 
nance  occurs  when  ka  =s  0.94.  Referring  to  the  scale  drawing 
of  the  sphere  in  Fig.  1 ,  it  is  clear  that  the  surface  layers  repre- 
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sent  a  small  fraction  of  the  total  sphere  radius  (about  10  per¬ 
cent).  One  would  expect  then  that  replacing  the  outer  portion 
of  the  homogeneous  sphere  by  a  layered  segment  which  on  the 
average  has  similar  dielectric  characteristics  and  whose  thick¬ 
ness  is  much  less  than  a  wavelength  in  the  material  will  have 
little  effect  on  the  geometrical  resonant  frequency.  This  is  the 
case  as  noted  in  Fig.  I .  The  layers  have  an  independent  reso¬ 
nant  frequency,  which  in  this  case  is  approximately  four  times 
the  geometrical  resonant  frequency.  The  frequency  at  which 
the  layers  resonate  is  a  function  of  their  thickness  and  constit¬ 
utive  parameters.  Given  the  electrical  characteristics  of  the 
layers,  the  layering  resonant  frequency  is  inversely  proportional 
to  the  thickness  of  the  layers.  The  fact  that  the  layering  reso¬ 
nant  frequency  is  well  removed  from  the  geometrical  resonant 
frequency  of  the  homogeneous  sphere  accounts  for  the  inde¬ 
pendence  of  the  two  resonances.  Numerical  calculations  show 
that  when  the  layers  are  a  larger  fraction  of  the  sphere  radius, 
i.e.,  the  frequency  of  surface  layer  resonance  approaches  the 
geometric  resonance  frequency  of  the  homogeneous  sphere, 
then  the  geometrical  resonance  of  the  inhomogeneous  sphere 
is  no  longer  the  same  as  the  geometrical  resonance  of  the 
homogeneous  sphere. 

Now  consider  the  second  question  as  to  why  the  planar 
models  can  predict  both  the  resonant  frequency  and  enhance¬ 
ment  over  the  homogeneous  case  for  the  layered  sphere. 
Looking  first  at  the  resonant  frequency,  we  note  that  in  the 
planar  case  the  incident  and  transmitted  waves  are  normal  to 
the  surface.  In  the  case  of  the  sphere,  the  only  statement  that 
can  be  made  with  certainty  is  that  the  incident  wave  is  in  gen¬ 
eral  not  normal  to  the  surface.  However,  the  absorption  is  due 
to  an  internal  interaction  and  the  behavior  of  the  fields  trans¬ 
mitted  into  the  sphere  must  be  considered. 

Snell’s  law  can  be  used  to  relate  the  angle  of  the  transmitted 
wave  to  the  angle  of  the  incident  wave  at  local  regions  on  the 
surface.  Snell’s  law  for  a  wave  propagating  from  free  space 
with  permittivity  eo  into  a  region  with  permittivity  e  is  given 
by 


sin  0,  =  (co/e)'^^  sin  6, 

where  0,-  and  0,  are  the  angles  of  incidence  and  transmission, 
respectively,  both  angles  being  measured  from  the  local  normal 
to  the  surface.  The  dielectric  characteristics  of  biological  tis¬ 
sue  are  characterized  by  a  relatively  large  complex  permittivity. 
For  example,  at  the  layering  resonant  frequency  of  2.06  GHz, 
the  dielectric  constant  of  the  outer  sk'  .  layer  is  47.5-jll.4. 
Substituting  this  into  Snell’s  law  for  0,  =  90°  (grazing  inci¬ 
dence),  it  is  found  that  the  angle  of  transmission  is  only  8°  from 
the  normal.  Since  the  angle  of  transmission  will  be  even  smaller 
for  other  angles  of  incidence,  it  is  clear  that  the  transmitted 
wave  propagates  almost  normally  into  the  tissue  regardless  of 
the  angle  of  incidence,  i.e.,  the  wave  transmitted  into  the 
sphere  interacts  with  the  surface  layers  in  essentially  the  same 
manner  as  it  does  in  the  planar  model,  as  shown  in  Fig.  4. 
Therefore,  the  resonant  frequency  resulting  from  this  inter¬ 
action  should  be  the  same  whether  the  layers  reside  on  the  sur¬ 
face  of  a  sphere  or  on  the  surface  of  a  plane. 


I'i);.  4.  The  local  behavior  of  the  transmitted  wave  is  almost  identical. 
I.e.,  normal  to  the  surface,  at  (a)  planar  and  (b)  nonplanar  interfaces 
for  materials  with  a  high  dielectric  constant. 


The  fractional  enhancement  of  the  absorption  due  to  the 
layering  can  best  be  understood  by  considering  the  surface 
layers  as  a  frequency-dependent  impedance  matching  device. 
At  frequencies  well  removed  from  the  layering  resonant  fre¬ 
quency  (either  far  above  or  far  below),  the  surface  impedance 
of  the  layered  model  is  the  same  as  for  the  homogeneous 
model,  i.e.,  at  the  lower  frequencies  the  layers  are  so  thin  as 
to  have  a  negligible  effect  and  at  very  high  frequencies  the 
depth  of  penetration  is  so  low  that  the  transmitted  power  is 
all  absorbed  in  the  surface  skin  layer  which  has  electrical  char¬ 
acteristics  almost  identical  to  those  of  the  brain  material  of 
the  homogeneous  model.  In  the  region  of  the  layering  reso¬ 
nant  frequency,  an  enhanced  power  transmission  occurs,  and 
this  enhancement  is  dependent  only  on  the  surface  layering 
configuration,  and  therefore  the  fractional  enhancement  in 
power  transmission  will  be  the  same  for  both  spherical  and 
planar  models.  In  the  layered  planar  model,  the  power  trans¬ 
mitted  into  the  innermost  layer  is  completely  absorbed.  The 
fractional  power  absorption  enhancement  in  the  sphere  can 
only  be  the  same  if  all  the  transmitted  power  is  also  absorbed. 
This  will  occur  if  the  diameter  of  the  innermost  brain  material 
is  many  skin  depths  thick.  A  calculation  at  the  layering  resonant 
frequency  of  2.06  GHz  shows  that  the  diameter  of  the  inner 
brain  sphere  is  greater  than  six  skin  depths,  insuring  almost 
total  absorption. 

The  planar-prediction  procedure  which  has  been  outlined 
here  for  the  spherical  head  model  is  equally  applicable  to  arbi¬ 
trary  nonspherical  bodies,  because  the  mechanism  of  absorp¬ 
tion  enhancement  depends  only  on  the  transmitted  wave  being 
nearly  perpendicular  to  the  surface  layers,  a  phenomenon  which 
is  independent  of  the  body  shape,  Fxtensive  calculations  for 
other  spherical  bodies  and  infinite  cylinders  (H.  Massoudi,  pri¬ 
vate  communication)  have  verified  the  validity  of  the  planar- 
prediction  method.  Summarizing: 

1)  The  layering  and  geometrical  absorption  resonances  for  a 
layered  three-dimensional  object  will  be  inde|)endcnl  if  the 
layering  resonant  frequency  is  well  above  the  geometrical 
resonance,  i.e.,  the  surface  layers  responsible  for  the  layering 
resonance  are  a  small  fraction  of  the  overall  size  of  the  object, 

2)  if  the  resonances  are  independent  (a  condition  which  can 
be  tested  by  comparing  the  geometric  resonance  of  the  homo¬ 
geneous  three-dimensional  object  with  the  layering  resonance 
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of  a  semiinfinite  planar  model),  then  a  planar  model  can  be 
used  to  predict  the  absorption  layering  resonance  for  non- 
planar  geometries  provided  that  the  complex  dielectric  constant 
of  the  object  is  large  so  that  the  transmitted  fields  are  normal 
to  the  layers  and  the  skin  depth  in  the  material  is  small  enough 
that  the  power  transmitted  into  the  interior  of  the  three- 
dimensional  object  can  be  assumed  to  be  completely  absorbed. 
Note  that  I)  must  be  satisfied  before  2)  can  be  applied,  i.e., 
tlie  geometric  resonance  must  be  independent  of  the  layering 
because  the  planar  model  cannot  predict  changes  in  the  geo¬ 
metric  resonance,  only  in  the  layering  resonance.  Conditions 
1 )  and  2)  may  appear  to  be  so  restrictive  as  to  have  little  appli¬ 
cation.  however,  they  are  in  fact  satisfied  by  most  biological 
models. 


Layerinc.  Absorption  Resonance  in  Man 

It  has  been  shown  that  a  semiinfinite  planar  model  can  accu¬ 
rately  predict  the  layering  resonance  in  a  nonplanar  biological 
model.  The  interest  here  is  in  determining  the  effect  of  layer¬ 
ing  on  the  whole-body  absorption  of  EM  waves  by  man.  This 
will  be  done  by  using  a  planar  model  to  determine  the  absorp¬ 
tion  enhancement  due  to  surface  layering  and  then  applying 
the  resulting  enhancement  factor  to  whole-body  absorption 
results  previously  obtained  for  a  homogeneous  prolate  spher¬ 
oidal  model  of  man. 

In  the  numerical  calculations  which  follow,  the  dielectric 
characteristics  of  the  various  tissue  types  are  taken  from  [12] 
for  frequencies  up  to  and  including  10  GHz.  Above  10  GHz, 
the  permittivity  and  conductivity  are  based  on  the  characteris¬ 
tics  of  electrically  polarizable  molecules  [10],  modified  to 
provide  continuity  at  10  GHz  with  the  lower  frequency  dielec¬ 
tric  characteristics. 

The  surface  layers  of  man  in  general  consist  of  skin-fat- 
muscle  or  skin-fat-muscle-bone-muscle  arrangements.  The 
surface  layering  information  required  for  the  multilayered 
planar  model  was  obtained  from  published  anatomical  cross- 
sectional  data  [13] ,  [14] .  Tissue  thicknesses  were  examined  in 
79  horizontal  cross  sections  of  man.  Average  thickness  of  the 
surface  layers  of  skin,  fat,  muscle,  and  bone  were  calculated 
over  the  front  half  of  each  cross  section.  It  was  found  that  the 
surface  characteristics  could  be  represented  by  three-layer 
skin-fat-muscle  configurations  in  37  of  the  cross  sections, 
while  the  remaining  42  cross  sections  were  better  represented 
by  a  five-layer  skin-fat-muscle-bone-muscle  configuration. 
The  mean  thickness  (and  standard  deviation)  of  the  skin  and 
fat  layers  for  79  cross  sections  are  2.25  mm  (0.73mm)  and 
2.62  mm  (1.47  mm),  respectively.  Corresponding  values  for 
the  muscle  and  bone  layers  for  42  cross  sections  are  9.74  mm 
(4.32  mm)  and  9.01  mm  (4.81  mm),  respectively.  Calculations 
similar  to  those  shown  in  Fig.  2  were  made  for  each  of  the  79 
cross  sections.  The  results  were  weighted  by  the  fractional 
frontal  surface  area  which  they  represented  and  then  averaged 
to  obtain  the  layered  curve  in  Fig.  5.  Note  that  there  is  an 
absorption  enhancement  over  homogeneous  muscle  tissue  of 
34  percent  at  1.8  GHz  and  a  reduction  of  24  percent  at  5.7 
GHz 
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I  it!.  5.  Predicted  whole-body  layering  resonance  in  man  utilizinj;  a 
planar  model. 


1  ig.  6.  Whole-body  absorption  versus  frequency  for  a  model  of  man 
exposed  to  a  broadside  incident  plane  wave  polarized  parallel  to  the 
long  axis  of  the  model.  The  dashed  curve  was  calculated  using  homo¬ 
geneous  prolate  spheroidal  and  circular  cylindrical  models  (a/6  = 
6.34,  a  =  0.875  m).  The  solid  curve  is  the  planar-predicted  absorp¬ 
tion  in  a  corresponding  multilayered  model  of  man. 

Fig.  6  shows  the  whole-body  absorption  efficiency  for  both 
a  homogeneous  and  a  multilayered  model  of  man.  The  homo¬ 
geneous  result  was  obtained  using  a  combination  of  prolate 
spheroidal  and  cylindrical  models  [15].  The  specific  result 
shown  is  for  the  maximum  absorption  case  which  occurs  w'hen 
the  model  is  illuminated  broadside  by  a  wave  polarized  parallel 
to  the  long  axis  of  the  model.  The  homogeneous  results  in 
[15]  were  obtained  using  the  dielectric  characteristics  for  a 
homogenized  mixture  of  fat,  bone,  and  muscle  rather  than 
muscle  tissue  alone.  Homogeneous  calculations  for  muscle  tis¬ 
sue  alone  would  result  in  a  more  accurate  multilayered  predic¬ 
tion;  however,  these  results  are  not  presently  available.  The 
multilayered  result  was  obtained  by  taking  the  ratio  of  the 
layered  to  homogeneous  absorption  efficiencies  for  the  planar 
model  in  Fig.  5  and  multiplying  the  layering  enhancement  factor 
by  the  homogeneous  result  in  Fig.  6.  Inasmuch  as  the  multi¬ 
layered  data  used  to  generate  Fig.  5  were  obtained  over  the 
front  surface  of  the  cross  sections  of  man,  the  multilayered 
curve  in  Fig.  6  is  applicable  to  frontal  illumination  of  a  man 
model.  The  net  effect  of  the  surface  layering  is  generally  to 
increase  the  absorption  in  the  postresonance  region.  The 
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('akulatioii>  bovi'iul  10  (ill/  show  lliat  llio  two  resulls  inciik- 
when  the  deplh  ol  peiiei ration  heeomes  so  small  that  the  total 
absoiption  in  tlie  imiltilav ered  model  oeenrs  m  the  imisele-like 
skin  tissue. 

Companne  Fpe  (i  to  l  ig.  1 .  it  is  elear  th.a  the  enhaneemeni 
ot  the  absorption  due  to  layering  is  not  nearly  as  great  for  the 
whole-body  model  as  it  is  for  the  isolated  head  model.  There 
are  two  reasons  for  this,  hirst,  for  the  polarization  of  the  inci¬ 
dent  wave  shown  here,  the  long  slender  shape  of  the  prolate 
spheroidal  whole-body  model  results  in  much  greater  absorp¬ 
tion  at  the  geometrical  resonance  (absorption  efficiency  =  4.43) 
as  compared  to  the  spherical  head  model  (absorption  effi¬ 
ciency  =  1.08).  Second,  the  layering  enhancement  for  the 
whole-body  case  is  much  less  than  for  the  spherical  head 
model  because  of  the  averaging  over  many  different  layering 
combinations  in  the  whole-body  calculations  rather  than  using 
a  single  layering  configuration  as  was  done  for  the  spherical 
head  model.  If  an  averaging  scheme  were  used  in  the  head 
model  calculations,  a  reduction  in  the  layering  enchancement 
would  occur  there  also. 


('i.)N<.Lt:SIONS 

It  has  been  shown  that  a  planar  model  can  accurately  predict 
the  layering  resonance  for  a  nonplanar  geometry.  Specifically, 
results  for  multilayered  and  homogeneous  seiniinfinite  planar 
models  determine  a  layering  enhancement  factor  which  can  be 
applied  to  whole-body  absorption  results  obtained  from  non- 
layered  three-dimensional  geometries  to  predict  the  absorp¬ 
tion  characteristics  of  three-dimensional  layered  geometries. 
The  technique  has  been  used  to  predict  the  power  absorption 
characteristics  of  a  multilayered  model  of  man. 
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ABSTRACT 

The  factors  of  microwave  radiation  frequency  and  orientation  of  an  animal 
in  the  microwave  beam  are  important  parameters  in  determining  dose  rate.  Dur¬ 
ing  microwave  irradiation  these  factors  may,  on  the  one  hand,  account  for  dif¬ 
ferential  temperature  rise  in  the  rat  and,  on  the  other  hand,  produce  differen¬ 
tial  disruption  of  behavioral  performance.  Fifteen  male  Long  Evans  rats,  in 
Experiment  I,  were  exposed  to  continuous  wave  (CW)  microwave  energy  in  the  220- 
500  MHz  region  In  three  body  orientations  to  determine  the  most  absorbant  fre¬ 
quency  and  body  orientation  within  the  microwave  field  as  evidenced  by  colonic 
temperature  rise.  The  most  absorbant  orientation  was  found  with  the  animal 
length  aligned  parallel  to  the  electric  field.  In  Experiment  II,  six  male  Long 
Evans  rats  were  trained  to  lever  press  for  dry  food  pellets.  Each  animal  was 
then  randomly  exposed  while  lever  pressing  to  three  microwave  frequencies  (360, 
480,  and  500  MHz)  in  two  body  orientations  (parallel  to  the  electric  field  or 
parallel  to  the  direction  of  microwave  travel).  The  differential  absorption 
rates  due  to  the  systematic  manipulation  of  frequency  and  orientation  variables 
did  selectively  disrupt  behavioral  performance  of  the  rat.  Greatest  disruption 
of  behavioral  performance  occurred  when  exposures  were  made  with  the  animal 
aligned  parallel  to  the  electric  field  at  500  MHz.  Less  disruption  of  perfor¬ 
mance  occurred  at  360  or  480  MHz  in  this  orientation.  Little  disruption  of 
performance  was  observed  when  the  animal  was  aligned  parallel  to  the  direction 
of  the  microwave  travel  at  360,  480,  or  500  MHz. 


INTRODUCTION 

In  recent  years  the  relationship  between  electromagnetic  phenomena  and 
biological  systems  has  generated  considerable  interest.  Of  particular  interest 
have  been  the  variables  which  determine  the  dose  rate  on  the  biological  target. 
A  variety  of  investigators  [1,  2,  3,  4]  have  noted  that  the  dose  rate  for  the 
biological  organism  depends  strongly  on  such  factors  as:  (1)  average  and 
instantaneous  field  intensity,  (2)  duration  of  exposure,  (3)  size  and  shape  of 
the  organism,  (4)  dielectric  properties  of  the  organism,  (5)  wavelength  of 
radiation,  and  (6)  orientation  ot  the  organism  in  microwave  fields.  Recently 
Gandhi  [5,  6]  has  noted  a  strong  resonance  phenomena  based  on  animal  size, 
wavelength  of  radiation,  and  animal  orientation  in  microwave  fields.  Since 
these  factors  should  account  for  differential  absorption  and  temperature  rise 
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in  the  rat  during  irradiation,  these  conditions  should  also  produce  differential 
disruption  of  behavioral  performance  in  the  rat  during  irradiation.  To  investi¬ 
gate  these  factors  further,  two  experiments  were  conducted  to  determine  the 
effect  of  these  factors  on  rat  heating  and  behavioral  performance  during  irradi¬ 
ation. 


EXPERIMENT  I;  MATERIALS  AND  METHOD 
Subjects 

Fifteen  male  Long  Evans  rats  350-380  grams  in  body  weight  were  used.  All 
rats  were  deprived  of  food  to  85  percent  of  their  free-feeding  weights. 

Microwave  Irradiation  Apparatus 

The  microwave  radiation  chamber  consisted  of  two  parallel  copper  plates. 

The  lower  plate  was  2  mm  in  thickness  and  measured  120  cm  x  198.2  cm.  The 
upper  plate  was  1  mm  in  thickness  and  measured  63.5  cm  x  198.2  cm.  Both  plates 
were  enclosed  in  a  rectangular  plywood  box  215  cm  x  88  cm  x  150  cm.  The  inte¬ 
rior  of  the  box  was  lined  with  sheets  of  microwave  absorbing  material  (11.8  cm 
thick  Eccosorb).  A  ventilating  fan  was  placed  at  one  end  of  the  chamber  to 
remove  air  from  the  chamber  at  the  rate  of  310  cubic  feet  per  minute.  Addi¬ 
tional  ventilation  openings  were  placed  near  the  position  of  the  animal  to 
allow  entry  of  fresh  air.  Two  small  doors  (63  x  63  cm)  were  constructed  in 
each  side  of  the  chamber. 

A  block  diagram  of  the  exposure  apparatus  is  shown  in  Fig.  1.  Microwave 
power  was  obtained  from  an  Epsco  signal  source  operating  in  the  200-500  MHz 
band  with  a  minimum  output  of  100  watts  power.  A  Philco  model  164B  bidirec¬ 
tional  power  meter  was  used  to  monitor  power  delivered  to  and  reflected  from  the 
parallel  plate  chamber.  Power  leaving  the  parallel  plate  radiation  chamber  was 
monitored  by  a  Hewlett-Packard  No.  434A  calorimetric  power  meter.  A  Heathklt 
cantenna  was  used  for  final  dissipation  of  power  from  the  parallel  plate  cham¬ 
ber.  A  Welnschell  No.  109L  double  stub  tuner  was  used  to  balance  the  irradia¬ 
tion  system  and  reduce  reflected  power  to  a  minimum.  Frequency  of  the  micro- 
wave  energy  delivered  to  the  parallel  plate  chamber  was  monitored  by  a  Hewlett- 
Packard  No.  5253B  frequency  counter.  Interconnection  of  the  various  components 
of  the  Irradiation  system  utilized  UG58/U  coaxial  cable. 

Determination  of  Absorbed  Microwave  Power 

The  Intensity  of  the  microwave  field  within  the  parallel  plate  chamber  at 
the  position  of  the  rat  was  measured  as  mW/cm^  by  a  General  Microwave  Corpora¬ 
tion  Raham  model  1  (300-18,000  MHz)  field  intensity  probe. 

A  calorimetric  method  [2]  was  used  to  determine  how  much  of  the  microwave 
power  was  actually  absorbed  by  the  animal  during  each  exposure.  The  amount  of 
absorbed  power  was  expressed  as 

AT 

mW/gra  >■  4180  x  — 


where 


AT  =  rectal  temperature  increase  of  rat  in  degrees  Celsius  during 
microwave  exposure 


t  =  time  of  exposure  in  seconds 


1 


4180  =  converaioa  constant 


Procedure 

Each  rat  was  habituated  for  thirty  minutes  a  day  over  a  four-day  period  to 
a  cylindrical  Plexiglas  rat  holder  (.32  cm  x  6.5  cm  x  24  cm).  The  holder  was 
constructed  with  1.27  cm  slots  on  both  upper  and  lower  surfaces  to  interfere 
little  with  normal  heat  dissipation  by  the  animal  and  to  ensure  release  of  urine 
and  fecal  material.  On  exposure  days  a  microwave  transparent  liquid  crystal 
fiber  optic  temperature  probe  [7]  was  inserted  colonically  (4  cm)  and  held  in 
place  by  masking  tape  applied  to  the  rat's  tail.  Each  rat  was  then  placed  in 
the  rat  holder  and  set  inside  the  radiation  chamber  such  that  the  length  of  the 
rat  body  was  parallel  to  the  electric  component  of  the  microwave  field  (E||  L) . 
Each  rat  was  allowed  approximately  twenty  minutes  before  microwave  exposure  to 
allow  for  colonic  thermal  equilibrium.  Exposure  to  microwave  irradiation  com¬ 
menced  when  colonic  temperature  had  stabilized.  In  each  case  microwave  power 
was  gradually  Increased  over  a  5-second  period  to  the  established  field  density 
of  25  mW/cm^.  This  was  done  to  avoid  peak  surges  of  power  above  25  mW/cm2  at 
radiation  onset.  Each  rat  was  then  exposed  to  ten  minutes  of  microwave  radia¬ 
tion  at  one  of  twenty  frequencies  in  the  220-500  MHz  range.  The  microwave  field 
intensity,  mW/cm^,  and  amount  of  power  absorbed  by  the  animal,  mW/gm,  were 
recorded  for  each  exposure.  All  rats  were  exposed  randomly  to  the  twenty  fre¬ 
quencies.  This  procedure  was  repeated  four  times.  Some  of  the  rats  were  addi¬ 
tionally  exposed  to  microwave  radiation  at  360,  480,  and  500  MHz  with  the 
length  of  the  rat  body  parallel  to  the  magnetic  field  (h||  L)  and  the  direction 
of  microwave  propagation  through  the  chamber  (cl|  L).  Relative  humidity  during 
the  radiation  exposures  varied  between  12  percent  and  40  percent.  Ambient 
temperature  of  the  radiation  chamber  was  controlled  by  room  temperature  (70-72° 
F).  All  exposures  were  CW  radiation. 

EXPERIMENT  I;  RESULTS 

Figure  2  shows  the  mean  relative  absorption  of  microwave  power  ((mW/gm)/ 
(mW/cm^jj  and  a  scale  indicating  the  animal  rise  in  body  temperature  (centigrade) 
per  minute  from  220-500  MHz  for  animals  aligned  with  the  long  axis  of  the  body 
parallel  to  the  electric  field  (e||  L).  Also  shown  are  the  mean  ((mW/gm)/(mW/ 
cm2))  values  of  absorbed  power  for  animals  aligned  parallel  to  either  the  mag¬ 
netic  field  (h||  L)  or  with  the  direction  of  wave  propagation  through  the  cham¬ 
ber  (c||  L)  at  360,  480,  and  500  MHz.  The  most  microwave  power-absorbant 
animal  orientation  in  this  frequency  range  clearly  was  with  the  long  axis  of  the 
animal  parallel  to  the  electric  field  (e]1  L) .  In  this  orientation  500  MHz 
appears  to  be  the  most  absorbant  frequency.  Analysis  of  variance  techniques 
[8]  indicates  that  there  is  a  differential  absorption  of  power  in  the  electric 
orientation  depending  on  frequency  (£  <  .01).  Use  of  Newman  Kuels  multiple 
comparison  techniques  further  indicates  that  exposures  to  500  MHz  irradiation 
in  the  electric  orientation  were  clearly  different  from  all  other  exposure  fre¬ 
quencies  in  this  orientation  (g  <  .05)  except  360  and  440  MHz.  Exposure  of 
rats  to  360  MHz  radiation  in  this  orientation  was  also  different  from  other 
frequency  exposures  (g  <  .05)  except  at  260,  440,  460,  480,  and  500  MHz.  These 
results  Indicate  that  for  the  passive  rat,  simple  body  heating  as  measured  by  a 
change  in  colonic  temperature  is' dependent  on  both  frequency  (g  <  .01)  and  body 
orientation  relative  to  the  microwave  fields  within  the  radiation  chamber. 


EXPERIMENT  II:  MATERIALS  AND  METHOD 


Subjects 


Six  male  Long-Evans  rats  350-380  grams  in  body  weight  were  used  as 
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Figure  2.  Mean  frequency  and  body  orientation  dependent  microwave 
absorption  and  colonic  temperature  rise  in  the  Long  Evans  rat. 
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subjects.  Al]  animals  were  deprived  of  food  to,  and  maintained  at,  85  percent 
of  free-feeding  body  weights. 


Apparatus 

The  parallel  plate  irradiation  system  and  rat  holder  were  identical  to 
that  of  Experiment  I.  The  rat  holder  was  modified  by  attaching  a  small  Plexi¬ 
glas  response  lever  and  small  1.25  cm  diameter  glass  food  cup.  Since  standard 
metallic  microswitches  cannot  be  used  in  a  microwave  field,  an  alternate 
system  of  monitoring  lever-pressing  behavior  was  devised.  This  was  accomplished 
by  passing  a  concentrated  light  beam  through  47  mil  plastic  Dupont  Crofon  fiber 
optics.  The  fiber  optics  were  arranged  in  such  a  way  that  each  lever  press 
made  by  the  rat  interrupted  the  passage  of  the  light  beam.  A  sensitive  photo¬ 
cell  arrangement  outside  of  the  radiation  chamber  could  then  detect  each  inter¬ 
ruption  of  the  light  beam  and  consequently  each  lever  press.  A  Lehigh  Valley 
Electronics  Corporation  pellet  feeder  was  mounted  outside  the  radiation  cham¬ 
ber  and  delivered  45  mg  Noyes  food  pellets  to  the  glass  cup  of  the  rat  holder 
inside  the  radiation  chamber  via  a  1.27  cm  x  1  m  section  of  Tygon  flexible 
plastic  tubing.  Scheduling  and  recording  of  behavioral  test  sessions  was  ac¬ 
complished  via  a  system  of  relay  circuitry,  digital  printing  counters,  and 
cumulative  response  recorders. 


Procedure 

All  six  rats  were  trained  to  lever  press  for  food  pellets  on  a  random 
interval  thirty-second  schedule  of  reinforcement.  The  clock  time  used  for  the 
scheduling  was  three  seconds,  with  a  probability  of  .10  food-pellet  delivery. 
With  this  arrangement  a  food  pellet  could  be  delivered  from  three  seconds  to 
every  180  seconds,  with  Che  average  time  of  each  food-pellet  delivery  being 
every  thirty  seconds.  Each  rat  was  given  daily  thirty-minute  lever-pressing 
sessions  until  a  stable  rate  of  responding  was  achieved.  Lever-pressing 
behavior  was  considered  stable  when  the  number  of  lever  presses  during  the 
first  five  minutes  of  a  session  varied  less  than  10  percent  from  a  similar 
measure  of  the  previous  daily  session.  On  the  average,  approximately  twenty 
dally  lever-pressing  sessions  were  required  to  achieve  stable  response  rates. 

Once  stable  responding  rates  were  achieved,  all  animals  were  repeatedly 
and  randomly  exposed  to  different  conditions  of  microwave  radiation.  For  all 
microwave  exposures,  Che  field  density  within  the  chamber  at  the  position  of 
the  rat  was  maintained  at  25  mW/cra2  as  measured  by  the  field  intensity  meter. 

All  rats  were  exposed  randomly  to  360,  480,  and  500  MHz  microwave  radiation, 
with  both  the  length  of  the  rat  body  parallel  to  the  electric  field  (e||  L)  and 
direction  of  wave  propagation  through  the  microwave  chamber  (c||  L) .  In  this 
fashion  each  rat  was  exposed  to  the  six  possible  combinations  of  three  frequen¬ 
cies  and  two  body  orientations  of  microwave  exposure  while  performing  a  stable 
lever-pressing  Cask.  Radiation  exposures  with  the  rat  body  parallel  to  the 
magnetic  field  were  not  done  due  to  the  difficulty  in  rotating  the  radiation 
chamber.  Retraining  sessions  were  given  each  rat  after  each  microwave  exposure 
to  maintain  stable  rates  of  responding  before  the  next  microwave  exposure  and  to 
check  for  possible  carryover  effects.  Microwave  irradiation  began  at  the  begin¬ 
ning  of  the  sixth  minute  of  each  exposure  session  and  terminated  at  the  end  of 
the  first  one-minute  Interval  during  which  responding  fell  below  one-third  of 
the  rat's  normal  response  rate  per  minute  determined  from  the  previous  daily 
session.  If  the  rat  responded  throughout  the  session  at  a  rate  above  one-third 
the  baseline,  irradiation  ceased  at  the  end  of  the  thirty-minute  session, 
resulting  in  a  maximum  exposure  of  25  minutes.  The  primary  measure  of  lever- 
pressing  disruption  was  time  to  work  stoppage  defined  as  the  time  from  radiation 
onset  until  the  end  of  the  minute  during  which  responding  fell  below  one-third 


id 


of  the  baseline  rate.  Each  rat  was  removeu  from  tie  radiation  cliamber  and 
animal  holder  at  work  stoppage  or  at  the  end  of  the  session. 

Colonic  temperature  measurements  were  made  at  tlie  beginniiig  of  the  radia¬ 
tion  sessions  and  again  at  work  stoppage  or  the  end  of  an  exposure  session  to 
determine  the  amount  of  microwave-induced  heating  for  each  of  the  exposure 
conditions.  Water  intake  during  the  24-hour  period  between  each  daily  session 
while  each  animal  was  in  its  home  cage  was  measured  and  recorded.  Relative 
humidity  during  the  radiation  sessions  varied  between  12  percent  and  40  per¬ 
cent.  Ambient  temperature  of  Che  radiation  chamber  was  controlled  by  room 
temperature  (70-72°  F) . 

Because  of  a  lack  of  closed  circuit  television  equipment,  visual  observa¬ 
tion  of  the  animal's  performance  during  microwave  radiation  was  not  possible. 

EXPERIMENT  II:  RESULTS 

Frequency  of  the  microwave  radiation  and  orientation  in  the  microwave 
field  both  proved  to  be  of  major  importance  in  determining  disruption  of  lever¬ 
pressing  behavior.  Figure  3  shows  the  time  to  work  stoppage  from  radiation  on¬ 
set  for  the  frequency  and  body  orientation  parameters.  For  microwave  exposure 
sessions  with  the  animal's  length  aligned  parallel  to  the  direction  of  wave 
propagation  through  the  chamber  (c||  L) ,  nearly  all  animals  lever  pressed  at 
near  normal  rates  throughout  the  session.  Only  one  animal,  and  this  at  500  MHz 
radiation  in  the  c||  L  orientation,'  reduced  responding  to  the  work-stoppage 
criterion.  With  the  length  of  the  rat  body  aligned  parallel  to  the  electric 
field  (e||  L),  however,  time  to  work  stoppage  was  much  shorter.  Average  values 
of  time  to  work  stoppage  and  standard  deviations  are  shown  in  Fig.  3.  Animals 
exposed  to  360  and  480  MHz  radiation  in  the  Ej]  L  orientation  appeared  nearly 
normal  upon  removal  from  the  microwave  chamber. 

Animals  exposed  to  500  MHz  radiation  in  the  e||  L  orientation  appeared 
flushed  and  remained  in  a  prone  position  with  the  body  extended.  The  condition 
of  the  animal  appeared  similar  to  the  description  of  what  Justesen  and  King  [9] 
termed  "flaccid  paralysis".  Four  of  the  animals  exposed  under  these  conditions 
had  urinated  profusely  in  the  animal  holder.  The  fur  on  the  underside  of  the 
animal’s  body  was  wet  in  spite  of  the  slots  made  in  the  holder  to  release  urine 

from  the  holder  and  away  from  the  microwave  fields.  The  rats  were  clearly  heat 

stressed.  Each  animal,  however,  quickly  regained  its  upright  position  within 
minutes  after  removal  from  the  radiation  chamber.  A  randomized  block  factorial 
analysis  of  variance  [8]  was  used  to  evaluate  the  effects  of  the  treatments. 
This  test  indicated  a  significant  difference  between  the  frequency  of  radiation 
during  exposure  (F  =  17.68,  df  =  2/25,  £  <  .001),  the  animal  orientation  in  the 
microwave  field  (F  =  119.24,  d£  =  1/25,  £  <  .001),  and  a  frequency  by  orienta¬ 
tion  interaction  (F  =  12.37,  df  =  2/25,  p  <  .001). 

A  Newman  Kuels  test  found  no  difference  in  the  times  to  work  stoppage  for 
the  exposures  with  the  length  of  the  animal  aligned  parallel  to  the  direction 
of  the  microwave  propagation  througii  the  chamber  (cj]  L)  at  360,  480,  and  500 
MHz.  Each  of  these  exposure  conditions  was,  however,  different  (p  <  .05  or 
greater)  from  the  exposures  with  the  animal  aligned  with  the  length  of  the  body 
parallel  to  the  electric  field  (e|(  L)  at  360,  480,  and  500  MHz.  With  the 
animal  aligned  parallel  to  the  electric  field,  times  to  work  stoppage  at  360 
and  480  MHz  were  not  significantly  different.  Exposure  to  500  MHz  radiation  in 
this  orientation  was,  however,  different  (£  <  .01)  from  all  other  exposure 
conditions. 

The  amounts  of  body  heating,  as  measured  by  colonic  temperature  change 
during  irradiation,  paralleled  the  behavioral  effects.  Figure  4  shows  the 
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relative  abKorptioii  ot  microwave  energy  as  determined  calorlmetrica] ]y  and 
given  as  mW/gra.  Also  shown  in  Fig.  4  is  a  scale  showing  the  animal  body 
temperature  fisc  (centigrade  per  minute).  For  animals  aligned  parallel  to  the 
wave  propagation  (c||l),  body  heating  was  minimal.  For  animals  aligned  paral¬ 
lel  to  the  electric  field  (e]|l),  however,  heating  was  significant  and  in¬ 
creased  as  the  frequency  of  radiation  increased.  A  randomized  block  factorial 
analysis  of  variance  [8]  was  used  to  evaluate  the  effect  of  the  treatments  on 
rat  whole  body  heating  as  measured  rectally.  This  test  indicated  a  difference 
between  the  frequency  of  radiation  during  exposure  (F  =  21.70,  df  =  2/25,  p  < 
.01),  the  animal's  orientation  in  the  microwave  field  (F  =  129. 8i,  df  =  1/25, 
g  <  .01),  and  a  frequency  of  radiation  by  animal  orientation  Interaction  (F  = 
8.92,  df  =  2/25,  g  -01). 

A  Newman  Fuels  test  indicated  that  rat  body  heating  during  exposure  was 
not  different  at  360,  480,  or  500  MHz  when  the  rat  body  was  parallel  to  direc¬ 
tion  of  wave  propagation  (c|l  L)  through  the  radiation  chamber.  With  the  rat 
aligned  parallel  to  the  electric  field  (Ej|  L) ,  however,  500  MHz  exposures  were 
different  (p  <  .01)  from  360  and  480  MHz  and  also  from  the  exposures  in  the 
wave  propagation  orientation  at  all  frequencies.  In  the  electric  field  orienta¬ 
tion,  rat  body  heating  at  360  and  480  MHz  were  not  different  from  each  other 
but  were  clearly  different  from  exposures  in  the  wave  of  propagation  orienta¬ 
tion  (g  <  .05). 

The  measures  of  time  to  work  stoppage  and  rat  body  heating  expressed  as 
mW/gm  followed  each  other  very  closely  in  a  reverse  fashion.  As  body  heating 
Increased,  time  to  work  stoppage  became  much  shorter.  The  correlation  coef¬ 
ficient  between  the  time  to  work  stoppage  and  mW/gm  measures  for  all  animal 
radiation  exposures  was  -.79  (g  <  .01). 

Individual  records  of  response  rate  are  shown  in  Figs.  5  and  6.  Respond¬ 
ing  during  radiation  in  the  e|1  L  orientation  at  500  MHz  is  shown  in  Fig.  5  with 
comparable  baseline  performance  and  recovery  responding  the  day  following  radia¬ 
tion  exposure  shown.  Responding  for  the  same  animal  during  radiation  in  the 
Cjl  L  orientation  at  500  MHz  is  shown  in  Fig.  6.  In  both  figures  radiation  on¬ 
set  is  shown  (RF^jj)  as  well  as  time  to  work  stoppage  and  radiation  termination 
(RFof f)* 

An  examination  of  the  average  number  of  responses  per  minute  during  micro- 
wave  irradiation  revealed  an  Interesting  finding.  A  comparison  of  each  animal's 
average  response  rate  during  microwave  irradiation  until  time  to  work  stoppage 
or  the  end  of  tVie  exposure  session  with  a  similar  measure  from  the  previous 
day's  baseline  performance  showed  no  difference.  Figure  7  shows  this  finding. 
Total  lever-press  rates  for  all  exposure  conditions  were  transformed  into  sup¬ 
pression  of  responding  ratios.  The  suppression  ratios  were  used  as  an  index  of 
the  effect  of  microwave  radiation  exposure  on  responding.  Each  suppression 
ratio  was  indexed  by  the  ratio  of  B/A  +  B  where  A  represents  the  total  number 
of  lever  presses  during  baseline  performance.  The  score  B  represents  the  total 
number  of  lever  presses  during  the  period  in  radiation.  Thus  a  ratio  of  0.50 
would  indicate  no  relative  change  of  lever  presses  during  microwave  irradiation 
compared  to  the  same  measure  from  baseline  responding  with  no  microwave  irradi¬ 
ation  present.  A  ratio  of  0.00  would  indicate  a  complete  cessation  of  lever 
pressing  during  microw’ave  radiation  exposure.  At  each  of  the  microwave  exposure 
conditions,  responding  remained  fairly  stable  with  an  abrupt  drop  or  cessation 
of  responding  d\i’'lng  the  one-minute  interval  when  the  work  stoppage  criterion 
was  met.  A  rando.nized  block  factorial  analysis  of  variance  [8]  revealed  that 
suppression  of  responding  up  until  the  work  stoppage  criterion  was  not  different 
for  either  the  frequency  of  microwave  radiation  (p  >  .10)  or  the  orientation  of 
the  animal  in  the  microwave  field  (g  >  .10). 


Similarly  a  measure  of  the  amount  of  water  consumed  by  the  rat  between 


Figure  5.  Cumulative  response  record  of  rat 
Mlv’-20,  500  MHz,  E||  L,  radiation  exposure. 
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Figure  7.  Mean  suppression  ratios  of  responding  during  microwave  exposure 


daily  lever-pressing  sessions  was  also  not  iniiuenced  by  exposure  to  microwave 
radiation.  The  amount  of  water  consumed  by  each  rat  after  microwave  exposure 
until  the  next  daily  session  was  measured.  This  amount  of  water  intake  (cc) 
was  then  compared  to  a  similar  measure  of  water  intake  before  microwave  expo¬ 
sure.  Figure  8  shows  the  mean  amount  and  standard  deviations  of  water  consumed 
after  microwave  radiation  exposure  as  a  percentage  of  each  animal's  baseline 
■jater  intake.  A  randomized  block  factorial  analysis  of  variance  [8]  revealed 
no  significant  difference  in  water  intake  as  a  function  of  frequency  of  micro- 
wave  radiation  (£  >  .10)  and  orientation  of  the  animal  in  the  microwave  field 
(£  >  .10). 

To  determine  whether  the  effects  of  a  particular  radiation  session  iiad  any 
effect  on  subsequent  responding  in  recovery  test  sessions,  each  animal's  per¬ 
formance  during  the  first  five  minutes  of  the  baseline  session  was  compared  to 
performance  on  the  day  after  or  two  days  after  microwave  radiation  exposure. 
Standard  suppression  ratios  of  B/A  +  B  were  used.  The  animal's  responding  dur¬ 
ing  the  first  five  minutes  of  the  baseline  sessions  is  represented  by  the  score 
A,  whereas  total  responding  during  the  first  five  minutes  of  the  recovery  ses¬ 
sion,  either  one  or  two  days  after  microwave  irradiation,  is  represented  as  the 
score  B.  Thus  a  ratio  of  .50  would  indicate  no  relative  change  in  responding 
during  the  recovery  sessions  compared  to  performance  during  previous  baseline 
measures.  Figure  9  shows  the  result  of  this  comparison,  suggesting  little  car¬ 
ryover  effect  from  radiation  sessions  into  recovery  sessions  one  or  two  days 
later.  A  randomized  block  factorial  analysis  of  variance  [8]  revealed  that 
there  were  in  fact  no  significant  differences  in  the  two  performance  measures 
due  to  either  the  frequency  of  radiation  factor  (£  >  .10)  or  the  orientation  of 
the  animal  in  the  microwave  field  (g  >  .10). 

DISCUSSION 

This  study  was  conducted  to  answer  two  questions.  First,  does  the  fre¬ 
quency  of  microwave  radiation  and  orientation  of  an  animal  in  a  microwave  field 
affect  tlie  absorption  of  microwave  energy  and  consequent  animal  body  heating? 
Second,  tan  these  factors,  by  producing  differential  heating,  significantly 
alter  the  behavior  of  an  animal  performing  in  a  microwave  field?  The  data 
reported  here  support  an  affirmative  answer  to  both  of  these  questions. 

The  results  of  both  experiments  clearly  indicate  a  strong  effect  of  micro- 
wave  frequency  and  animal  body  orientation  in  determining  microwave  absorption, 
rat  body  leating,  and  consequent  suppression  of  behavioral  responding.  Animals 
exposed  S',  the  most  absorbant  frequency  and  body  orientation  combination  (500 
MHz,  e||  L,  were  clearly  heat  stressed.  Vasodilation  and  assumption  of  a  prone 
body  position  in  these  animals  after  radiation  exposure  were  marked.  Time  to 
work  stoppage  for  animals  exposed  under  these  conditions  was  short.  For  animals 
exposed  to  less  absorbant  frequencies  and  body  orientations,  the  rate  of  body 
heating  was  less  and  consequent  disruption  of  behavior  was  not  evident.  The 
high  correlation  between  rate  of  body  heating  and  time  to  work  stoppage  further 
supports  the  hypothesis  that  differential  microwave  absorption,  body  heating, 
and  behavioral  suppression  can  be  governed  by  the  factors  of  microwave  fre¬ 
quency  and  animal  body  orientation  in  the  microwave  fields. 

The  question  may  be  asked  as  to  how  mild  hyperthermia  may  be  respcnsible 
for  suppression  of  a  well-trained  bp>iavior  in  the  rat.  It  is  interesting  to 
note  that  the  lever-pressing  animal  maintained  near  normal  average  response 
rates  during  microwave  irradiation.  There  seemed  to  he  an  abrupt  cessation 
of  lever  pressing  during  microwave  exposure,  under  niost  absorbant  conditions, 
rather  than  a  gradual  decline  in  behavior  below  the  work-stoppage  criterion. 

This  effect  was  unexpected,  but  in  agreement  with  the  eoriier  report  by 
Justesen  and  King  [9]  of  a  similar  finding  in  rats  exposed  to  microwave 
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Figure  9.  Carryover  effects  suppression  ratios. 


radiation.  Roberts,  et  at.  ,  [10]  have  shown  that  rats,  heated  by  raising  the 
environmental  temperature,  .show  characteristic  changes  in  behavioral  thermo¬ 
regulation.  Rats  heated  in  this  fashion  show  an  initial  increased  level  of 
body  grooming  and  activity  during  the  first  ten  minutes  of  exposure  to  ambient 
temperatures  of  37.7*  C.  .Signs  of  body  relaxation  and  extension,  characteris¬ 
tic  of  the  heat-stressed  rat,  occurred  only  after  45  minutes  exposure  to  37.7®C 
ambient  temperature.  The  possibility  then  exists  that  efferent  mechanisms 
responsible  for  grooming  behaviors  as  a  thermoregulatory  response  may  also  be 
elicited  in  the  microwave  heated  rat.  Mild  hyperthermia  Induced  by  microwave 
exposure  may  elicit  thermoregulatory  behavior  in  the  rat  such  as  grooming  which 
conflict  with  lever-pressing  behavior.  Thus  time  to  work  stoppage  measured  in 
this  experiment  may  in  fact,  for  some  animals  at  least,  be  measuring  time  to 
elicitation  of  grooming  behavior.  Since  microwave  exposure  may  be  a  more  ef¬ 
ficient  and  rapid  method  of  Inducing  mild  hyperthermia,  several  classes  of  the 
classic  thermoregulatory  behaviors  may  be  elicited  simultaneously  or  in  rapid 
succession.  This  latter  possibility  may  be  true  for  the  animals  exposed  under 
the  most  absorbant  radiation  conditions.  Since  direct  observation  of  each 
animal  during  Irradiation  could  not  be  done,  further  Investigation  of  this 
hypothesis  Is  required. 

It  is  possible  that  other  factors  may  be  responsible  for  behavioral  sup¬ 
pression  than  just  mild  hyperthermia  induction  and  elicitation  of  behavioral 
thermoregulation.  One  may  speculate  that  evaporative  water  loss  during  micro- 
wave  exposure  was  enhanced  and  chat  consequent  thirst  induction  may  contribute 
to  behavioral  suppression.  A  finding  of  no  change  in  normal  water  intake  over 
the  24-hour  period  after  microwave  exposure,  however,  seems  contrary  to  this 
hypothesis.  The  lack  of  any  evidence  of  proactive  effects  due  to  microwave 
exposure  suggests  that  the  most  reliable  effect  of  such  exposure  was  body  heat¬ 
ing  sufficient  to  suppress  behavior  but  not  intense  enough  to  produce  actual 
physical  thermogenic  damage  Co  Che  animal.  This  may  suggest  Chat  the  primacy 
effect  of  moderate  intensity  microwave  exposure  is  to  induce  mild  hyperthermia 
and  consequent  suppression  of  behavior,  both  of  which  are  transient  in  nature. 

It  is  important  to  note  in  Experiment  I  that  even  though  each  animal  was 
allowed  sufficient  time  for  colonic  thermal  equilibration  before  exposure  to 
microwave  radiation,  each  animal  was  alert  during  the  exposure  period.  It  is 
well  known  that  emotional  reactions  produce  psychogenic  effects  in  the  rodent 
and  that  temperature  regulation  can  be  altered  by  such  factors.  It  is  difficult 
to  rule  out  such  reactions  to  the  environment  when  using  colonic  temperature  in 
the  alert  animal  as  the  dependent  variable.  The  use  of  habituation  techniques 
such  as  familiarizing  the  animal  to  the  animal  holder  and  radiation  chamber  may 
help  minimize  such  reactions.  Nevertheless  the  animal  cannot  be  readily 
habituated  to  the  microwave  exposure  without  long-term  adaptation  to  induced 
heating  effects.  Determination  of  absorbed  microwave  power  as  described  in 
this  report  must  be  distinguished  as  a  dynamic  measure,  since  the  rat  is  cap¬ 
able  of  active  thermoregulation  during  irradiation. 
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I.  INTRODUCTION 

In  electromagnetics,  discretization  for  transforming  an  inte¬ 
gral  equation  to  a  matrix  equation  is  commonly  accomplished 
using  pulse  functions  as  a  basis  [  1  |-|4I .  More  elaborate  bases 
such  as  modal  fields  may  be  justified  in  a  particular  problem 
but  are  avoided  in  the  matrix  formulation  for  the  general 
scattering  body. 

Discretization  with  pulse  functions  requires  that  the  scatter¬ 
ing  body  be  partitioned  into  a  number  N  of  cells,  where  N  is 
large  enough  that  complex  permittivity  and  the  complex  time- 
independent  electric  field  may  be  assumed  constant  within 
each  subvolume.  Procedures  for  improving  convergence  with 
the  general  scattering  body  are  of  interest  since  the  cost  of 
computation  is  a  rapidly  increasing  function  of  N. 

This  communication  describes  two  new  procedures  which 
have  been  found  to  improve  convergence  in  solution  of  the 
two-dimensional  problem  of  transverse  magnetic  (TM)  exci¬ 
tation  of  an  infinite  cylinder  of  arbitrary  cross-section  shape 
Convergence  is  improved  because  the  calculations  include  part 
of  the  variation  of  the  field  within  each  cell. 


II.  REVIEW  OF  EARLIER  SOLUTIONS  USING 
PULSE  FUNCTIONS 

For  TM  excitation  of  an  infinite  cylinder,  both  incident  and 
scattered  fields  have  E  directed  parallel  to  the  infinite  dimen¬ 
sion  (z  axis).  The  electric  field  is  described  by  a  scalar  integral 
equation: 

£z(x,  y)  =  y)  -  ~  j  j  lfr(A'.  -  1 1 

'  £■/(•* ’.  y')Mo^'‘HkoP)dx'  dy\  (1  I 
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where  (x,  y)  and  (x',  y')  are  coordinates  of  the  observation 
point  and  source  point,  respectively;/?  is  the  distance  between 
the  two  points;  E/  and  E^  are  incident  and  total  electric  field 
intensity,  respectively,  with  time  dependence;  e^fx',  y')  is 
complex  permittivity  at  the  source  point  relative  to  free  space; 
*0  =  snd  Hfy^^HkoP)  is  the  Hankel  function  of  zero 

order. 

The  discrete  analog  of  (1)  consists  of  the  A'-by-A' system  of 
linear  equations; 

N 

tw  =  l,2,"*,  Af.  (2) 

n  =  1 

Richmond  performed  the  discretization  using  pulse  func¬ 
tions  [  1 J .  Each  square  cell  was  approximated  by  a  circle  of 
radius  a  having  an  equal  area  to  simplify  integration  of  the 
Bessel  functions.  His  expressions  for  the  matrix  elements 
follow: 

(3) 

=  '— ;^(ern  -  >)^l(*0«Vfo‘*H*oPmn).  m  n. 

(4) 


=  I  +  (Cr, 


1) 


inkna 
1  + - 
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hig.  1.  Nolalion  used  in  evaluation  of  off-diagonal  mainx  clemenls 


In  this  ."nethod,  the  electric  field  intensity  is  considered  to 
be  constant  over  the  area  of  a  cell,  an  approximation  that 
results  in  significant  error  if  the  cells  are  not  very  small.  Larger 
cells  can  be  used  if  the  variation  of  the  electric  field  is  included 
in  the  calculations.  A  method  of  approximating  the  variation 
of  the  electric  field,  called  "plane-wave  correction,”  is  described 
in  the  next  section. 

III.  PLANE-WAVE  CORRECTION 

Variation  of  is  initially  unknown,  but  we  may  make  the 
approximation  of  representing  the  fields  within  each  cell  by  a 

superposition  of  plane  waves.  AU  member  plane  waves  are  /.  a,  ,i,e  center  of  the  mth  cell,  as  shown  in  Fig.  1 ,  Lei  p„„ 

required  to  have  propagation  vectors  of  magnitude  k  corre-  g^d  r  be  distances  from  the  center  of  the  nth  cell  to  the 

spending  to  the  frequency  of  the  incident  wave  and  local  observation  point  and  source  point,  respectively.  If  the  angle 

complex  permittivity  of  the  dielectric  body.  If  Oj  and  repre-  between  p^n  and  r  is  6  ,  then  the  distance  between  the  source 

sent  orientation  of  propagation  and  amplitude  of  the  ith  com-  pojn,  and  observation  pomt  is 

ponent,  respectively,  the  superposition  of  piece  waves  may  be 
written  as 


Pmu^  "  yPmn  COS  7’. 


CM 


E^(r'.  C')  =  2 


(5) 


If  e,  is  constant  within  the  nth  cell  and  the  approximation  of 
circular  cells  is  used  in  (1),  then  the  off-diagonal  matrix 
If  complex  permittivity  is  constant  within  a  cell  and  the  elements  are 
approximation  of  circular  cells  is  us' in  1 1 ),  then  the  diagonal 
matrix  elements  are 


Moslem  - 


=  1  + 


4£,(0, 

•  2tr 
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(6) 


Using  (5)  for  E^ir  .  0  ncglecMr.g  terms  higher  than  second 
order  .‘n  kr'  in  the  expansion  of  the  exponential,  and  perform¬ 
ing  the  O'  integrations. 


=  1  -I-  -- 


4  2 


4£'^(0,  O') 

•  I  EAr  .0  . 

Jo 

Using  Graf’s  addition  theorem  1 5 1  in  ( 1 0), 

4  =  _ _ _ 

^tnn  , 

4/;,(0.0  ) 


(101 


1(2/2 


fJo^^’ikor')r  dr'. 


f  Jiikgr'V' dr'  f  E;{r'.  O' 
Jo  Jo 


)  cos  10' dO'.  Ill) 


(7) 


Performing  the  integrations  over  the  Hankel  functions  gives, 
for  the  diagonal  elements. 


Using  (5)  for£^(r’.  O')  and  neglecting  terms  higher  than  second 
order  in  kr'  in  the  expansion  of  the  exponential. 


^mm  ^  i^rm  ~  0 
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cos  10'  dO'. 


In  evaluating  the  off-diagonal  matrix  elements,  suppose  that 
the  source  point  P'  is  in  the  nth  cell  and  the  observation  point 


(12) 
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Performing  the  d'  integrations  in  (12),  and  neglecting  terms  If  quadratic  and  higher  order  terms  in  (6)  are  neglccieJ.  then 
with  I  /  I  >  1 , 


ikoHern  -  I)Wo<*>(A.oPmn) 


ka 


t{a.  0)  -  t{{).  0) 
ti0,0) 
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Performing  the  integrations  in  ( 14), 
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•  j  B,rr*r'  cos  a, 


Since  only  the  last  membei  in  il'i  leper  :>  i  ■ 
unknown  coefficients  in  the  plane  wave  expansi.  i.  ai  ■ 
delete  the  term  and  use  the  first  three  niemt  er-  :  •  r  . 
diagonal  matrix  elements. 


jnkoa 


Urn  I  I 


Jiikor'ydr'.  (13) 


^rn  \  ,  , 

1  {kf^Q)  I  J  ^ik^tj ) 

4  ' 


The  summations  in  all  three  integrands  of  (13)  contain 
factors  independent  of  the  summation  index  that  may  be 

removed  from  the  summations.  Also  changing  to  a  dummy  The  use  of  (18)  is  justified  only  if  the  la^t  member  of  ( 1  7)  is 
variable,  small.  The  fractional  error  in  A „„  due  to  neglecting  the  last 

_ _ _  — _  member  of  ( 1 7)  is 


I  2(^0^) 


Ho^^KkoPmn) 


JlUoo)  ~  l^o<^J2Uo‘‘) 


tl'M 
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rkQO 
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Since  the  magnitude  of  the  ratio  of  the  Hankel  functions  is 
greatest  for  a  small  argument,  we  set  p„„  =  s/na  for  adjacent 
cells.  Expanding  the  Bessel  functions  for  small  arguments  and 
using  ( 1 6)  to  estimate  7,  we  obtain 


8 


+  -Urr,  -  iW^nUl^^HkoPmn) 


2  Bj  cos  a, 


2^1 


I6ma  l<- 
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4\/7r  In 
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£(fl,  0)  -  £(0.  0) 
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(20) 


Define 


rkoa 

I  Ji(x)x^  dx. 
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,  Bi  cos  Oj 

/ 


where  F  =»=  1 .78107. 

In  the  plane  wave  correction  method,  (8)  and  (18)  are  used 
for  the  matrix  elements,  and  the  solution  corresponds  to  the 
(14)  value  of  at  each  cell  center.  Numerical  considerations 
suggest  that  deletion  of  the  last  member  in  (1 7)  is  the  principal 
source  of  error.  Equation  (20)  requires  that  if  4\/7r" 
In  (2/(r\/ffko‘*ll  exceeds  the  condition  number  of  the  matrix, 
then  the  error  in  E^  calculated  for  the  cell  centers  will  be  less 
than  the  variation  of  E^  within  each  cell  relative  to  values  at 
the  cell  centers.  Hence  the  plane  wave  correction  will  account 
for  much  of  the  variation  of  E^  within  each  cell.  We  must  have 
a  <  0.0876\o  fhe  criterion  will  not  be  satisfied  even  for  a 


(15) 
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condition  number  of  unity.  Similar  convergence  criteria  are 
not  available  in  the  method  using  pulse  functions  that  is 
described  in  Section  11.  A  second  method  of  approximating 
the  variation  of  the  electric  field,  called  "cylindrical-cell 
correction."  is  described  in  the  next  section. 


IV  CYLINDRICAL-CELL  CORRECTION 


In  the  methods  described  in  Sections  II  and  III,  the  scatter¬ 
ing  body  IS  divided  into  cells  that  are  approximated  by  circular 
cylinders  in  calculation  of  the  matrix  elements.  The  method  of 
cylindncal-cell  correction,  which  is  presented  for  the  first  time 
in  this  communication,  emphasizes  the  properties  of  such  a 
model  of  the  original  body. 

Let  the  cross  section  of  the  scatterer  be  divided  into  square 
cells,  each  having  side  i.  For  the  model,  replace  each  cell  with 
a  cucular  cylinder  having  radius  u  and  assume  properties  of 
free  space  between  the  cylinders.  For  the  chosen  geometry  (6) , 
the  effective  relative  complex  permittivity  of  the  model 
and  the  complex  permittivity  e/  of  each  cylinder  are  related  by 


«reff  -  • 


1 1. 


(21  1 


11  we  require  that  the  effective  relative  complex  permittivity 
of  the  model  equal  that  of  the  original  scatterer,  then 


ITtl 


used  to  calculate  the  diagonal  matrix  elemenis: 

i'”  . 

m  “  “  [  ^  irT  j  ( A:  u  )//y  ^  ^  ^  j 

-  .  (  ;(,i 

From  (10)  and  (24),  the  off-diagonal  matrix  elements  are  given 
by 
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Using  Grafs  addition  theorem  [.‘'1  and  the  notation  of  Fig.  1. 
jk  ^  * 

"'Imn  ~  (trn  "  ^  1  n  ^ 

4  1=  » 

•  f  Ji{kQr'\lQ(k‘r')r'  Jr' 

■>0 

Performing  the  6'  integration,  which  collapses  the  summation, 
and  using  a  relationship  from  [7]  for  the  /'  integration,  wc 
obtain  the  following  expression  used  to  calculate  the  ofl- 
diagonal  matrix  elements: 


/ 


cos  !fl'  dO 


I  2S  I 


If  the  electrical  size  of  the  cylinders  is  sufficiently  small, 
the  inhomogeneities  of  the  model  will  have  no  significant 
effect,  so  that  scattering  from  the  model  satisfying  (22)  will 
duplicate  that  of  the  physical  scatterer.  The  model  is  useful 
since  we  may  approximate  the  variation  of  the  electric  field 
within  each  cylinder.  A  circular  cylinder  with  TM  excitation 
will  have  internal  fields  given  by 

=  2  I’n-lnil^’')  cos  (nd  +  Cn),  (23) 

n  =  0 


rimn  =-  HQ^^HkoPmn)[k'aJi(k'a)Jo(koa) 

—  kQaJo{k'a)Ji{koa)\ .  m  ^  n. 

In  the  cylindrical-cell  correction  method,  (26)  and  (2'^)  are 
used  for  the  matrix  elements,  and  the  solution  corresponds  to 
the  value  of  at  each  cell  center.  A  solution  is  defined  for 
any  ratio  of  a/s.  but  only  two  cases  have  been  studied :  1)  tan¬ 
gent  cylindrical  cells 


where  the  b„  and  c„  are  determined  by  the  excitation.  If  the 
cylinder  is  sufficiently  small,  the  zero-order  term  will  domi¬ 
nate,  so  we  may  use  the  approximation 


£'2  =  (24) 

From  (6)  and  (24),  the  diagonal  matrix  elements  are  given 


=  1  + 


1) 


/■ 

Jo 
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s  ,  4 

“  =  -  ,  €r  =  1  +  -  (Cr  -  1  ) 

2  TT 


and  2)  overlapping  cylindrical  cells 


s 


For  very  large  or  very  small  values  of  a/s.  the  structure  of  the 
model  differs  significantly  from  that  of  the  scatterer,  so  that 
the  solution,  which  actually  corresponds  to  scattering  from  the 
composite  of  circular  cylinders,  will  differ  from  that  for  the 
desired  scatterer.  Advantages  of  the  method  will  be  seen  in  the 
example  in  Section  VI. 

V.  CALCULATION  OF  POWER  ABSORPTION 


where  k'  =->^er„  *0. 

Performing  the  d'  integration  and  using  a  relationship  from 
[7]  for  the  r'  integration,  we  obtain  the  following  expression 


With  both  plane-wave  correction  and  cylindncal-cell 
correction,  the  solution  corresponds  to  the  value  of  at  each 
cell  center.  The  plane-wave  correction  leaves  the  variation  of 
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/\.  within  each  cell  undefined,  but  with  cylindrical-cell 
correction,  (24)  describes  the  approximate  variation  of 
within  each  circular  cylinder  of  the  model. 

If  the  value  of  power  absorption  is  desired,  then  1/2 
may  he  computed  w  ithin  each  cell  using  the  values  of  E^  found 
using  pulse  functions  or  plane-wave  corrections.  In  this  section 
the  approximate  variation  of  E^  within  each  cell  is  used  for  an 
improved  expression  for  power  absorption  if  cylindrical-cell 
correction  is  used. 

Total  power  absorption  in  one  circular  cylinder  is  given  by 


E,E,*rilrM. 


t30) 


From  t  221  the  conductivity  of  the  cylinder  is  related  to  that  of 
the  scatterer  by 


(31) 


Substituting  ( 24)  and  ( 3  1 )  in  1 30)  gives 


~  ^(1^1(1*  /  J r)r  dt\ 

a-  Jo 


(32) 


I'sing  a  relationship  from  [  7 1  for  the  /  integration  and  simpli- 
lying.  the  average  absorbed  power  density  within  the  cell  is 
given  by 


I-..  I-./ 


Im  (/;'(/y„*(A'ul/i(A''a)j 
Re  ik'a)  Im  (k'a) 


(33) 


where  is  the  calculated  value  of  /•._ 
Using  cyluidrical-cell  correction. 


found  for  the  cell  center 


Fig.  2.  Variation  of  average  specific  absorption  rate  with  frequency 
for  infinite  cylinder  model  of  man.  Incident  power  density  is 
I  niW,cm2. 


VI.  M  Ml  Kl(  Al.  FXAMPLl  :  TllF  INFINITE  CYLINDER 
MODEL  OF  MAN 

I  he  infinite  cylinder  has  been  suggested  for  use  as  a  model 
tor  body  extremities  or  the  chest  cavity  in  the  evaluation  of 
biological  hazards  from  electromagnetic  radiation  [81.  Approx¬ 
imation  of  an  extremity  or  the  whole  body  by  an  infinite 
cylinder  may  he  made  using  the  ratio  of  volume  to  length  for 
the  cross-sectional  area. 

I  he  model  considered  in  this  example  is  an  infinitely  long 
cylinder  with  square  cross  section  of  20  by  20  cm,  having  TM 
excitation  hy  a  plane  wave  incident  perpendicular  to  one  of 
the  lour  congruent  faces.  The  body  is  homogeneous  with 
values  of  permittivity  and  conductivity  that  are  typical  for 
muscle,  skin,  and  other  tissues  having  high  water  content  [dj . 

Itie  methods  described  in  this  communication  have  been 
used  to  calculate  average  specific  absorbed  power  for  the 
model  Suitahilily  of  the  model  can  only  he  justified  at  high 
trequencies  where  end  effects  may  be  neglected,  but  calcu¬ 
lations  have  been  made  over  a  wide  range  of  frequencies  to 
allow  comparison  of  convergence  with  the  different  methods. 

Fig  2  illustrates  the  freijuency  dependence  of  average 
specific  absorption  rate  (SAR)  for  the  model.  The  values 
obtained  using  the  accepted  method  of  pulse  functions  with 
100  cells  are  the  standard  for  comparison.  Values  obtained 
using  pulse  functions  with  81  cells  were  found  to  differ  from 


those  with  100  cells  by  less  than  1  percent  for  frequencies  up 
to  200  MHz.  Values  found  using  9  cells  with  pulse  functions, 
plane  wave  correction,  and  cylindrical-cell  correction  with 
tangent  and  overlapping  cells  are  also  given  in  the  figure.  For  9 
cells,  all  four  calculations  are  m  good  agreement  with  the 
standard  at  low  frequencies,  where  the  cells  have  small  elec¬ 
trical  size  Calculations  with  the  new  methods  appear  to  have 
comparable  accuracy  and  appear  to  offer  a  significant  improve¬ 
ment  in  convergence  when  compared  to  pulse  functions. 

The  criterion  u  <  0.0876Xo  suggests  that  the  plane  wave 
correction  with  9  cells  will  cease  to  be  useful  for  frequencies 
exceeding  700  MHz,  but  the  exact  limit  depends  on  the  con¬ 
dition  number  as  mentioned  in  Section  HI.  Increasing  error  is 
observed  for  frequencies  exceeding  200  MHz  which  is  -on- 
sistent  with  the  convergence  criteria. 

VII.  CONCLUSIONS 

Two  methods  have  been  found  which  improve  convergence 
in  solution  of  the  two-dimensional  problem  of  TM  excitation 
of  an  infinite  cylinder  of  arbitrary  cross-section  shape.  The 
extension  to  the  three-dimensional  problem  has  not  yet  been 
completed. 

Convergence  is  improved  because  the  calculations  include 
approximation  for  the  variation  of  the  field  within  each  cell. 
The  new  methods  give  expressions  for  the  matrix  elements 
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that  are  different  from  those  found  using  pulse  functions,  but 
values  of  the  matrix  elements  are  in  good  agreement  for  elec¬ 
trically  small  cells. 

The  matrix  elements  found  with  the  new  methods  do  not 
require  a  substantial  increase  in  computational  effort.  For  a 
homogeneous  scatterer,  if  equal  size  cells  are  used,  the  new 
methods  require  one-time  computation  of  several  additional 
Bessel  functions,  but  time  for  calculation  of  the  matrix 
elements  is  dominated  by  recalculation  of  the  zero-order 
Hankel  function  which  is  needed  with  or  without  the 
corrections  made  in  the  new  methods.  If  tiie  scatterer  is  not 
homogeneous  or  if  different  cell  sizes  are  used,  then  time  for 
calculation  of  each  matrix  element  is  approximately  doubled 
in  the  new  methods.  Time  spent  in  calculation  of  the  matrix 
elements  is  proportional  to  A'*,  whereas  time  spent  in  solving 
the  matrix  equation  is  proportional  to  A'^ .  so  for  large  numbers 
of  cells  there  is  still  no  significant  increase  in  computation 
time  with  the  new  methods. 
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Abstract — Two  InterpoUnLs  are  described  that  may  be  used  to  ccarect 
the  results  of  a  mooient-metbod  solutiou  using  pulse  fuDctioas  as  a  basis 
and  delta  functions  for  testing.  Tbe  interpolants  allow  for  some  of  the 
variatioo  of  tbe  fields  within  each  cell  and  thereby  increase  accuracy  and 
improve  convergence.  The  interpolants  are  usable  for  a  general  scatterer 
and  typically  retpiire  about  1  percent  of  the  cost  of  the  initial  numerical 
solution. 

I.  Introduction 

OMENT-METHOD  solutions  for  electromagnetic 
scattering  with  a  general  dielectric  object  often  use 
pulse  functions  as  a  basis  and  delta  functions  for  testing 
[l]-(5].  The  object  of  this  paper  is  to  show  that  an  inter- 
polant  may  be  used  to  allow  for  some  of  the  variation  of 
the  fields  within  each  cell  and  thereby  increase  accuracy 
and  improve  convergence.  The  computational  expense  of 
using  the  interpolant  is  typically  about  1  percent  of  the 
cost  of  the  initial  numerical  solution. 

A  numerical  solution  using  a  pulse  function  basis  re¬ 
sults  in  a  single  value  representing  £  within  each  cell.  The 
delta  functions  used  for  testing  enforce  the  integral  equa¬ 
tion  at  the  center  of  each  cell  so  that  the  calculated  E 
values  are  most  representative  of  the  cell  centers.  Experi¬ 
mental  tests  have  shown  that  the  error  in  E  calculated  for 
the  cell  centers  is  relatively  small  even  when  adjacent  cells 
base  values  which  may  differ  by  an  order  of  magnitude 
[0).  Since  the  values  of  E  may  be  assigned  to  points  in 
space,  interpolation  is  possible. 

The  present  study  has  been  restricted  to  piece-wise 
interpolation  in  which  the  value  of  the  corrected  specific 
absorption  rate  (SAR)  is  calculated  for  one  cell  at  a  time. 
In  this  paper,  two  interpolants  that  have  been  developed 
by  the  authors  are  described.  The  triquadratic  interpolant 
is  usable  when  cell  centroids  are  on  a  Cartesian  product 
mesh.  The  NEWSUD  interpolant  is  useful  for  problems  in 
which  the  cubical  cells  may  have  different  sizes  and 
general  locations. 
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11.  Te.st  of  a  One-Di.mensionai-  Interpolant 

Moment-method  solutions  with  any  subsectional  basis, 
such  as  pulse  functions,  require  that  the  scatterer  be 
approximated  by  a  composite  of  cells  [l]-(5].  Since  ana¬ 
lytical  solutions  for  such  composites  are  unknown  in 
three-dimensional  problems,  a  test  of  convergence  has 
been  made  for  a  one-dimensional  problem. 

Consider  the  one-dimensional  problem  of  a  plane  wave 
polarized  in  the  x  direction  incident  upon  a  dielectric  slab 
which  extends  from  z  =  0  to  z  =  a.  The  electric  field  is 
described  by  a  scalar  integial  equation; 

EJz)  =  E;(z)-^^f\(,(z')-l)EJz  )e  (1) 

where  z  and  z'  are  coordinates  of  the  observation  point 
and  source  point,  respectively,  p  is  the  distance  between 
the  two  points,  and  £,  are  the  incident  and  total 
electric-field  intensity,  respectively,  with  e'“'  time  depen¬ 
dence,  t,(z')  is  the  complex  permittivity  at  the  source  point 
relative  to  free  space,  and  kQ  =  uV 
The  discrete  analogue  of  (1)  consists  of  the  N-hy-N 
system  of  linear  equations: 

A 

"Z  E',„.  W=l,2 . A  (2) 

n  “  1 

where  the  dielectric  is  partitioned  into  ,V  cells  which  are 
thin  slabs,  and  n  and  m  are  indices  for  the  cells. 

The  matrix  elements  are  readily  evaluated  using  ( I )  with 
the  procedure  in  [2]: 

/l„„=l-E(€„„-IKl-e  (3) 

l)e  sin  j,  (4) 

where  A  is  the  thickness  of  each  cell  and  is  the 
distance  between  the  centers  of  the  mth  and  nth  cells. 

The  analytical  .solution  for  a  homogeneous  dielectric 
slab  IS  given  by 

£.(c) 

£'(0)[  2  \  t,  cos  \  k„{  a  ~  z)  +  2j  sin  \  t,  A„( «  -  z )  1 

2\  t,  cos  \  E,  A||<i  -Ey(  I  -E  t,)  sin  \  t, 
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Fig.  1.  Calculated  average  specific  absorption  rate  (or  a  2-cm  thick 
layer  of  muscle  at  2450  MHz.  Incident  power  density  is  1  MW/cm^, 
«ri  =  50,  0  =  2.2  mho/m. 


Fig.  I  gives  the  average  SAR  for  a  2-cm  thick  layer  of 
muscle  at  2450  MHz  for  three  different  computational 
procedures.  The  analytical  value  was  found  by  integration 
of  1/2  oE  E*  using  (5).  Numerical  solutions  were  ob¬ 
tained  using  (2)-(4),  Numerical  values  without  the  inter- 
polant  were  found  by  averaging  1  /2  a  E  E*  for  each  cell. 
Numerical  values  were  also  obtained  by  using  a  piece-wise 
quadratic  interpolant  [7)  to  obtain  an  expression  for  £ 
within  each  cell  and  then  integrating  the  expressions  for  a 
volume  average  of  1/2  <7  E  E*. 

It  is  readily  seen  in  Fig.  1  that  the  interpolant  causes  a 
significant  improvement  in  convergence  to  the  known 
analytical  solution.  We  may  relate  the  improvement  in 
convergence  to  the  well-known  preference  of  Simpson’s 
rule  to  the  trapezoidal  rule  for  numerical  quadrature. 
Simpson's  rule  requires  fitting  a  piecewise  quadratic 
through  the  data  points  and  generally  gives  a  better  ap¬ 
proximation  of  the  integrand  function  for  greater  ac¬ 
curacy  than  the  trapezoidal  rule. 

111.  Triqi  ADRATir  Interpoi.ant 

We  have  developed  a  triquadratic  interpolant  which 
may  be  used  when  the  cell  centroids  are  on  a  Cartesian 
product  mesh.  As  the  name  suggests,  quadratic  interpola¬ 
tion  is  used  parallel  to  each  of  the  three  Cartesian  axes. 

A  ,jtal  of  27  cells  is  used  in  the  stencil  for  each 
calculation.  Fig.  2  shows  the  location  of  the  26  cells 
surrounding  the  cell  in  which  the  correction  is  made.  The 
Cartesian  product  is  shown  for  X  =  -  1.  0.  1.  y'=  —  I,  0,  I. 
Z  =  -  1.  0,  1,  but  scaling  is  readily  used  to  adjust  to  a 
specified  cell  size. 

Quadratic  interpolation  may  be  used  with  univariate 
data  for  the  function  /-(.V  )  at  the  three  points  ,V  =  -  1, 0,  I 
by  the  rule 

£(.V)=  £(  -  l)y(.V-  l)-(-£(0)(l  -,V) 

•(l-F.V)-F£(l)y(A  +  l).  (6) 
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Fig.  2.  Location  of  cells  in  slencil  for  the  Inquadraiic  inierpolanl 

Triquadratic  interpolation  may  be  performed  on  the  sten¬ 
cil  shown  in  Fig.  2  by  using  products  fi^rmed  from  three 
univariate  interpolanis.  For  example,  the  contribution  of 
the  function  value  at  .\'=  -  1.  >'=  -  1.  Z  =  0  (cell  number 
2)  is 

£(-l.-l.0)y(.V-|)y()  l)(l  /)(l-l-/|. 

Triquadratic  interpolation  as  described  in  the  last  para¬ 
graph  is  only  usable  for  correcticins  in  a  cell  surrounded 
by  26  other  cells,  as  shown  in  Fig.  2,  It  is  possible  to  use  a 
single  interpolation  rule  based  on  the  full  27-cell  stencil  if 
£  values  are  estimated  for  iinc'ccupied  positions  by  means 
of  a  series  of  fill-in  rules.  We  have  used  the  following 
series  of  rules  in  which  the  £  values  are  filled  in  the  order 
of  increasing  distance  from  the  central  cell  by  averaging  E 
values  known  in  adjacent  cells. 

1)  Fill  in  the  six  closest  cells  (numbers  5.  II.  I.l.  15.  17, 
23).  If  any  are  unoccupied,  use  the  value  of  £  in  tlie 
central  cell. 

2)  Fill  in  the  twelve  next  closest  cells  (numbers  2,  4.  fi.  8. 
10.  12,  16.  18.  20.  22.  24,  26).  If  anv  are  unoccupied,  use 
one-half  the  sum  of  f  for  the  pair  from  the  six  closest  cells 
of  part  I  which  share  a  face  with  the  uiuKcupied  cell. 

3)  Fill  in  the  eight  corner  cells  (numbers  I,  3.  7,  9,  |9, 
21.  25.  27).  If  any  are  unoccupied,  use  ime  third  the  sum 
of  £  for  the  three  cells  sharing  a  face  with  the  uniKcupied 
cell 

Application  of  triquadratic  interpolation  to  the  general 
scatterer  requires  the  following  calculations  be  niadc  for 
each  cell.  First,  a  sieve  is  used  to  find  which  of  the  26 
surrounding  cells  are  occupievi  Next,  the  fill-in  rules  are 
used  as  required  to  complete  the  stencil.  Finally,  the 
triquadratic  interpolant  is  used  to  evaluate  the  integral  of 
1/2  o  E  E*  over  the  central  cell.  If  the  cells  have  different 
complex  permittivities,  both  the  fill-in  and  interpolatum 
rules  must  be  modified  by  suitably  multiplying  bv  the 
ratio  of  permittivities  so  that  the  interpolant  has  both  D 
normal  and  E  tangential  to  each  boundary  continuous. 
We  have  made  such  modifications  and  have  used  inter¬ 
polation  in  solutions  for  inhomogeneous  dielectric 
scatterers. 

Since  analytical  solutions  are  not  available  for  arrays  of 
cubes,  a  standard  for  evaluating  results  obtained  by 
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numerical  methods  is  not  easily  obtained,  but  one  conver¬ 
gence  test  has  been  made  for  the  tnquadratic  mterpolant. 
big.  3  gives  the  average  specific  absorption  rate  calculated 
without  interpolation  for  a  12-cm  cube  of  muscle  with 
plane  wave  at  I  MHz.  normally  incident  on  one  face. 
Linear  convergence  is  demonstrated  in  the  figure.  L.x- 
trapolation  using  the  values  for  216  and  ,‘il2  cells  gives  an 
S.AR  of  0.3275  x  10  W /kg  for  an  incident  power  density 
of  I  MW/cm^.  Table  I  gives  the  estimated  errors  in  SAR 
without  interpolation  and  with  the  triquadratic  interpolant 
found  by  compari.son  with  the  extrapolated  value  of  SAR. 
Results  in  Table  I  suggest  that  the  triquadratic  interpolant 
causes  a  significant  improvement  in  convergence. 

IV.  NEW3UD  Imerpoi  AM 

We  have  developed  the  NEWSUD  interpolant  which 
may  be  used  when  the  cells  have  different  sizes  and/or 
arbitrary  placement.  The  NEWSUD  interpolant  is  more 
general  than  the  triquadratic  interpolant  and  does  not 
require  fill-in  rules  but  is  limited  to  linear  rather  than 
quadratic  correction.  The  name  NliWSUD  is  an  acronym 
of  the  words  “north,"  “east,"  “west."  “south,"  "up,"  and 


f'lg  4  Confjguratii»n  for  linear  interp<>laiion  in  ihe  NliWSl’D 
mcthixi 


■'down"  since  the  stencil  contains  only  one  cell  near  each 
of  the  six  faces  of  a  central  cube  for  calculations  within 
that  cube. 

When  the  NEiW.Sl  D  interpolant  is  used  for  calcula¬ 
tions  within  a  cell,  first,  one  imaginary  plane  is  passed 
through  each  edge  of  the  cube  and  the  veil  centroid.  The 
planes  divide  the  cube  into  six  congrucm  pyramids  as  well 
as  dividing  the  exterioi  into  six  corresponding  solid  an¬ 
gles.  Next,  each  solid  angle  is  searched  to  find  a  cell 
centroid  as  close  as  possible  to  the  interpolated  cube. 
Interpolation  is  performed  sep.irately  in  each  of  the  six 
pyramids  using  only  the  centroids  of  the  interpolated  cube 
and  corresponding  external  cell.  If  a  solid  angle  does  not 
contain  an  external  cell  centroid  within  a  reasonable  dis¬ 
tance,  say  one  side  of  the  central  cell  or  less,  the  uninter¬ 
polated  value  of  E  is  used  within  the  corresponding 
pyramid. 

Ihe  configuration  for  linear  mterpohition  within  ;i 
pyrtimid  in  the  Nl.W.Sl  1)  methorl  is  shr'wn  in  I  ig.  4.  The 
point  at  which  interpolation  is  desired  (  \’.  >.  / 1.  the 
centroid  of  the  external  cell  (  A,.  )  ,.  /, ).  and  the  centroid 
v>f  the  interpolated  cell  which  is  used  as  the  (vrigin  de¬ 
termine  a  plane  which  is  used  for  the  figure. 

The  equation  of  the  plane  perpendicular  to  the  line 
connecting  (.V,.  Y/.  7.,)  to  the  origin  and  passing  through 
the  point  (,V.  T,  Z)  is 

.V-V,-l-  YY,  +  //.^  =  ayl\}+  T,’-i-Zr  (7) 

where  o  is  the  distance  from  the  origin  to  the  plane.  The 
linear  interpolation  used  in  the  NEWSUD  methrxl  when 
an  externa!  cell  is  found  corresponding  to  a  pyramid  of 
the  interpolated  cell  is  accomplished  by  the  approxima¬ 
tion  of  identifying  E(X,  Y.  Z)  with  t'f.V',  T  ,  Z  )  so  that 


}IxJ+  ) ;  +  Z  |- 


'  r  +  ^  r  + 


E(O.O.O).  (8) 


l.inear  interpolation  defined  by  (7)  and  (8)  is  used  to 
evaluate  the  integral  of  1/2  a  E-E*  in  a  pyramid  for 
which  a  corresponding  external  cell  is  found.  If  the  cells 
have  different  complex  permittivities,  (8)  must  be  nuxi- 
ified  by  suitably  multiplying  by  the  ratio  of  permittivities 
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Fig  5.  An  impro'  ed  model  of  man  for  numerical  calculations. 


Fig.  6.  Whole-bixly  SAR  for  a  homogeneous  model  of  man.  k 

fronl-to-back.  incident  intensity  =  I  MW/cm^ 


SO  that  both  D  normal  and  E  tangential  to  a  boundary  are 
continuous. 

Fig.  5  illustrates  a  model  of  man  which  has  been  u.sed 
for  numertcal  calculations  [8|.  Moment-methttd  solutions 
have  been  made  for  the  model  using  pulse  functions  as  a 
basis  and  delta  functions  for  testing.  The  NEWSUD  inter- 
polant  has  been  used  in  calculations  since  the  model  has 
cells  which  are  cubes  of  different  sizes.  Fig.  6  shows  the 
frequency  dependence  of  the  average  SAR  calculated  for 
the  model  of  man  with  a'  vertically  polariz.ed  plane  wave 
having  frontal  incidence.  Values  found  both  with  and 
without  interpolation  are  given. 

If  the  local  E  values  are  exact  samfdes  and  if  there  is 
appreciable  variation  between  adjacent  cells,  then  statisti¬ 
cal  procedures  may  be  used  to  calculate  confidence  limits 
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Fig.  7.  Correction  to  whole-body  SAR  for  homogeneous  model  of 

man 


for  values  of  the  average  SAR  found  without  aid  of 
interpolants  [9].  The  80-percent  confidence  half  widths  for 
average  SAR  values  have  been  calculated  and  plotted  in 
Fig.  7.  The  percent  correction  which  the  NEWSUD  inter- 
polant  makes  is  also  shown  in  Fig.  7. 

The  frequency-dependent  vanation  of  the  true  E  must 
cause  increasing  error  in  calculations  of  local  E  at 
frequencies  greater  than  200  MHz  for  the  cell  sizes  used 
with  the  model  of  man  [10].  The  80-percent  confidence 
half  widths  shown  in  Fig.  7  are  not  valid  at  frequencies 
above  200  MHz  since  their  calculation  requires  the 
assumption  that  the  local  E  values  are  exact  samples.  Fig. 
7  shows  that  the  NEWSUD  interpolant  makes  relatively 
little  correction  to  the  SAR  at  low  frequencies  where 
relatively  little  error  is  expected,  but  the  correction  in¬ 
creases  sharply  at  frequencies  for  which  the  error  is  ex¬ 
pected  to  increase. 


V.  CONCIISIONS 

The  interpolative  correction  to  calculated  values  of  the 
average  SAR  has  been  tested  in  a  one-dimenstonal  prob¬ 
lem  in  which  the  analytical  solution  is  available  for  com¬ 
parison  with  numerical  solutions.  A  test  was  also  made  in 
a  three-dimensional  problem  in  which  the  true  solution 
was  estimated  by  extrapolation.  In  both  tests  the  inter¬ 
polant  makes  small  corrections  when  few  cells  are  used 
since  sampling  is  insufficient  to  allow  proper  estimation  of 
the  variation  of  the  fields.  The  fractional  correction  is  also 
observed  to  decrease  as  many  cells  are  used  since  the 
interpolated  and  uninterpolated  solutions  must  converge 
to  the  same  answer.  In  both  tests  the  interpolant  was 
found  to  provide  a  significant  reduction  in  error  for  a 
negligible  increa.se  in  computational  expense. 

It  is  hoped  that  the  two  interpolants  presented  in  this 
paper  are  only  the  first  steps  in  the  development  of  a 
procedure  which  will  find  general  usage  as  a  follow-up  to 
moment-method  solutions.  Use  of  the  interpolant  for 
calculation  of  the  external  scattered  fields  and  develop¬ 
ment  of  a  general  interpolant  consistent  with  the  wave 
equation  appear  to  be  the  next  necessary  steps. 
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